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It  MawKi 

Since  1977,  the  U.S.  Coast  Guard  has  been  conducting  studies  of  the  suitability 
of  Loran-C  as  a  precision  aid  to  navigation  in  the  harbor-harbor  entrance  (HHE) 
areas  of  the  continental  U.S.  The  final  phase  of  this  effort  involves  an 
assessment  of  the  stability  of  the  signals  of  the  existing  Loran-C  system  along 
with  an  examination  of  stability  improvement  methods.  The  final  efforts  were 
begin  early  in  1981  with  the  deployment  of  lor  an  data  collection  sets  (the 
so-called  "harbor  monitors")  in  select  harbor  areas.  In  cooperation  with  the 
St.  Lawrence  Seaway  Development  Corporation,  the  St.  Lawrence  Seaway  Authority 
and  the  Canadian  Coast  Guard,  these  efforts  were  extended  to  the  St.  Lawrence 
Seaway  area  in  late  1981.  Resulting  data  is  analyzed  in  the  report  and 
indicates  the  performance  of  the  existing  Loran-C  system  is  inadequate  to 
support  precisian  applications  in  the  St.  Lawrence.  The  performance  of 
Differential  Loran-C  is  hypothesized  and  found  to  be,  at  best,  marginal  due  to 
suboptimal  system  geometry.  The  performance  which  would  result  with  the 
addition  of  another  transmitting  station  is  also  hypothesised.  With  this 
improved  geometry,  "raw”  Loran-C  performance  is  still  found  inadequate  but 
Differential  Loran-C,  yielding  maximum  cross-track  errors  of  27  meters  (99.9% 
probability),  leaves  23  meters  error  margin  for  the  largest  vessels. 
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EXECUTIVE  SUM4ARY 


In  mid-1981,  the  St.  Lawrence  Seaway  Development  Corporation.  St. 
Lawrenoe  Seaway  Authority,  and  Canadian  and  U.S.  Coast  Guards  began  a  joint 
effort  to  explore  the  capability  of  Loran-C  to  meet  precision,  all-weather 
navigation  requirements  in  the  Seaway.  Loran-C  surveys  of  major  portions  of 
the  Seaway  have  been  made  and  guidance  equipment  demonstrations  are  planned 
for  late  summer  1982  and  beyond. 

Critical  questions  of  concern  relate  to  the  stability  of  the  existing 
Loran-C  grid  in  the  Seaway  and  the  feasibility  of  methods  to  improve  the 
stability.  To  begin  to  address  these  questions,  Loran-C  data  collection 
equipment  sets  were  deployed  at  several  sites  along  the  Seaway  in  late  1981. 
The  U.S.  Coast  Guard,  as  part  of  an  ongoing  R&D  project,  had  'eveloped 
equipment,  the  so-called  Harbor  Monitor  Set,  which  was  ideally  e  ed  to  a 
study  of  this  nature.  Unfortunately,  most  setB  were  already  deploy  at  sites 
throughout  the  U.S.  Originally,  therefore,  only  one  Harbor  Monit  Get  was 
deployed  in  the  Seaway  -  at  the  Eisenhower  Lock,  Massena,  N.Y.,  October 
1981. 

The  Loran-C  receivers  to  be  used  with  future  Harbor  Monitor  Sets  were 
available  for  use  so,  along  with  some  teleprinters/cassette  data  recorders, 
these  were  made  into  so-called  "Ad  Hoc"  equipment  sets  and  deployed  at  three 
additional  sites  along  the  Seaway  in  December  1981.  Additionally,  the 
Canadian  Coast  Guard  was  in  the  process  of  developing  Loran-C  data  collection 
equipment  for  use  in  a  buoy-auditing  project.  Development  was  accelerated  and 
these  sets  were  deployed  at  two  sites  along  the  Seaway  late  in  the  period  of 
this  report.  In  early  February  1982,  USCG  RiD  Center  personnel  installed 
another  Harbor  Monitor  Set  at  Tibbetts  Point  N.Y.  A  preliminary  analysis  of 
the  data  collected  with  these  equipment  sets  is  provided  in  this  report. 

Data  obtained  by  the  Harbor  Monitor  Sets  has  been  of  excellent  quality  - 
far  exceeding  the  quality  obtained  in  the  St.  Marys  River  Stability  Study  for 
which  earlier  versions  of  the  equipment  were  developed.  Conversely,  except 
for  a  critically  important  month  long  period  during  which  the  Canadian  Coast 
Guard  was  able  to  obtain  data  at  Iroquois  Lock,  continual  problems  have  been 
encountered  trying  to  get  reliable  data  from  the  other  sites.  Hie  key  feature 
of  the  Harbor  Monitor  Set  that  leads  to  the  higher  reliability  is  the  remote 
control  capability.  via  phone  line  access,  project  personnel  at  the  RtD 
Center  are  able  to  closely  monitor  the  status  of  the  equipment  on  a  daily 
basis.  Problems  which  threaten  to  cause  loss  of  data,  primarily  a  result  of 
electrical  power  interruption,  can  be  detected  and  corrected  in  a  timely 
fashion.  Unfortunately,  similar  problems  are  not  detected  at  the  other  sites 
until  considerable  data  has  been  lost.  Thus,  a  major  conclusion  of  the  study 
to  date  is  that  Harbor  Monitor  Sets  should  be  deployed  at  all  sites. 

A  major  stumbling  block  in  the  execution  of  a  study  such  as  this  is 
identified  at  the  beginning  of  the  report:  there  are  no  standard  guidelines 
for  determing  navigation  system  requirements  for  restricted  waterways.  Thus, 
the  recommendations  of  a  report  describing  navigation  system  accuracy,  or  ways 
to  improve  it,  lack  a  foundation.  The  report  acknowledges  this  problem  and 
offers  an  assumed  set  of  requirements.  Hie  assumptions  are  based  on  similar 
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experiences  to  date  and  the  results  of  ongoing  "requirements"  studies. 
Wherever  possible,  care  has  been  taken  to  tabulate  the  results  of  the 
stability  study  in  a  method  amenable  to  update  when  more  concrete  requirements 
can  be  stated. 

Prior  to  considering  any  of  the  data  obtained  along  the  Seaway,  the 
report  presents  an  extensive  series  of  predictions,  ltie  predictions  are  based 
on  a  model  of  the  nature  of  Loran-C  time  difference  variations  which  has  been 
developed  over  the  years  along  with  an  extrapolation  of  data  obtained  at 
Harbor  Monitor  sites  located  along  the  northeast  U.S.  coast.  The  methodology 
involved  in  the  prediction  process  provides  a  framework  in  which  the  data  can 
be  viewed  and  allows  an  extensive  set  of  conclusions  to  be  drawn  from  the 
limited  data  base. 

The  predictions  indicate  trie  stability  of  "raw"  (unaugmented)  Loran-C  is 
inadequate  for  precision  use  throughout  the  shipping  season  in  the  Seaway. 
Data  from  Massena  and  Tibbetts  Point  confirm  this.  Extending  the  use  of  the 
prediction  model,  Differential  Loran-C  performance,  using  shore  stations  at 
Massena  and  Wellesley  Island,  is  hypothesized.  Adequate  performance  is 
expected  except  for  several  reaches  in  the  region  from  the  Brockville  Narrows 
to  the  Thousand  Island  Bridge.  Here,  adequate  performance  is  defined  to  mean 
the  99.9%  probability  cross-track  error  estimate  leaves  at  leaBt  10  meters  for 
guidance  error  for  a  vessel  with  a  beam  of  34  meters.  Further  investigation 
shows  the  major  problem  is  the  suboptimal  geometry  associated  with  the  loran 
grid  in  the  Seaway.  Fixes  obtained  by  Loran-C  over  the  course  of  a  year,  or 
just  the  shipping  season,  will  form  an  elliptical  pattern  at  any  point  of 
measurement.  In  the  Seaway,  the  longer  axis  of  the  ellipse  is  nearly 
perpendicular  to  the  predominant  course.  The  geometry  is  bad  enough  to  mean 
that  if  the  prediction  regarding  the  magnitude  of  the  short  term  variations 
were  either  true  or  optimistic.  Differential  Loran-C  could  n^.t  produce 
satisfactory  performance  in  these  reaches. 

In  spite  of  the  sub-par  data  base,  an  elaborate  series  of  analysis 
techniques  allows  us  to  determine  that  the  predictions  regarding  the  short 
term  variations  were  pessimistic.  In  the  process,  we  discover  that  the 
predictions  regarding  the  magnitude  of  the  long  term  variations  were 
optimistic.  The  resulting  conclusion,  although  more  data  is  needed  to 
"harden"  the  estimates,  is  that  Differential  Loran-C  appears  capable  of  the 
desired  performance  -  barely- 

This  having  been  said,  vc  caution  that  Differential  Loran-C  is  merely  a 
concept  -  never  having  been  implemented.  Results  of  this  study  to  date 
indicate  a  flawlessly  operated  Differential  Loran-C  system  would  produce 
satisfactory  performance  with  a  moderate  amount  of  room  left  for  imperfect 
surveying  and  other  practical  limitations.  Again,  the  magnifying  effect  of 
poor  loran  geometry  in  the  source  of  the  problem.  Differential  Loran-C  tests 
are  encouraged  so  that  an  appreciation  for  the  "other  practical  limitations" 
can  be  gained.  Pesulting  insights  will  allow  the  statement  of  the  "final 
word"  on  the  subject.  Tbe  concept  to  bear  in  mind  is  that  we  are 
contemplating  the  first  operational  use  of  Differential  Loran-C  in  an  area 
where  there  is  scant  room  for  error. 


As  a  final  note,  and  as  a  raault  of  tha  conclusions  regarding  poor 
existing  geoawtry*  we  introduce  a  consideration  of  what  would  result  were 
another  Loran-C  transmitting  station  added  to  the  existing  chain.  With  a 
station  at  North  Bay*  Ontario*  predictions  show  "raw"  Loran-C  still  will  not 
produce  satisfactory  performance  throughout  the  shipping  season.  The  use  of  a 
single  differential  monitor  station*  in  conjunction  with  this  new  transmitting 
station*  however*  allows  performance  which  easily  exceeds  all  requirements. 
This  suggests  the  addition  of  the  new  transmitting  station  would  allow  enough 
room  for  error  that  a  "milder"  alternative  (e.g.*  dally  corrections)  to  full 
Differential  Loran-C  could  be  considered.  The  alternatives  which  require  a 
new  transmitting  station  Involve  considerable  expense  and,  therefore,  risk. 
At  present,  they  are  offered  as  subjects  for  careful  scrutiny  in  future 
studies. 


1. 


INTRODUCTION 


1-1  USOG  Harbor/Harbor  Entrance  Loran-C  R6D  Studies. 


Since  1977  the  U.S.  Coast  Guard  has  been  sponsoring  Research  and 
Development  studies  to  assess  the  capability  of  the  Loran-C  radionavigation 
system  to  support  precision  navigation  requirements  typically  encountered  in 
piloting  situations.  TTie  studies  have  centered  on  four  areas  of 
investigation : 

a.  Measuring  the  year-round  stability  of  the  loran  signals  in 
selected  harbor/harbor  entrance  ( HHE)  areas. 

b.  Investigating  methods  of  enhancing  system  performance  for 
areas  in  which  system  stability/geometry  is  not  adequate  for  high  precision 
use . 

c.  Developing  survey  techniques  to  precisely  determine  the  loran 
coordinates  at  critical  "waypoints"  in  channels. 

d.  Developing  modern  guidance  equipment  which  displays  loran 
information  in  a  manner  most  useful  to  a  mariner  in  piloting  situations. 

Two  guidance  equipment  systems/ developed  for  the  USCG  by  the  Johns 
Hopkins  Applied  Physics  Lab  (APL)  under  this  program,  have  proven  to  be 
effective  in  exploiting  the  high  precision  of  Loran-C  in  stable  signal 
areas.  One  such  system,  called  PILOT  (Precision  JUitracoastal  iDran 
Translocator),  features  a  graphical  display  of  the  HHE  area  being  transitted 
and  is  intended  for  permanent  installation  on  a  vessel.  The  other  system, 
called  PIAD  (Portable  Loran  Assist  Device)  features  an  alpha-numeric  display 
of  ship's  position  relative  to  channel  waypoints,  and  is  designed  to  be 
carried  aboard  a  vessel  by  pilots.  Both  devices  are  described  in  more 
detail  in  reference  1. 


1.2  St  Lawrence  Seaway  Development  Corporation  Interests. 


In  February  1981,  the  Office  of  Comprehensive  Planning  of  the  Saint 
Lawrence  Seaway  Development  Corporation  (SLSDC)  asked  for  an  APL/USCG 
simulator  demonstration  of  PILOT.  As  a  result  of  the  simulator 
demonstration,  the  SLSDC  concluded  Loran-C,  perhaps  only  as  a  part  of  an 
integrated  navigation  package,  had  the  potential  to  satisfy  navigation 
requirements  on  the  St.  Lawrence  Seaway.  Thus,  SLSDC  asked  for  USCG 
assistance  in  sponsoring  a  PILOT /PLAD  demonstration  on  the  St.  Lawrence 
Seaway. 

In  August  1981,  SLSDC  and  USCG  personnel  conducted  a  trackline  survey 
of  the  St.  Lawrence  Seaway  from  the  vicinity  of  Ogdenburg,  N.Y.  to  Cornwall, 
Ontario.  PILOT  chartlets  were  prepared  for  the  surveyed  area  and  a 
verification  run  was  made  in  late  October  1981. 
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Meanwhile,  with  a  reasonable  assurance  that  the  survey  and  PILOT 
chartlet  preparation  efforts  would  be  a  success,  agency  heads  from  the 
SLSDC,  the  Saint  Lawrence  Seaway  Authority  and  the  Canadian  Coast  Guard 
(CCG)  met  with  the  Commandant  of  the  USCG  to  discuss  further  plans.  It  was 
determined  that  a  stability  study  would  be  initiated  and  that 
representatives  from  all  four  agencies  would  form  a  "Precision  All-Weather 
Navigation  System"  (PAWNS)  Steering  Committee.  The  committee  would 
initially  be  concerned  with  insuring  the  lor an  studies  were  underway  and 
then  proceed  to  study  matters  such  as  precision  requirements  and 
requirements  for  an  integrated  systems  package. 

By  late  October  1981,  USCG  and  CCG  personnel  began  deploying  several 
types  of  loran  data  collection  units  along  the  St  Lawrence  Seaway.  It 
should  be  noted  that  these  diverse  units  were  deployed  on  a  hurried  basis  so 
that  conditions  during  the  end  of  the  shipping  season  could  be  monitored. 

At  a  subsequent  meeting  of  the  Steering  Committee  it  was  determined  that  an 
early,  initial  review  of  the  data  should  be  made  so  that  any  equipment  or 
data  compatability  problems  could  be  identified  and  rectified  before  the 
next  critical  period  of  the  season  (i.e.,  the  Fall  of  1982)  was  entered. 

This  report  summarizes  data  collection  results  to  date  and  addresses  the 
problem  areas  of  concern  to  the  Steering  Committee. 

1.3  Applicability  of  Harbor  Monitor  Program  to  Seaway 


It  should  be  emphasized  at  the  start  that  there  is  considerable  concern 
as  to  the  suitability  of  the  USCG  Harbor  Monitor  Program  for  answering  all 
of  the  questions  of  concern  in  the  St.  Lawrence  Seaway.  The  USCG  Stability 
Studies  are  directed  at  providing  a  general  characterization  of  Loran-C 
accuracy  throughout  all  coastal  and  harbor  areas  of  the  continental  U.S. 

With  such  a  broad  objective  in  mind,  the  experiment  strategy  recognizes  it 
is  impossible  to  obtain  year-round  data  in  every  area  of  interest.  This  is 
not  just  a  result  of  inability  to  procure  enough  equipment  -  the 
availability  of  skilled  technicians  to  support  the  equipment  is  considerably 
more  of  a  limitation.  Thus  the  program  attempts  to  follow  a  strategy  which 
features  spatial  sampling  adequate  to  allow  interpolation  of  the  results  to 
all  areas  of  interest. 

Additionally,  the  program  features  time  sampling,  i.e.,  data  is  not 
collected  on  a  continuous  basis  throughout  the  year.  A  key  consideration 
here  is  that  whereas  recent  advances  in  computer  technology  h«ve  made  it 
possible  to  amass  and  maintain  huge  data  bases  at  reasonable  costs,  the 
analysis  will  only  be  as  "believable"  as  the  data.  To  maintain  data  base 
integrity,  unfortunately,  requires  skilled  technical  review  of  the  data 
being  accepted  as  valid  (e.g.,  not  produced  by  faulty  instrumentation). 

Thus,  in  spite  of  the  computer  revolution,  the  cost  of  the  data  base 
remains,  essentially,  directly  proportional  to  size  and  a  reasonable 
quantity /quality  tradeoff  must  be  sought. 

Drawing  on  the  results  of  numerous  experiments  conducted  over  the  years 
plus  the  benefits  of  many  years  of  Loran-C  system  operational  observations, 
the  USCG  adopted  the  two  hour  per  day  sampling  strategy  as  the  tradeoff 
commensurate  with  the  harbor  monitor  effort.  Whereas  there  is  agreement  on 


the  Steering  Committee  that  tha  strategy  la  a  correct  on#  for  the  general 
Loran-C  characterisation,  there  la  concern  it  nay  not  be  adequate  for 
special,  in-depth  applications  such  aa  the  Seaway.  Thus,  considerable  care 
will  be  given  in  t/iis  report,  mtd  in  further  work,  to  discuss  and  teat  the 
strategy. 


1.4  Treatment  of  Systaa  Mequireaenta  in  This  Report 

Ideally,  the  goal  of  a  atudy  such  as  this  should  extend  beyond  siaply 
stating  what  the  stability  of  the  Loran-C  systaa  is  or  could  be  aade  to  be. 
It  should  assess  the  existing  or  potential  stability  with  respect  to 
navigation  systaa  accuracy  retail raaenta.  At  present,  however,  there  are  no 
standard  aethodo logics  for  quantifying  accuracy  requirements  in  restricted 
waterways.  The  consequence  is  that  although  systaa  s tab ility/accu racy  can 
be  described  in  highly  quantified  terns,  the  ultiaate  question  "is  it  good 
enough?"  aust  typically  be  avoided. 

Reference  11  documents  results  of  ongoing  D.S.  Coast  Qiard  sponsored 
studies  to  address  this  "requlraaMnts"  problem.  Although  results  to-date 
are  encouraging,  it  aust  be  conceded  that  the  questions  are  non-trivial  and 
prospects  for  definitive  answers  in  the  near  future  are  not  good.  It  should 
also  be  emphasised  that  this  is  not  a  "Loran-C  per  se"  problea:  the  same 
unanswered  questions  apply  to  all  precisian  navigation  systems  be  they 
conventional  visual  or  radionavigation  systems.  A  aajor  task  of  the 
Steering  Committee  will  be  to  provide  some  answers  to  the  "requirements” 
question. 

Having  recognised  that  system  requireamnts  cannot  yet  be  conclusively 
stated,  the  report  also  recognizes  a  statement  of  existing  stability,  or 
ways  to  improve  stability,  without  reference  to  any  requirement,  is  of 
marginal  utility.  Conversely,  the  report  recognizes  a  tutorial  on  the 
nature  of  Loran-C  instabilities  and  "enhancement"  techniques  can  be  of 
significant  value.  A  meaningful  statement  of  requirements  probably  cannot 
be  made  without  an  understanding  of  such  issues.  Thus,  the  report  proceeds 
with  an  assumed  set  of  requirements.  Since  these  assumed  requirements  may 
someday  prove  to  be  in  need  of  revision,  considerable  effort  is  made  to 
present  the  results  in  a  fora  amenable  to  future  update. 

In  developing  the  assumed  requirements,  past  experiences  are  drawn 
upon.  As  indicated  in  reference  4,  the  originally  stated  requirement  for 
the  St.  Marys  River  navigation  systea  was  10  feet.  When  it  became  apparent 
that  no  practical  systea  could  provide  such  accuracy  (but  that  it  probably 
was  not  required),  Loran-C  was  offered  as  a  "100  foot  systea.”  This  report 
demonstrates  why  requirements  statements  of  this  nature  are  to  be  avoided. 

In  deference  to  the  people  involved  in  drawing  up  the  St.  Marys  River 
requirements  statements,  they  actually  meant  to  say  that  10  feet  (later  100 
feet)  was  the  requirement  in  the  "tightest  reaches."  A  key  concept  for 
anyone  who  would  state  requirements  to  recall  is  that,  at  least  for  fixed 
location  earth-bound  systems,  both  measureaent  accuracy  and  the  "geometry 
factors”  which  transform  the  measurement  errors  to  position  errors  vary  from 
point  to  point  in  the  service  area,  (for  satellite  or  celestial  systems, 
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the  geometry  also  varies  froa  tiae  to  tlae.)  A  system  is  best  designed  when 
this  varying  aocuracy  is  made  to  be  the  best  in  the  tightest  channels.  In 
the  case  of  the  St.  Marys  River,  the  100  foot  accuracy  is  only  required  in  a 
region  ccaprising  about  10  %  of  the  river  length.  Over-specification  by  an 
order  of  magnitude ,  as  was  done  several  times  in  the  St.  Marys  River,  will 
typically  render  a  system  *non-cost-benef icial.* 

Another  important  consideration  demonstrated  in  this  report  is  that 
aocuracy  is  a  two-dimensional  quantity  for  marine  applications.  As  an 
example,  we  can  imagine  a  situation  in  which  the  true  requirement  is  only 
that  cross-track  error  be  less  than  100  feet.  If  this  is  not  explicitly 
stated,  m  could  inoorrectly  declare  inadequate  a  system  which  features  a  SO 
foot  cross-track  aocuracy  in  the  area  of  concern  because  the  total 
(root -sum-squared  of  cross-track  and  along-track)  error  exceeds  100  feet. 

With  this  background  as  to  motives,  the  report  makes  the  following 
assumptions  about  St.  Lawrence  8eaway  requirements: 

(1)  Because  of  the  orientation  of  the  Lor an -C  error  ellipses  and  the 
predominant  courses  of  the  Seaway  (to  be  describe),  cross-track  error  is 
the  key  parameter.  Along-track  errors  (indeed,  actual  error  ellipses)  for 
all  reaches  are  furnished  in  the  Appendices  in  case  future  work  indicates 
the  need  for  examination  of  these  quantities. 

(2)  Cross-track  error  requirements  are  assumed  to  be  direct  functions 
of  channel  width  and  vessel  width.  Specifically,  the  total  allowable 
cross-track  error  is  assumed  to  be  defined  as  one-half  the  channel  width 
minus  one-half  the  beast  of  the  largest  vessel  using  the  Seaway.  This 
relationship  causes  the  requirements  to  be  sore  stringent  in  the  narrower 
reaches  -  as  desired. 

(3)  Of  the  total  allowable  error,  some  portion  must  be  allocated  to 
the  navigation  system  and  soam  to  guidance  error  (i.e.,  trackkeeping 
ability).  The  Loran-C  errors  are  specified  at  the  99.9%  probability  level. 
The  95%  probability  levels  are  also  provided  as  an  example  of  how  to  use  the 
mathematics  of  Appendix  A  to  extend  the  analysis  to  any  probability. 

(4)  Max  liman  guidance  errors  of  10  meters  are  assumed  as  the  basis  for 
requiring  Loran-C  chain  enhancements.  As  with  the  loran,  there  is  no  reason 
to  expect  the  guidance  error  should  be  independent  of  the  reach 
characteristics  but  an  assumption  of  some  sort  is  required.  The  10  meters 
represents  a  stringent  demand  on  trackkeeping  skills  but  is  consistent  with 
the  stringent  99.9%  demands  being  placed  on  the  Loran-C  and  is  assumed 
consistent  with  normal  demands  imposed  by  the  narrowness  of  the  channels  for 
which  it  is  an  issue.  Bctensive  use  of  tables  is  made  so  that  the  results 
can  be  extended  to  whatever  guidance  errors  future  studies  may  indicate  to 
be  appropriate. 


1.5  Report  Overview. 


Section  2  begins  with  a  description  of  the  various  equipment  sets  used 


to  collect  data  thus  far.  Next,  a  ay nope la  of  data  collection  efforts  to 
date  is  presented.  Section  2  concludes  with  the  presentation  of  a  newly 
developed  data  collection  set  and  a  discussion  of  compatabllity  for  all  of 

the  sets. 

Section  3  begins  the  discussion  of  the  adequacy  of  the  data  collection 
strategy  by  Introducing  a  nodal  for  Loran-C  variations.  The  need  for  the 
model,  and  how  it  relates  to  both  the  systeai  requirements  and  the  scope  of 
the  experiment  is  developed.  Due  to  late  starts  and/or  site  and  equipment 
problems,  we  do  not  yet  have  a  sufficient  data  base  to  say  the  final  word  on 
the  adequacy  of  Loran-C  in  the  Seaway.  We  can,  however,  use  the  model  to 
make  some  preliminary  predictions.  We  draw  upon  data  obtained  from 
Northeast  U.S.  chain  operational  and  harbor  monitor  sites  to  make  the 
predictions  which  oonclude  Section  3.  The  adequacy  of  these  predictions  can 
be  tested  somewhat  with  data  collected  thus  far  in  the  Seaway  and  with  data 
collected  from  future  efforts.  The  model  will  probably  need  future 
refinement,  but  can  serve  to  steer  future  efforts  in  the  right  direction. 

In  Section  4  we  present  the  results  of  the  data  collection  efforts. 

The  results  are  encouraging,  both  in  regards  to  the  stability  they  indicate 
and  the  agreement  with  the  predictions.  What  is  not  encouraging,  however, 
is  the  problem  we  have  experienced  in  obtaining  usable  data.  Of  particular 
concern  is  the  inadequacy  of  the  so-called  "Ad  Roc"  data  sites  which,  it  is 
shown,  threatens  to  compromise  future  efforts.  We  confirm  the  preliminary 
conclusion  of  Section  3  that  an  additional  monitor  site  may  be  needed. 

Section  5  provides  scam  brief  comments  regarding  the  implementation  of 
differential  Loran-C  -  a  non-trivlal  undertaking.  An  alternative  is 
introduced  for  initial  consideration. 

Section  6  summarises  the  conclusions  -  the  prime  one  being  that  data 
collection  efforts/equipment  must  be  upgraded  before  this  fall.  Conclusions 
about  the  sampling  strategy  are  also  presented.  Finally,  the  significant 
effect  that  existing  chain  geometry  has  on  the  error  budget  and  the 
viability  of  Loran-C  as  a  precise  navald  for  the  Seaway  is  discussed.  The 
report  ends  with  a  tabulation  of  the  the  recommendations  that  seem  indicated 
from  the  results  to  date.  In  short,  further  efforts  with  upgraded  equipment 
are  recommended  for  the  stability  study  and  a  test  implementation  of 
Differential  Loran-C  is  urged. 
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2. 


STABILITY  STUDY  EQUIPMENT /SITES 


2.1  USOS  Harbor  Monitor  Seta 


As  discussed  in  paragraph  1.1,  the  08CG  has  developed  a  so-called 
"Harbor  Monitor  Set”  for  use  in  the  HHE  Loran-C  Stability  Study  project 
eleaent.  The  set  consists  of  a  survey  grade  Loran-C  receiver,  the  Internav 
404,  a  snail  micro-computer,  and  alscellaneous  support  equipment  allowing 
battery  back-up  and  resete  (via  telephone)  access  to  the  collected  data. 

The  equipment  is  designed  to  accomodate  a  "low  density  data  analysis”  of 
lor an  signals  as  described  in  reference  2.  Simply  stated,  this  means  data 
is  collected  during  only  two  hours  of  the  day  -  generally  from  11  to  12,  AM 
and  PM.  Figure  2-1  shows  the  prime  components  of  the  system  in  block 
diagrasi  form. 


Figure  2-1  Harbor  Monitor  Set  Block  Diagram 


Simulator  tests  have  shown  the  Loran-C  receiver  has  a  "servo  loop  time 
constant”  of  about  6-8  seconds  under  conditions  typically  encountered  at 
harbor  monitor  sites.  Thus  the  receiver  output  is  sampled  every  40  seconds 
so  that  the  samples  can  be  treated  as  statistically  independent.  The 
micro-computer  uses  a  real-time  clock  to  begin  the  sampling  period  at  the 
prescribed  time.  At  the  end  of  the  sampling  period,  the  mean,  standard 
deviation,  and  minimum  and  amximum  values  of  each  time  difference  are 


recorded.  Depending  upon  the  number  of  signals  being  observed,  available 
memory  will  hold  from  10  to  20  days  of  data.  Phone  line  access  to  the 
micro-computer  allows  retrieval  of  the  stored  data.  It  also  allows  a 
remotely  located  operator  to  proaipt  the  computer  to  exercise  any  receiver 
command  which  a  local  operator  could  enter  via  the  front  panel  controls. 
Finally,  the  entire  micro -computer  program  can  be  changed  via  the  phone 
line.  The  concern  about  the  adequacy  of  this  type  of  data  for  the  study  at 
hand  was  introduced  in  Section  1.3.  Reference  2  addresses  the  "adequacy” 
question  to  a  certain  extent  and  more  comments  will  be  provided  in  the  body 
of  this  repot  *  . 

Regardless  of  the  question  of  "inadequacy,"  a  Harbor  Monitor  Set  was 
installed  at  Eisenhower  Lock,  Massena,  N.Y.  on  21  October  19S1  because  it 
met  one  vital  criterion:  it  was  immediately  available  for  deployment.  Data 
suitable  for  low  density  analysis  has  been  obtained  from  the  Massena  Harbor 
Monitor  set  from  October  1991  to  date  and  is  presented  in  subsequent 
sections  of  this  report. 

On  1  February  1982,  an  additional  Harbor  Monitor  Set  was  installed  on 
USCG  property  at  Tibbetts  Point  (Cape  Vincent)  N.Y.  Subsequently,  the 
Harbor  Monitor  Set  was  relocated  to  a  nearby  facility  in  the  area.  Data 
from  the  site(s)  from  2  February  1982  to  date  is  available  and  is  discussed 
in  this  report.  Figure  2-2  shows  the  harbor  monitor  set  installed  at 
Massena . 


Figure  2-2 

t 


Harbor  Monitor  Set  at  Massena,  N.Y. 
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2.2  COG  Lor an -C  Buoy  Auditing  Bquipment 


The  Canadian  Coaat  Guard  had  baan  in  tha  process  of  davaloping  protoypa 
equipment  aata  to  support  a  progrmi  of  auditing  tha  position  of  buoys  via 
Lor  an -C.  The  equipment  is  oemprised  of  an  Zntarnav  404  Lor an -C  racaivar 
(tha  saaa  racaivar  usad  with  tha  Harbor  Monitor  8ats) ,  a  Motorola 
aicroprocassar  and  a  Texas  Instruments  talapr inter /data  recorder.  By  using 
data  recording  cartridges,  tha  sets  are  capable  of  storing  significantly 
sore  data  than  tha  Harbor  Monitors  -  thereby  allowing  a  true  "high  density" 
analysis  of  tha  lor an  data.  Figure  2-3  shows  the  prime  components  of  the 
system  in  block  diagram  form. 


Figure  2-3  Buoy  Auditing  Set  Block  Diagram 


1*o  receivers  are  usad  because  of  a  buoy  auditing  project  requirement. 
The  micro-computer  obtains  smaplaa  of  tha  data  from  tha  receivers 
approximately  once  every  2  seconds.  After  100  samples  are  obtained,  the 
micro-computer  sends  the  computed  mean  and  standard  deviation  of  the  time 
differences  to  the  cassette  recorder.  Such  data  is  recorded  throughout  the 
day  mid  the  cassettes  must  be  changed  about  every  fourth  day. 

The  only  major  drawback  of  this  equipment  set  is  that  it  cannot  be 
remotely  accessed.  Consequently,  monitoring  of  the  equipment  status  by 
experienced  project  peraonnel  is  not  possible  on  a  day-to-day  basis  (nor  is 
it  necessary  -  for  the  buoy  auditing  project).  An  additional  drawback  is 
that  the  wilts  were  still  in  the  development  stage  when  they  had  to  be 


pressed  into  service  for  the  purposes  of  this  study.  Thus,  this  is  actually 
the  first  true  field  test  and  evaluation  of  the  units.  The  first  set  was 
deployed  at  Fort  Wellington  on  9  February  1982.  Due  to  electrical  grounding 
problems  at  the  sites,  no  useful  data  was  obtained.  The  site  was  moved 
several  hundred  yards  to  where  a  good  equipment  ground  could  be  established 
on  11  March  1982  and  good  results  were  achieved.  Soon,  however,  power  and 
equipment  problems,  causing  loss  of  data  from  2  to  19  April  1982,  were 
encountered.  On  20  April  1982,  the  equipment  was  relocated  to  the  nearby 
Prescott  CCG  Base.  A  second  set  was  installed  at  Iroquois  lock,  Ont.  on  21 
April  1982.  Although  data  suitable  for  "high  density"  analysis  has  been 
obtained,  only  that  subset  of  the  total  data  base  which  is  compatible  with 
the  Massena  and  Tibbets  Point  data  base  is  analyzed  in  this  report.  Figure 
2-4  shows  the  Buoy  Auditing  Equipment  at  Iroquois  Lock. 


Figure  2-4 


Buoy  Auditing  Set  at  Iroquois  Lock 


2.3  "AD  HOC"  Loran-C  Data  Collection  Units. 


Because  of  the  Halted  availability  of  either  Harbor  Monitor  Sets  or 
Buoy  Auditing  Bets,  along  with  the  need  to  get  data  froa  several  sites  as 
soon  as  passible ,  extraordinary  aeans  were  taken.  USCG  had  additional 
Internav  404  Loran-C  receivers  available.  Texas  Instruaent  teleprinters 
were  also  available  for  recording  data  directly  froa  the  receivers  (without 
an  intervening  aicro-ccaputer) .  Figure  2-5  shows  the  prime  components  of 
the  systea  in  block  diagraa  fora. 


Figure  2-5  Ad  Hoc  Equipment  Set  Block  Diagram 


Via  front  panel  entry,  the  receiver  can  be  made  to  automatically  send  a 
aesaage  to  the  teleprinter  at  a  rate  of  some  integer  multiple  of  10  GRI. 
Originally,  the  setting  was  every  6020  GRI  which,  for  rate  9960,  resulted  in 
a  aesaage  about  every  10  ainutes.  At  that  rate,  cassette  tapes  had  to  be 
changed  about  every  10  days.  In  late  April  1982,  the  sampling  rate  was 
changed  to  one  sample  about  every  3  ainutes. 

The  output  aesaage  consists  of  an  instantaneous  sample  of  the  output  of 
each  servo  loop  (including  the  aaster  servo  loop)  along  with  the  receiver's 
current  eatiaate  of  the  signal-to-noise  ratio  for  each  station  being 
tracked.  The  entire  aessage  is  recorded  on  the  cassette.  As  a  general 
coament  on  the  receiver,  special  notice  should  be  made  of  the  aaster  servo 
loop  output  as  shown  in  the  column  labeled  "GRI”  in  the  figure*  of  Appendix 
F  (e.g.,  figures  F-l,  F-2) .  For  the  9960  chain,  the  reading  would  ideally 


be  99600.00.  A  deviation  from  this  reading  provides  the  operator  with  an 
indication  of  the  internal  oscillator  stability  (adjustment  recommended  if 
the  reading  differs  by  more  than  0.20  from  the  ideal).  More  importantly, 
the  receiver  software  applies  a  portion  of  this  offset  (computed  as  the 
ratio  of  the  TD  to  99600)  to  correct  the  measured  time  difference  for 
oscillator  error. 

Off-line  averaging  of  these  samples  can  be  accomplished  to  make  the 
data  more  nearly  compatible  with  the  HMS  data.  As  with  the  data  obtained 
from  the  buoy  auditing  sets,  only  the  subsets  of  the  total  data  base  which 
correspond  to  HMS  data  are  obtained  and  analyzed  in  this  report.  Figure  2-6 
shows  the  "Ad  Hoc"  Equipment  at  Beauharnois  Lock. 


Figure  2-6  Ad  Hoc  Equipment  Set  at  Beauharnois  Lock 


The  Ad  Hoc  equipment  approach  to  data  collection  is  suboptimal  for 
numerous  reasons  which  will  be  discussed  in  later  sections.  Primarily,  it 
requires  periodic  operator  assistance  and,  without  the  ability  to  remotely 
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access  the  data,  the  day-to-day  operation  of  the  system  cannot  be  properly 
monitored.  Nevertheless,  the  sets  were  deployed  at  CCG  sites  at 
Beauharnois,  and  Cote  Ste.  Catherine,  and  at  the  USCG  station  at  Tibbetts 
Point  in  late  December  1981.  Similar  CCG  equipment  was  installed  at 
Iroquois  Lock  on  17  March  1982.  Failure  of  the  teleprinter  at  Tibbetts 
Point  resulted  in  the  loss  of  data  at  that  site.  For  this  reason,  the 
equipment  was  replaced,  as  discussed  in  Section  2.1,  as  soon  as  a  harbor 
monitor  set  became  available.  The  equipment  at  Iroquois  Lock  was  replaced  on 
21  April  1982  as  discussed  in  paragraph  2.2  above.  A  series  of  problems  was 
encountered  at  Cote  Ste.  Catherine  so  this  site  was  shut  down  on  5  May  1982 
and  a  nearby  site  was  established  at  the  St.  Lawrence  Seaway  Authority 
facility  at  Brcmsard. 


The  locations  of  these  and  all  other  monitor  sites  are  indicated  in 
figure  2-7.  On  a  smaller  scale,  figure  2-8  shows  the  relation  of  the  Seaway 
to  the  stations  of  the  Loran-C  triad  which  provides  the  signals  examined  in 
this  report. 
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2.4  Improved  OSOG  Harbor  Monitor  Sets. 


As  mentioned  In  Section  1.3,  ths  Harbor  Monitor  Set  described  in 
paragraph  2.1  is  deemed  adequate  for  original  OSOG  Harbor  Monitor  project 
needs.  During  the  past  year,  however,  several  additional  applications,  for 
which  the  equipment  is  not  ideally  suited,  have  become  apparent.  Por  reasons 
such  as  this,  OSOG  has  taken  steps  to  develop  a  sore  general  purpose  set. 

The  set  features  an  off-the-shelf  micro-computer  (both  of  the  sets  described 
in  paragraphs  2.1  and  2.2  above  require  user  hardware  development),  will 
allow  for  high  density  data  storage  and  can  be  accessed/controlled  from  a 
remote  location  via  telephone.  Perhaps  more  importantly,  the  sets  have  been 
developed  to  incorporate  receiver  automatic  control  features  which  over  a 
year  of  harbor  monitor  experience  has  shown  desirable. 

Bigineering  development  of  these  units  was  completed  by  DSCG  Research 
and  Development  Center  personnel  in  early  May  1982  and  a  test  "field 
deployment"  was  started  in  mid-May.  These  units,  called  "Type  D  Harbor 
Monitors,”  will  be  available  for  other  deployments  in  the  August  1982 
timeframe  and  are  mentioned  later  in  this  report  for  consideration. 


2.5  Data  Collection  Set  Campatability . 


As  this  report  will  show,  it  is  possible  to  make  "low  density  analysis” 
comparisons  among  the  data  collected  by  the  myriad  equipment  deployed  thus 
far.  Hie  basic  approach,  however,  is  extremely  time-consuming  and  yields 
statistically  sub-optimal  results.  A  secondary  purpose  of  this  report  will 
be  to  address  the  suitability  of  each  type  of  set  and  make  recommendations 
for  future  data  collection  efforts. 


3. 


PERFORMANCE  PREDICTIONS. 


3.1  Need  For  A  Prediction  Model. 


In  conducting  a  Loran-C  stability  study,  one  is  initially  faced  with 
several  experimental  design  problems.  For  the  portion  of  the  St.  Lawrence 
Seaway  of  interest  -  nearly  two  hundred  miles  in  length  -  the  first 
important  question  involves  spacing  of  monitor  sites.  Although  no  formal 
deterministic  proof  can  be  offered,  it  is  reasonable  to  assume  that  if  a 
monitor  site  were  located  at  every  mile  along  the  river,  adequate  knowledge 
of  the  spatial  nature  of  the  variations  could  be  obtained.  As  a  practical 
matter,  and  besides  involving  prohibitively  expensive  site  installation 
costs,  it  would  not  be  possible  to  obtain  suitable  sites  in  all  cases  at 
such  a  spacing.  Additionally,  site  operation  and  maintenance  would  become 
an  unmanageable  effort.  Thus,  it  becomes  necessary  to  settle  for  a  much 
smaller  nianber  of  monitor  sites.  Even  with  only  six  sites,  the 
establishment  of  suitable  sites  -  particularly  in  a  rushed  timeframe  -  and 
operating  and  maintaining  the  site  equipment  -  become  non-trivial  tasks. 

(The  brief  history  of  the  sites  provided  in  section  2  givet  an  indication  ot 
the  scope  of  the  problems.  More  examples  will  be  shown  throughout  the 
report).  The  availability  of  data  from  six  sites  is  therefore  the  best  we 
can  hope  for  and  it  is  no  longer  reasonable  to  Mindly  assume  we  have  an 
adequate  number.  Clearly,  the  need  for  a  site  spacing  methodology,  such  as 
could  be  obtained  if  air  adequate  model  were  available,  exists. 

Another  experimental  design  problem  involves  the  time  between  samples. 
Again,  although  no  formal  proof  is  demanded,  all  involved  will  feel 
comfortable  with  samples  obtained  every  second  of  the  year  (including  leap 
seconds  where  appropriate)  from  each  site.  Even  with  only  six  sites, 
sampling  at  the  rate  mentioned  above  involves  the  accumulation  of  well  over 
1  billion  bytes  of  "data"  over  the  period  of  a  year.  In  these  days,  such  i 
data  base  can  be  constructed  at  a  small  cost.  This  having  been  said,  a 
disclaimer  is  in  order.  Whereas  all  involved  feel  comfortable  if  such  a 
data  base,  consisting  of  "known-to-be-valid-data"  is  presented  to  them, 
nobody  cares  to  assume  the  responsibility  to  assure,  or  pay  someone  to 
assure,  all  of  the  data  is  valid .  This  "data  base  integrity"  problem  (as 
opposed  to  data  base  storage  or  data  base  manipulation)  is  an  elementary 
concept  introduced  early  in  management  information  system  courses /seminal- 
fin  the  60's  the  buzzwords  were  different,  the  concept  was  known  as 
*garbage-in-garbage-out."  Only  the  size  of  the  problem  has  changed  -  hv 
orders  of  magnitude).  Numerous  examples  of  invalid  data  will  be  presented 
later  in  this  report  along  with  suggested  corrective  action  for  future  data 
collection  efforts.  It  will  suffice  at  this  point  to  simply  note  that  th- 
"once  per  second"  sampling  strategy,  given  available  resources,  exceeds  our 
analysis  (including  editing)  capabilities  by  several  orders  of  magnitude.  A 
less  intense  sampling  strategy  is  needed.  Again,  this  translates  to  a  ne°d 
for  a  model  of  the  variations. 


3.2  System  Requirements. 

Although  numerous  models  of  lor an  system  variations  have  been  developed 
over  the  years,  none  are  perfect.  Additionally,  most  are  not  directly 
applicable  to  the  types  of  questions  that  must  be  addressed  in  the  St. 
Lawrence  Seaway  study.  For  this  reason,  this  report  will  involve 
significant  discussions  of  modeling  techniques  which  the  authors  feel  are 
most  useful  for  our  purposes. 

It  is  note!  that  early  discussions  of  the  Steering  Committee 
concentrated  heavily  on  questions  of  sampling  strategy.  A  key  source  of 
disagreement  centered  on  the  fact  that  the  navigation  system  requirements 
were  not  properly  stated  (or  understood).  For  example,  as  late  as  February 
1982,  the  SLSDC  received  a  report  from  a  contractor  identifying  "Loran-C 
Requirements"  in  various  reaches  of  the  river.  For  many  of  the  reaches, 
(i.e.,  those  with  "half-channel  widths"  of  100  feet  or  less)  it  was 
indicated  that  Lor an -C  should  produce  a  cross-track  error,  99.9%  of  the  time 
less  than  about  7  meters.  At  about  the  same  time  (specifically,  at  the  17 
February  1982  meeting  of  the  Steering  Committee)  USCG  presented  the  plot 
shown  in  figure  3-1.  It  is  noted  here  that  figure  3-1  shows  no  half-channel 
width  of  less  than  200  feet  (61  meters) . 


Figure  3-1  Preliminary  Description  of  St.  Lawrence  Seaway 

Half-Channel  Widths 
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It  should  also  be  noted  that  USCG  was  neither  trying  to  "sell  Loran-C* 
nor  "kill  Loran-C."  It  was  simply  assumed  in  February  (and  will  be  shown 
herein)  that  Loran-C  cannot  give  7  aster  accuracy  99.9%  of  the  tiae  in  the 
St.  Lawrence  Seaway.  If  7  asters  (99.9%)  is  indeed  a  valid  Loran-C 
requirement,  then  this  report  will  serve  the  purpose  of  proving  Loran-C  will 
not  work.  If  however,  other  means  can  be  anployed  in  certain  reaches 
(street  lights  can  be  used  in  several  of  the  reaches  for  which  Loran-C  was 
stated  as  being  required),  then  the  job  becomes  more  difficult  and  we  are 
forced  to  continue. 

The  report  will  proceed  under  the  assumption  that  it  is  worthwhile  to 
consider  the  question:  "what  will  it  take  for  Loran-C  to  provide  accuracy  on 
the  order  of  25-30  aeters?”  Older  this  assumption,  vice  the  7  meter 
accuracy  criterion,  it  becoaies  possible  to  construct  a  reasonably  useful 
error  aodel  and  both  the  sampling  strategy  and  the  number  of  monitor  sites 
becomes  appropriate  for  answering  the  questions  being  asked. 


3.3  The  Nature  of  Loran-C  Positional  Errors. 

We  have  concluded  the  previous  section  without  explicitly  stating  the 
Loran-C  performance  requirement.  The  need  to  take  this  approach  can  be  seen 
by  examining  early  Steering  Committee  attempts  at  accuracy  requirement 
statements.  Specifically,  very  early  in  the  study,  considerable  time  was 
spent  trying  to  find  a  suitable  performance  metric.  Such  figures  of  merit 
as  "CEP"  or  "d-rms"  were  scrutinised  for  applicability.  As  part  of  the 
presentation  made  at  the  17  February  Steering  Committee  meeting,  USCG 
emphasized  the  concept  of  the  error  ellipse  -  the  pattern  of  positional 
errors  in  the  fixes  provided  by  Loran-C  and  most  other  positioning  systems. 
The  concept  is  well  illustrated  in  figure  3-2  which  ms  obtained  by 
processing  harbor  monitor  set  TO  data  using  the  math  of  Appendix  A.  The 
data  was  obtained  at  Groton,  Ct.  during  the  period  1  January  to  31  March 
1982. 


AVERY  POINT  RLDC  (X-V) 

CPI  9960  US  NORTHEAST 

JULIAN  DAY  Ml  TO  MS 

SI  JAN  88  TO  31  HAP  St 

AVC  RADIAL  ERROR  (MET EPS)  -  30.33 

RMS  RRDXAL  ERROR  (METERS)  -  37. 78 

MAX  RADIAL  ERROR  (METERS)  -  80.38 


Figure  3-2 


Fix  Scatter  Plot  of  Loran-C  Fixes  at  Avery 
Point  (Groton) 


Figures  of  aerit  such  «•  CEP  or  d-ras  ere  very  useful  in  considering 
large  scale,  general  accuracy  questions  such  as  those  involved  in  general 
lor an  chain  planning.  For  questions  relating  to  use  of  existing  loran 
coverage  in  the  St  Lawrence  Seaway,  however,  the  requireaents  are  too 
deaanding  for  general  aetrics.  Me  Bust  concentrate  on  specific  requireaents 
and,  for  the  Seaway,  this  aeans  cross-track  error.  As  can  be  seen  from  a 
quick  glance  at  figure  3-2,  the  error  pattern  at  Groton  is  predominantly 
oriented  at  a  bearing  of  about  000°  T.  If  the  course  of  a  particular 
channel  in  the  Groton  area  were  000°  T  (it  is!),  cross-track  errors  less 
than  about  10  asters  could  be  expected  alacet  all  the  tiae.  Alternatively, 
for  a  course  of  090°  T,  cross-track  errors  as  large  as  90  aeters  can  be 
reasonably  expected.  This  order  of  aagnitude  swing  in  error  as  a  function 
of  course  is  an  extreme  example  but  illustrates  the  significant  contrast 
with  the  "one-dimensional”  aetrics  such  as  CEP  or  drms  which  do  not  change 
as  a  function  of  course. 

In  considering  existing  loran  coverage  errors,  along  with  alternative 
configurations  (e.g.,  as  obtained  by  the  use  of  differential  corrections), 
we  cannot  afford  to  be  so  imprecise  as  to  consider  only  d-rms  or  CEP. 
Detailed  application  of  this  concept  will  be  explored  in  later  sections. 

For  now  we  simply  note  that  error  ellipses  and  associated  cross-track  error 
performance  should  be  estimated  for  each  reach  requiring  Loran-C.  This 
having  been  said,  we  can  return  to  a  discussion  of  system  requirements. 

In  1978,  the  SLSDC  received  a  contractor  report  (reference  3)  outlining 
the  characteristics  (including  course  and  minimum  channel  width)  of  the  103 
reaches  of  the  8eaway  from  Montreal  harbor  to  Tibbetts  Point.  In  March 
1982,  SLSDC  received  a  report  from  another  contractor  which  integrated  the 
reach  characteristic  table  with  a  preliminary  statement  by  the  Steering 
Committee  regarding  navigation  system  performance.  The  statement  was  that 
the  cross-track  error  of  the  navigation  system  should  be  5  aeters  less  than 
20%  of  the  channel  width  -  99.9%  of  the  time.  The  resulting  tabular 
presentation  of  requirements  is  included  in  Appendix  E.  At  the  May  1982 
Steering  Committee  meeting,  the  Seaway  Authority  presented  a  further 
refinement  of  this  table.  This  version  approaches  that  of  figure  3-1,  with 
more  detail,  and  is  also  provided  in  Appendix  E.  The  table  is  converted  to 
a  "half-channel-width"  plot  and  shown  in  figure  3-3. 
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Figure  3-3  Updated  Description  of  St.  Lawrence  Seaway 

Half-Channel  Widths 


The  two  "spikes"  shown  just  south  of  Beauharnois  represent  bridges,  in 
this  report  we  will  assume  loran-C  is  not  required  to  navigate  under  these 
structures . 

At  this  point  in  the  discussion  of  system  requirements,  i.e.,  having 
identified  cross-track  error  and  channel  width  as  key  ingredients,  we  move 
to  a  consideration  of  the  Loran-C  system  to  complete  the  discussion. 


3.4  Basic  Model  -  The  Double  Range  Difference. 

As  developed  in  Appendix  A,  we  can  generate  error  ellipses  for  any 
point  along  the  river  as  soon  as  we  can  obtain  a  measure  of  the  standard 
deviations  of  the  TD's  produced  by  the  two  major  Loran-C  baselines  in  the 
area  (9960-W  and  9960-X)  along  with  a  measure  of  the  correlation  coefficient 
between  the  two  baselines.  In  spite  of  the  fact  that  the  literature  is 
replete  with  referents  to  these  three  quantities,  it  is  unfortunately  true 
that  they  are,  strictly  speaking,  totally  fictitious.  In  obtaining  data 
from  the  monitor  sites,  we  are  obtaining  samples  of  random  processes.  Over 
the  period  of  a  year,  these  processes  cannot  be  considered  stationary  -  ever 
in  the  widest  sense.  Specifically,  this  means  we  are  collecting 
observations  of  random  variables  whose  means,  standard  deviations  and 
correlation  coefficient  change  from  time  to  time.  Thus  the  math  of  Appendix 
A,  or  any  number  of  similar  developments,  cannot  be  directly  applied. 

Moreover,  the  development  of  Appendix  A  and  related  treatments  assume 
the  variations  are  Gaussian  in  nature.  Again,  strictly  speaking,  this  is 
not  true.  Thus,  again,  direct  application  cannot  be  made. 

All  of  this  simply  means  that  we  must  be  clever  in  our  modeling  and 
analysis  techniques.  The  Gaussian  assumption  is  popular  because  it  very 
nearly  approximates  the  nature  of  variations  encountered  in  many 
applications.  As  will  be  developed,  we  need  to  break  the  total  data 
variation  model  down  into  its  constituent  parts  and  on  a  case-by-case  basis, 
determine  ways  to  approximate  the  Gaussian  model.  Although  numerous 
statistical  methods  exist  to  check  for  the  validity  of  the  Gaussian 
assumption,  we  will  shy  away  from  these  more  rigid  treatments  in  this  report 
and  simply  use  examples  extracted  from  the  data  base  to  illustrate 
applicability  of  the  model. 

With  proper  attention  to  details  of  the  model,  along  with  use  of  the 
concepts  provided  in  Appendix  B,  considerable  headway  can  be  made  in 
examining  the  suitability  of  Loran-C  to  satisfy  some  St.  Lawrence  Seaway 
navigation  requirements. 

The  major  constituents  of  the  Loran-C  time  difference  variations  are  as 
follows: 


a.  Seasonal  Component.  This  component,  in  most  cases,  is  the  largest 
contributor  to  the  net  error  budget.  It  has  a  period  of  1  year.  Examples 
of  thi3  type  of  variation  are  discussed  in  Appendix  B  and  shown  in  figure 
3-4.  If  the  time  series  of  figure  3-4. a  is  "low  pass  filtered"  by  use  of  a 
'-week  sliding  average,  a  "smoothed  TD"  data  plot,  as  in  figure 


3-5 


Bvery  Pt 


. . i . ! . i . . . 1 . * . - 


(a)  Actual  Data  From 
9960- W  Baseline 


. r  . t . t . I . r 

.—.4 . 4 . \ . 4 . .  — 


tt#  IM  IS#  III  •«#  |l|  M|  IM  MS 

Julian  Dsy  # 


flvery  Pt 


(b)  Sliding  21 -Day 
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(c)  "Residuals" 
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Figure  3-4  9960-W  Loran-C  Data  at  Avery  Point,  Ct. 


3-*.b  results.  The  peak-to-peak  value  (or  any  related  measure  such  as 
zero-to-peak  or  ras  value)  is  found  to  closely  follow  the  "modified" 
double -range  difference  (DRD)  model  which  is  described  in  detail  in  Appendix 
B.  Basically,  DRD  (measured  in  kilometers)  is  a  direct  measure  of  what 
hyperbolic  LOP  the  user  is  "on"  with  respect  to  the  hyperbolic  LOP  of  the 
System  Area  Monitor  (SAM).  If  the  observer  is  on  the  same  LOP  as  the  SAM, 
the  DRD  is  zero.  If  the  observer  is  on  an  LOP  which  is  closer  to  the  master 
section  than  the  SAM's  LOP  is,  the  DRD  is  positive.  As  developed  in 
Appendix  B,  we  have  found  it  necessary  to  not  count  saltwater  paths  in  the 
calculation  of  DRD.  The  regression  analysis  plotted  in  figure  3-5 
illustrates  this  "goodness  of  fit"  of  the  seasonal  variations  as  a  function 
of  DRD  for  the  9960-W  and  -X  baselines.  Given  the  slowly  varying  nature  of 


3 


.12 


SAMPLE  STATS 

>-«m  taa.  an 

X-Var  4TS4M.M 
X  STO  OEV  MSII 

T-VAM  4T888 1.667 

y  std  oev  was 


REGRESSION  RESULTS 
Slaf*  SMS 
•TO  term-  SMI 
Iw%»r.«^>-44. 7S8 
■TO  In-ar  18.617 
Na  *>  Pair*  4.M 
Can-.  Caa».  6.888 


»-  I L  242  ♦  1.M0X 


this  error  component,  the  twice  a  day  sampling  strategy  used  in  the  low 
density  analysis  is  entirely  adequate  for  measuring,  and  hypothesizing 
correction  schemes  for,  this  type  of  error. 

b.  Medium-Term  Component.  Here  we  are  talking  about  variation 
occurring  over  periods  ranging  from  several  days  to  several  weeks.  Examples 
can  be  seen  in  figure  3-4. c  in  which  the  seasonal  component  has  been  removed 
from  the  total  data  record  (accomplished  by  subtracting  the  data  in  figure 
3-4. b  from  the  data  in  figure  3-4. a).  It  can  be  seen  that  these  variations 
are  typically  smaller  than  the  seasonal  variations  and  that  the  "magnitude" 
of  the  variations  differs  from  season  to  season.  The  variations  are  most 
significant  in  the  winter  months  -extending  from  late  October  to  late  April 
in  northern  regions.  They  are  minimal  in  the  summer.  In  a  statistical 
sense,  these  variations  are  found  to  follow  the  modified  double  range 
difference  model  very  well.  Twice  a  day  sampling  is  adequate  to  measure  the 
variations  and  form  the  basis  for  a  correction  strategy. 

c.  Short-Term  Component.  These  variations  occur  over  periods  ranging 
from  several  hours  to  a  few  days.  The  fact  that  these  variations  can  occur 
as  rapidly  as  within  a  few  hours  is  the  prime  reason  for  the  concern  which 
has  been  expressed  by  the  Steering  Committee.  It  is  agreed  that  the  low 
density  sampling  strategy  is  not  adequate  for  a  deterministic  statement  of 
the  size  of  these  variations.  It  is  noted,  however,  that  significant 
variations  over  a  period  of  several  hours  are  rare  occurrences.  Moreover, 
the  variations  are  smaller  than  the  seasonal  or  medium  term  components. 
Finally,  since  most  of  the  examples  of  variations  assigned  to  this  category 
occur  over  the  period  of  a  day  or  90,  the  lew  density  analysis  provides  an 
adequate  statistical  measure  of  the  size  to  be  expected.  Concentration  on 
the  larger  components,  augmented  with  this  statistical  measure  of  the 
short-term  component,  is  argued  to  be  the  correct  strategy  for  proper  use  of 
the  resources  presently  available. 

d.  Near -Instantaneous  Variation  Component.  These  are  considered  to  be 
any  variations  which  occur  over  the  period  of  up  to  an  hour  or  so.  These 
variations,  for  practical  purposes,  are  considered  to  be  equipment  related 
and  not  associated  with  changes  in  signal  propagation  characteristics. 
Examples  can  be  seen  when  equipment  is  changed  at  transmitting  stations, 
when  corrections  are  made  to  chain  timing  by  the  control  station,  or  when 
there  are  receiving  equipment  problems.  These  variations  are  small  but 
cannot  be  considered  negligible  when  considering  the  capability  of  Loran-C 
for  precise  applications.  In  general,  no  attempt  need  be  made  to  measure 
these  types  of  variations  by  use  of  remote  monitors.  Observations  of  the 
effects  can  be  obtained  by  examining  chain  operating  records  and  by 
laboratory  measurements  or  simulations.  The  lab  measurements  should  be 
merged  with  the  field  measurements  to  obtain  the  total  error  description. 

FOr  the  purposes  of  this  study,  we  will  be  concerned  with  these  type 
variations  only  to  guard  against  their  effects  creeping  into  our 
.instrumentation  (i.e.,  make  sure  our  receivers  are  ins  tailed /opera  ted 
properly) . 


j.5  Application  of  Model  Using  A  Priori  Estimates 


Reference  4  featured  a  direct  application  of  the  double  range 
difference  model  described  in  Appendix  B.  An  ultimate  goal  of  the  St. 
Lawrence  Seaway  studies  will  be  to  duplicate  the  analyses  as  well  as  adding 
refinements  such  as  the  modified  range  calculation  presented  in  Appendix  B. 
Unfortunately,  for  reasons  already  alluded  to  and  to  be  illustrated  in 
section  4,  we  simply  do  not  yet  have  the  data  base  to  allow  direct 
application  of  the  model.  Thus  we  will  present  a  crude,  "first  cut" 
analysis  based  on  USCG  harbor  monitor  program  data  obtained  at  sites  outside 
the  St.  Lawrence  Seaway  region.  Since  these  performance  predictions  are 
|  made  without  considering  any  St.  Lawrence  Seaway  data,  we  call  them  "a 

priori"  estimates. 

The  presentation  will  be  primarily  heuristic  and  serves  the  purpose  of 
illustrating  the  role  of  key  parameters  that  must  be  estimated  before  the 
data  collection  effort  can  be  called  a  success.  An  auxiliary  purpose  is  to 
show  how  much  can  be  done  with  no  St.  Lawrence  data.  How  very  little  we  can 
presently  go  beyond  this  level  illustrates  hew  subpar  the  data  collection 
results  to  date  have  been. 

We  begin  by  noting  Appendix  A  allows  us  to  state  the  performance  at  any 
Seaway  reach  if  we  can  obtain  estimates  of c  w, 0  x  and  p  ^  -  the 
standard  deviations  of  the  two  prime  TD*s  in  the  area  and  the  correlation 
coefficient  between  the  TD's  -  at  any  point  of  interest.  We  have  previously 
noted  that  these  quantities  vary.  For  our  first  cut  at  estimating 
performance,  however,  we  will  ignore  this  technicality. 

Appendix  C  contains  the  details  of  the  estimation  procedure.  In  it  we 
first  look  at  the  seasonal  component  of  the  time  difference  variations.  We 
argue  that  this  component  is  so  predominant  in  the  St  Lawrence  Seaway  region 
that  we  can  closely  approximate  the  year-round  standard  deviations  by 
considering  only  this  component.  We  compute  estimated  peak-to-peak 
variations  for  each  TD  at  the  center  of  each  reach.  We  argue  the  seasonal 
variation  pattern  is  approx imatedly  sinusoidal  and  equate  the  rms  value  of  a 
sinusoid  with  the  indicated  peak-to-peak  value  to  the  appropriate  standard 
deviation.  Thus  we  obtain  estimates  of  ow  andOx  for  each  reach. 

To  estimate  the  correlation  coefficients,  we  argue  the  seasonal 
components  have  a  correlation  coefficient  approaching  unity.  Using 
arguments  suggested  by  Ligon  in  reference  5,  we  argue  the  "shorter-term" 
components  have  a  correlation  coefficient  of  about  0.5.  We  shew  how  a  total 
correlation  coefficient  is  formed  from  a  weighted  combination  of  these 
constituents.  A  total  correlation  coefficient  value  of  at  least  0.95  is 
moRt  likely  but  we  choose  a  more  conservative  estimate  (results  in  larger 
predicted  errors)  of  0.90  and  use  this  value  for  all  reaches.  Resulting 
"95%  probability  contour"  ellipses  are  generated  for  each  reach  -  and  listed 
m  Appendix  D.  The  projections  of  these  contours  on  corresponding 
cross-track/along-track  axes  are  computed  as  explained  in  Appendix  A.  Th“ 
resulting  cross-track  errors  are  listed  in  Table  3-1. 


•SLi 

WlfrMg 

01  f- 

000 

CTT  9*0  »1 

mo  # 

m 

017- 

000 

r 

ft/4  . 

. 

0 

at 

m-b.i/s-m.i 

a 

m 

M*T 

a  • 

Ill  * 

0 

w 

44-0 .0/70-14.0 

114 

214 

m 

a 

10 

0 

044 

44-Vl .4/79-14.0 

107 

Ml 

MM 

s-a.i/n-B.i 

itt 

10 

0 

Mil 

44-0 .0/79-17.2 

10 

10 

Ml 

0-19.1/ 73-M.4 

at 

10 

47 

04 

TO 

02 

_ 

ft/4  - 

0 

- 

ft/4  - 

...  - 

m.i 

M-1I.I/7VM.0 

a 

10 

0 

04 

44-40.2/79-19,4 

107 

at 

m 

a 

10 

0 

07 

m-m.Mn-m* 

a 

210 

tv 

ikivom  via  ii 

a 

10 

11 

(Ml 

44-41.4/74-0.4 

a 

07 

i/4  - 

B 

B1-I 

u-m.vn-H.t 

a 

m 

M.I 

0-11.0/74-0.7 

a 

14! 

0 

00 

M -44.1/S. M.I 

a 

07 

m 

a 

10 

0 

» 

IkM  «ltk  B 

114 

07 

ft/4  - 

. 

m 

B0 

M-41.I/S-M.I 

0 

BO 

07.1 

«.u.«m40.r 

a 

a 

0 

ta 

IS 

100 

M.I 

0-11.2/  34-12.2 

a 

ia 

47 

221 

44. ».»/*-»■* 

10 

8* 

w.f 

0-12.9/74-11.4 

in 

10 

0 

01 

•4-M.T/S-M.I 

IB 

223 

m 

44-11.4/74-10.4 

117 

172 

0 

B1 

44- ».  2/ 79-40. 4 

0 

m 

01.1 

*4  vn-nt 

0 

171 

a 

20 

44.M.4/S-4I.I 

H 

m 

19 

m.i 

44- 07. 4/ 74- M.  f 

73 

m 

41 

04 

M-S.I/S-tt* 

01 

217 

a 

m 

44-07.0/74-17.9 

0 

179 

0 

20 

44-0.1/79-41.2 

•1 

217 

V 

07.1 

im/H-ni 

0 

10 

0 

a? 

44.0.4/71-44.1 

0 

HI 

n 

M 

44-O4.4/W-H.0 

0 

191 

0 

ia 

44-0.1/79-0. ft 

0 

tu 

o 

n.< 

m-oi.i/m-».i 

0 

172 

« 

pi 

•S 

ID 

M 

w 

m 

113 

0 

m 

MM  «it*  M 

0 

... 

n 

01 

IHI.I/M-P.I 

a 

m 

n 

M 

04*0.0/71-0.4 

0 

Ml 

» 

09 

M-OD.  ?/>*-».  7 

a 

10 

m 

IW 

MlM  **«»  ct 

0 

Ml 

27 

ta 

44-0.1/74-0.4 

107 

10 

m 

20 

44-8.I/S-M.4 

0 

02 

a 

m 

0-0.2/74-41.4 

0 

•7 

TO 

20 

WM  alt*  » 

0 

nt 

» 

ta 

44-0.7/74-41.2 

10 

103 

71 

no 

44-S.I/S.M.I 

91 

142 

a 

07 

ikiihm  —  u  r 

a 

127 

n 

214 

laclaM  «1M  H 

01 

242 

ji 

M7 

M.8.4/M-4M 

107 

127 

71 

no 

IKlaM  alt*  n 

0 

Ml 

» 

-  ft/4 

74 

IS 

44-14.1/79-0.9 

0 

20 

n 

«f7 

44-r. 9/74-0. 2 

0 

147 

71 

04 

lacM  a«t»  x 

0 

242 

a 

07 

0-47.4/0-0.4 

91 

10 

0 

01 

IKIM4  «1t»  s 

0 

242 

* 

IB 

44-47.0/71-0.1 

a 

200 

77 

09 

M-U.t/S-8.1 

0 

no 

0 

09.1 

0-0.4/71-0.4 

no 

177 

a 

00 

44- 1ft.  4/79*0.  ft 

a 

20 

» 

m 

0-0.0/74-0.0 

1  JO 

10 

0 

04 

0-13.9/79^9.4 

in 

10 

a 

03 

0-0.4/71-0.0 

111 

104 

0 

M 

44-12.9/79-14.1 

IB* 

10 

a 

m 

0-0.1/74-07.1 

107 

10 

01 

10.4 

0-19. 9/ 71- 19. 4 

1»* 

27ft 

a 

BO 

0-0.4/71-0.4 

117 

219 

0 

ni 

44-0.I/79-B.1 

10 

27ft 

41 

07.1 

0-D.  4/74-10.4 

91 

10 

0 

04 

0-0.4/71-0.9 

01 

Ml 

0 

m 

0-0.0/71-12.2 

a 

10 

Table  3-1 


"Raw"  Loran-C  (9960-W/X)  Cross  Track  Error  (Preliminary)  95% 
Probability  Predictions  for  St.  Lawrence  Seaway  Reaches 
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A  comparison  of  the  predicted  cross-track  errors  and  corresponding 
half-channel  widths  shows  that  "un-augmented  Loran-C"  simply  cannot  yield 
satisfactory  year-round  performance:  the  error  ellipse  extends  well  outside 
the  channel  in  95*  of  the  reaches.  Thus  the  need  for  differential 
corrections  of  some  sort  is  indicated. 

Before  proceeding  to  a  consideration  of  differential  corrections,  we 
should  make  a  statement  about  the  probability  contour  used.  If  the  random 
variables  were  truly  Gaussian  (normal),  we  would  not  be  satisfied  with  the 
95%  contour  -  exceeding  the  channel  boundaries  5%  of  the  time  is  clearly  an 
unacceptable  performance  goal.  However,  we  have  already  noted  we  have  a 
distribution  which  follows  a  nearly  sinusoidal  pattern  -  i.e.,  one 
distinctly  different  from  the  Gaussian  distribution,  particularly  regarding 
extremes  of  variation.  We  can  illustrate  the  difference  by  comparing  a 
one-dimensional  sinusoid  to  a  Gaussian  random  variable.  A  gaussian  random 
variable  remains  within  two  standard  deviations  of  the  mean  about  95%  of  the 
time.  It  remains  within  three  standard  deviations  of  the  mean  about  99.75% 
of  the  time.  Thus,  if  a  time  difference  were  truly  Gaussian  distributed,  we 
would  expect  year-round  experimentation  to  show  that  if  the  deviation  from 
the  mean  exceeded  some  value,  say  X,  5%  of  the  time,  it  would  exceed  the 
value  1.5  X  0.25%  of  the  time. 

For  the  sinusoid-type  variation  we  expect  to  see  in  the  St.  Lawrence 
Seaway,  this  will  not  be  the  case.  For  a  pure  sinewave,  for  example,  simple 
trigonometry  shows  the  deviation  from  the  mean  which  is  not  exceeded  95%  of 
the  time  is  about  99.7%  of  the  zero-to-peak  value.  Thus,  if  X  again 
represents  the  deviation  from  the  mean  which  is  only  exceeded  5%  of  the 
time,  1.003  X  represents  the  value  never  exceeded.  We  simply  mention  all 
this  in  passing  -  along  with  an  observation  from  the  harbor  monitor  program 
that,  in  situations  wherein  the  seasonal  component  is  predominant,  the  "95%" 
error  ellipse  closely  approximates  the  maximum  error  (i.e.,  100% 
probability)  contour.  Thus  the  cross  track  errors  listed  in  Table  3-1  are 
good  estimates  of  the  maximum  expected  error.  The  sinusoidal  vice  Gaussian 
nature  of  the  variations  is  the  reason  why. 

We  should  now  move  the  discussion  to  predictions  of  differential 
Loran-C  performance.  As  outlined  in  Appendix  C,  we  hypothesize  a 
differential  monitor  station  at  the  Eisenhower  Locks  near  Massena,  N.Y. 

We  show  hew  it  is  necessary  to  re-compute  the  double  range  differences  and 
corresponding  standard  deviations  for  the  seasonal  component  of  the 
variations.  The  "differential  action,"  in  most  reaches,  removes  a 
considerable  portion  of  the  seasonal  component  of  the  variation  so  that  it 
is  no  longer  appropriate  to  ignore  the  shorter  term  components  in  estimating 
the  total  variation.  Additionally,  it  becomes  necessary  to  compute  a 
separate  estimate  of  the  total  correlation  coefficient  for  each  reach.  We 
hypothesize  a  representative  short  term  noise  standard  deviation  (20 
nanoseconds)  and,  using  the  formulation  of  Appendix  C,  devise  appropriate 
total  standard  deviation  and  correlation  coefficient  estimates.  Resulting 
error  ellipses  are  presented  in  Appendix  D  and  "95%"  probability  cross-track 
errors  are  listed  in  Table  3-2  in  the  column  labelled  "'Diff'  CTE  Pred  #  1." 
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Table  3-2 


Differential  Loran-C  Cross  Track  Error  (Preliminary)  95% 
Probability  Predictions  for  St.  Lawrence  Seaway 
(Differentia)  Corrections  From  Massena  ) 


(Note:  Thus  far  in  this  section,  as  well  as  in  Appendix  C,  we  have 
used  the  terms  "seasonal  component*  and  "shorter  term  component"  as  the 
basis  for  the  differential  Loran-C  predictions.  This  choice  of  words 
implies  a  similarity  with  the  categorization  of  the  "frequency"  components 
defined  at  the  end  of  Section  3.4.  For  the  purposes  of  the  discussions  in 
this  section,  there  is  no  reason  to  refute  this  similarity.  In  reality, 
there  is  a  subtle  but  important  distinction  to  be  made.  Hie  importance  of 
the  distinction  does  not  affect  the  discussions  until  Section  4.  Thus,  we 
will  avoid  pursuing  this  matter  until  Section  4.2  where  the  introduction  of 
actual  data  into  the  discussions  helps  illustrate  the  concepts  and  ease  the 
presentation. ) 


From  a  first  glance  at  Table  3-2  it  appears  the  Massena  differential 
monitor  solved  all  problems:  the  cross-track  errors  are  smaller  than  the 
corresponding  half-channel  widths  in  all  cases  -  with  considerable  error 
margin.  However,  with  the  large  seasonal  component  of  the  variations 
removed,  we  expect  the  "residuals"  to  more  closely  fit  the  Gaussian  mode1 
(although  still  not  perfectly).  Thus,  we  can  no  longer  use  the  "95%  is 
essentially  100%"  sinusoidal  case  argument  and  should  develop  a  more 
conservative  Gaussian  error  estimate. 

We  choose  the  figure  99.9%  -  a  somewhat  arbitrary  figure  -  for  a  number 
of  reasons.  First  we  should  note  we  still  have  a  non-Caussian  component  so 
that  we  really  expect  the  results  from  a  "99.9%  Gaussian  analysis"  to  be 
conservative.  Next  we  note  the  navigation  season  only  extends  from  about  1 
April  to  about  15  December  -  thus  excluding  most  of  the  worst  time  of  the 
"loran  year."  Whereas  use  of  the  95%  probability  figure  produces  overl.- 
optimistic  predictions,  use  of  some  figure  such  as  99.99+%  would  produce 
overly  pessimistic  predictions  (besides  being  a  number  which  exceeds  the 
availability  percentage  of  Loran-C  signals)  .  Neither  extreme  is  called  f 
at  this  point  in  the  study  and  the  99.9%  figure  is  offered  as  a  reasonable 
compromise.  Finally,  we  should  note  we  are  talking  about  a  single  AID  to 
navigation  -  not  an  input  to  a  robot  navigator.  The  99.9%  figure  leaves 
room  (about  9  hours  per  year,  6  hours  per  shipping  season)  for  the  navigat  r. 

Using  the  formulas  of  Appendix  A,  we  find  the  "size"  of  the  ellipse  >.s 
increased  by  a  factor  of  1.52  in  going  from  a  probability  of  95%  to  a 
probability  of  99.9%.  Cross-track  errors  for  this  predicted  probability  are 
listed  in  Table  3-3  in  the  column  called  " 'Conservative '  Differential  CTE 
Pred  #  1."  Hie  values  in  the  column  labelled  "error  margin"  are  computed  by 
subtracting  the  conservative  cross-track  error  predictions  from  the 
half-channel  widths. 

In  examining  the  table,  we  see  the  "conservative"  error  prediction 
never  exceeds  the  half-channel-width.  Thus,  the  center  line  of  a  vessel 
following  the  loran  with  no  guidance  error  never  goes  outside  the  channe1  - 
i-  the  99.9%  probability  level.  We  must  note,  however,  that  the  "half-wid*-h" 
of  the  largest  vessels  can  be  17  meters  in  the  Seaway.  Thus,  17  meters  i*=  a 
nunber  of  concern.  We  place  three  asterisks  in  the  table  for  any  reach  wot 
in  error  irargin  less  than  this  value.  Three  asterisks  thus  indicate 
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Table  3-3  Differential  Loran-C  Cross-Track  Error  Predictions  - 

95%  and  99.9%  Probability 


4  situations  wherein  if  the  conservative  error  estimate  is  exceeded,  part  * 

vessel  with  a  half-width  of  17  meters,  following  the  loran  with  no  guidance 
error,  is  outside  the  channel.  We  cannot  .  ~>f  course,  expect  zero  ship 
guidance  error.  We  assign  two  asterisks  a  reach  in  which  the  error 
margin  is  less  than  22  meters.  Thus  two  ar -■’risks  indicate  situations  in 
which  if  the  conservative  error  margin  estimate  is  exceeded,  part  of  a 
4  vessel  with  a  "half-width"  of  17  meters,  fol  ~>wing  the  loran  with  5  meters 

or  more  guidance  error,  is  outside  the  channel.  One  asterisk  indicates  a 
situation  in  which  the  conservative  error  es  mate  leaves  less  than  10  meters 


I 
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for  guidance  error  for  a  17  meter  "half-width"  vessel.  We  do  not  offer 
these  criteria  as  the  ultimate  design  goal  of  the  study  -  just  as  some 
figures  of  merit  for  the  initial  examination  of  predicted  differential 
Loran-C  performance. 

In  Table  3-3  we  see  two  areas  where  asterisks  are  clustered  -  both 
being  considerable  distances  frcm  Massena.  This  suggests  we  consider 
adding  another  monitor  station  -  in  the  vicinity  of  reach  74  at  Wellesley 
island.  If  we  re-compute  the  predicted  errors  based  on  differential 
corrections  from  this  station  (for  appropriate  reaches  only),  we  obtain  th 
results  indicated  in  Table  3-4. 


"Conservati **" 


Retch  # 

Course 

Htlf- 

Wldth 

CTl  Prtd  #1 

*01  ff* 

Crt  Pred  #2 

Differential 
CTt  Pred  #2 

Error 

Hergln 

56 

238 

152 

208 

28 

43 

109 

57 

223 

152 

229 

27 

41 

111 

58 

231 

152 

223 

27 

41 

111 

59 

223 

83 

235 

27 

41 

42 

60 

233 

76 

222 

28 

43 

33 

61 

236 

61 

217 

27 

41 

20  ** 

62 

222 

61 

237 

29 

44 

17  ** 

63 

217 

99 

241 

28 

43 

56 

64 

194 

93 

233 

25 

38 

55 

65 

218 

93 

237 

29 

44 

49 

66 

216 

93 

49 

67 

209 

93 

243 

28 

43 

50 

68 

207 

93 

50 

69 

225 

93 

232 

29 

44 

49 

70 

228 

93 

49 

71 

218 

91 

242 

29 

44 

47 

n 

214 

61 

17  ** 

73 

219 

68 

24  * 

74 

235 

68 

220 

28 

43 

25  • 

75 

235 

68 

25  • 

76 

231 

68 

25  • 

77 

239 

83 

216 

28 

43 

40 

78 

228 

91 

238 

29 

44 

47 

79 

246 

111 

199 

29 

44 

67 

80 

263 

122* 

148 

26 

39 

83* 

81 

193.5 

122* 

275 

28 

43 

79* 

82 

213 

165 

275 

32 

49 

116 

83 

234 

381 

242 

32 

49 

332 

Table  3-4 


Revision  to  Table  3-3  Based  on  an  Additional  Differential 
Monitor  at  Wellesley  Island 


Comparing  these  results  with  those  shown  in  Table  3-3  suggests  it  is 
appropriate  to  switch  from  Massena  to  Wellesley  Island  for  corrections  after 
reach  55.  Table  3-5  is  constructed  using  the  two  station  differential 
control  strategy  thus  obtained.  We  see  we  no  longer  have  "three  asterisk" 
situations  so  the  attention  must  focus  on  how  much  room  must  be  allowed  for 
guidance  error.  Note  that  except  for  reaches  61,  62,  and  72,  at  least  7 
meters  is  allowed  -  a  figure  which  could  be  argued  to  be  very  close  to 
acceptable . 
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Table  3-5 


Differential  Loran-C  Cross-Track  Error  Predictions 
Monitors  at  Massena  and  Wellesley  Island 


Before  we  decide  to  try  to  locate  another  differential  Monitor  station 
in  the  problematic  reaches  we  should  consider  the  source  of  the  problem.  As 
Table  C-3  in  Appendix  C  shows,  with  differential  corrections  being  applied 
frem  Wellesley  Island,  there  is  a  slight  seasonal  component  which 
contributes  to  the  total  error  budget  in  reaches  61,  62,  and  72  through  76. 
Without  trying  to  find  a  particular  site  which  removes  this,  we  should 
immediately  move  to  consideration  of  an  "optimal"  monitor  site  -  one  which 
results  in  the  removal  of  all  seasonal  variations.  Error  ellipses  for 
reaches  61  and  62  under  such  conditions  are  shown  on  page  D-32  which  show 
the  predicted  95%  cross-track  error  is  28  meters.  This  translates  to  a 
99.9%  probability  prediction  of  43  meters  -  leaving  an  error  margin  of  only 
18  meters. 

Thus,  under  our  assumptions  regarding  the  short-term  standard 
deviations,  we  "bottom  out"  at  the  99.9%  level  in  these  tight  reaches.  Our 
assumptions,  of  course,  nay  be  wrong  so  we  cannot  draw  final  conclusions  at 
this  point.  It  is  important,  however,  to  emphasise  the  source  of  the 
problem  which  becoams  clear  from  examination  of  the  error  ellipses  which  are 
reproduced  in  figure  3-6  below.  A  general  note  about  the  error  ellipse 
plots  which  appear  throughout  the  remainder  of  the  report  is  in  order.  The 
vertical  axes  of  the  plots  have  a  true  north-south  orientation  and  the 
horizontal  axes  have  a  true  east-west  orientation.  The  straight  line  which 
extends  from  the  origin,  shown  below  at  a  heading  of  236°T  for  Reach  #61, 
indicates  the  channel  course.  Characteristics  of  the  ellipse  and  the 
statistics  from  which  it  was  generated  are  provided  below  each  plot. 


bench  ti  optimrl  north 


Prfltotbilitr  of  93  '*  6'gr>*  M  -  .099  u**c 
•ho  -  .300  Stgro  X  •  .Ml  u**c 
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Hoe  RTC  •  11  Motoro  H%x  CTE  ■  29  Motor* 


Figure  3-6  Differential  Loran-C  Error  Ellipses  at  Reaches 

61  and  62 
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The  orientation  of  the  ellipse  with  respect  to  the  course,  in  these 
critically  tight  reaches,  is  essentially  "worst  case."  If  the  orientation 
were  rotated  90°,  we  would  have  cross-track  errors  of  11  or  12  meters  and 
there  would  be  no  problem. 

The  ultimate  conclusion  may  be  that  Loran-C  will  not  work. 
Alternatively,  we  may  decide  to  deploy  some  other  system  for  these  special 
reaches.  He  should  note,  however,  the  criticality  of  system  geometry  as 
well  as  the  importance  of  obtaining  good  data  in  the  vicinity  of  Wellesley 
Island. 


3.6  Summary  of  Performance  Predictions. 

In  Section  3  we  have  introduced  a  model  for  Loran-C  time  difference 
variations  for  use  in  the  St.  Lawrence  Seaway  study.  He  argued  that  the 
variations  can  be  broken  down  into  four  "frequency"  components.  Two  of 
these,  it  is  claimed,  can  be  adequately  measured  -  deterministically  -  by 
the  "low-density"  sampling  strategy  used  in  the  harbor  monitor  program.  For 
the  third  component,  it  is  argued,  the  low-density  sampling  strategy  allows 
an  adequate  statistical  measure.  The  remaining  component  is  called 
"near-instantaneous"  and  it  is  claimed,  since  this  is  not  "propagation 
phenomenon  related,"  it  is  not  necessary  to  deploy  spatially  separate 
monitors  in  the  service  area  to  monitor  this.  If  the  need  is  indicated, 
estimates  of  the  "rapid"  effects  can  be  obtained  from  chain  operational 
monitors.  Hhereas  it  has  never  been  demonstrated  that  these  "rapid"  effects 
are  perfectly  monitored  at  some  location  other  than  the  point  of  interest, 
that  also  is  something  which  should  be  pursued  only  if  the  need  is  indicated 
-  and  in  an  extremely  cautious  and  careful  manner.  For  now,  we  must  try  to 
avoid  "swallowing  any  camels  while  straining  for  gnats."  As  the  next 
section  will  show,  there  appear  to  be  many  "camels"  in  the  St.  Lawrence 
Seaway . 

To  consider  the  spatial  nature  of  the  time  difference  variations  we 
introduce  the  so-called  "double-range-difference"  model  which  has  been  use^ 
in  similar  studies  of  this  nature.  The  model  is  described  in  Appendix  B. 

He  use  data  obtained  from  sites  outside  the  St  Lawrence  area  to  apply  the 
model.  In  other  applications,  a  more  empirical  use  of  the  model  is  made  - 
and  nothing  presented  here  should  be  construed  to  imply  we  do  not  believe 
this  is  the  desired  approach.  Unfortunately,  we  do  not  yet  have  the  data 
base  necessary  to  support  this  type  of  approach  and  had  to  turn  to  more 
"heuristic"  methods.  We  did  this  to  emphasize  there  is  a  theory  for  t  e 
types  of  variations  we  expect  to  see  and  that  its  presentation  is  too 
important  to  wait  for  an  adequate  data  base.  In  applying  the  model,  we  did 
not  have  a  convenient  way  to  estimate  all  four  frequency  components  so  we 
considered  only  two  types  of  components.  In  Section  4  we  will  further 
explore  the  subtle  differences  in  the  tv«  different  types  of  classification 
of  variations.  It  is  noted  that  in  ignoring  the  distinction  for  the  moment 
we  do  not  compromise  the  validity  of  the  predictions  of  this  section.  The 
presentation,  crude  as  it  may  be,  serves  to  focus  our  attention  on  the 
critical  parameters  and  areas  we  must  concentrate  on  for  the  rest  of  the 
study . 
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The  predictions,  which  we  call  "a  priori*  since  they  are  aade  without 
the  benefit  of  observations  in  the  region  of  interest,  show  that  the 
performance  to  be  obtained,  year-round,  with  "raw*  Loran-C  is  clearly 
inadequate  for  our  requirements.  Although  we  do  not  have  a  simple  method  of 
separating  the  shipping  season  fraai  the  entire  year,  the  year-round  errors 
are  large  enough  for  us  to  easily  argue  "raw"  Loran-C  will  be  inadequate 
throughout  the  shipping  season.  Similarly,  the  errors  are  large  enough  to 
argue  a  "simple”  solution  (e.g.,  the  application  of  monthly  or  weekly 
corrections)  does  not  appear  promising.  At  this  point  it  is  not  productive 
to  speculate  on  the  viability  of  such  schemes  as  daily  corrections  so  we 
move  directly  to  consideration  of  "full  Differential  Loran-C,"  i.e.,  we 
assume  corrections  frcm  the  monitor  are  available  to  the  user  on, 
essentially,  a  real-time  basis. 

Hypothesising  a  differential  control  site  at  Eisenhower  Locks,  we  find 
adequate  performance  -  at  the  99.9%  probability  level,  is  predicted  for  most 
sites.  A  cluster  of  problematic  reaches  is  seen  in  the  vicinity  of 
Wellesley  Island.  We  hypothesize  a  differential  monitor  at  Wellesley  Island 
and  find  the  errors  are  still  "too  large  for  comfort"  in  three  of  the 
reaches.  Moreover,  we  find  the  errors  are  predicted  to  be  too  large  no 
matter  where  we  locate  a  differential  correction  station.  This  finding  is 
based  on  an  arbitrary  assunption  regarding  the  short  term  statistics  of  the 
variations  so  that  we  cannot  rule  out  Loran-C  at  this  time.  We  can, 
however,  be  aware  of  the  "nearly  worst  case"  effects  of  chain  geometry  which 
is  identified  as  the  prime  culprit. 

If  data  shows  our  assumptions  about  the  time  difference  statistics  were 
not  overly  pessimistic,  the  situation  is  not  all  that  bad.  Fortuitously, 
the  areas  of  concern  are  located  in  reaches  in  which  the  proximity  of  small 
islands  leads  to  consideration  of  a  "hybrid  system"  solution.  RACONS,  for 
example  could  be  located  on  the  islands  to  "mark”  the  channel  edges.  The 
entire  exercise  shows,  however,  that  consideration  should  be  given  to 
locating  a  monitor  in  the  Wellesley  Island  vicinity. 

Now  that  we  have  made  all  possible  predictions,  we  should  turn  to  a 
consideration  of  available  data.  In  our  considerations,  we  should  first 
attempt  to  verify  that  the  predicted  "raw"  Loran-C  error  ellipses  are 
approximately  representative  of  what  actually  happens.  If  this  can  be 
accomplished,  we  can  verify  our  conclusion  that  with  existing  chain 
geometry,  "true"  Differential  Loran-C  must  be  employed.  Next,  we  should  try 
to  verify  the  differential  Loran-C  predictions  -  particularly  in  the 
critical  reaches  near  Wellesley  Island.  Although  we  do  not  expect  available 
data  to  allow  us  to  say  the  final  word  at  this  time,  examining  how  much  of 
the  predictions  of  this  Section  can  be  tested  provides  an  assessment  of  the 
effectiveness  of  the  experiment  to  date. 
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4. 


DATA  COLLECTION  RESULTS 


4.1  Data  Availability 


Data  collection  for  the  St  Lawrence  Seaway  Loran-C  stability  study 
began  at  Eisenhower  Locks  (Massena,  N.Y.)  on  21  October  1981.  As  described 
in  Section  2.1,  data  at  Nassena  is  collected  every  12  hours.  During  each 
sampling  period  we  measure  and  record  the  hourly  average  value  and  standard 
deviation  of  the  9960-W  and  9960-X  time  differences.  We  also  record  the 
minimum  and  maximum  value  observed  during  the  sampling  period.  All  of  this 
data  is  stored  along  with  notations  as  to  the  date  and  time  of  the 
observations.  In  writing  this  report  we  have  used  data  up  to  11  June  1982 
for  Massena  and  thus  have  a  total  of  over  450  sample  periods.  Although 
equipment  problems  have  caused  the  loss  of  some  data,  we  find  we  have 
obtained  98%  of  the  desired  data  at  Massena. 

A  98%  complete  record,  over  a  period  of  almost  8  months,  should  allow 
an  adequate  "low  density"  data  analysis  -  particularly  if  the  period 
includes  all  of  the  winter  as  the  Massena  record  does.  Unfortunately, 
Massena  is  the  only  site  in  the  St  Lawrence  for  which  we  have  so  much  data. 
Thus,  we  can  make  fairly  well  substantiated  statements  about  "raw  Loran-C," 
but  have  problems  in  trying  to  verify  differential  Loran-C  performance 
predictions.  The  problems  are  beat  seen  by  examining  Table  4-1  which  shows 
the  "%-usable"  (low  density)  data  available  to  the  authors  at  the  time  of 
this  writing  for  the  period  from  21  October  1981  to  11  June  1982. 


Site 

Total 

%-usable 

Massena 

98% 

Tibbetts  Point 

52% 

Iroquois 

20% 

Pr escott 

7% 

Beauharnois 

4% 

ote  Ste  Catherine 

0% 

Broesard 

0% 

Table  4-1 


Percent  Usable  Data  from  St.  Lawrence  Seaway  Loran-C  Sites 


A  reason  for  low  data  availability  at  some  sites  (e.g.,  Brossard)  is 
the  fact  that  the  installations  took  place  late  in  the  October  '81  -  June 
'82  period.  The  fact,  for  example,  that  we  expect  a  much  higher  performance 
record  from  the  Brossard  site  over  the  next  few  months  does  nothing  to 
alleviate  the  plight  of  the  authors  who  are  trying  to  report  -  now  -  on  the 
results  of  the  8-month  effort. 

Another  reason  for  the  low  data  availability  reflected  in  Table  4-1  for 
sites  other  than  Tibbetts  Point  and  Massena  is  the  substantial  "time  delay" 
involved  in  obtaining  the  data.  As  an  example,  we  should  note  that  as  of 
"close  of  business"  on  21  June  1982,  data  from  Massena  and  Tibbetts  Point  up 
to  and  including  the  first  hourly  sample  taken  on  21  June  1982  was  obtained 
and  examined  by  project  personnel  at  the  USCG  Research  and  Development 
Center.  This  data  was  available  to  the  authors  that  same  day  (it  was  not 
used  due  to  a  self-imposed  deadline  regarding  the  preparation  of  plots) . 

This  situation  contrasts  significantly  with  that  associated  with  the  data 
from  other  sites.  As  of  21  June  1982,  no  data  beyond  18  May  1982  has  been 
received  and  reviewed.  The  great  question  regarding  the  data  from  these 
other  sites  is:  "suppose  'something  broke'  on  19  May  1982?"  This  is  not 
merely  a  hypothetical  question:  situations  such  as  this  have  already 
occurred  several  times  during  the  course  of  this  project.  If  this  has 
happened,  we  have  lost  more  than  a  month  of  data  unless  something  has  "fixed 
itself."  Those  with  a  mild  technical  background  (e.g.,  those  who  own  TV 
sets)  can  appreciate  the  unlikeliness  of  this  happening. 

A  final  major  reason  for  low  data  availability  is  a  feature  found  only 
at  the  "Ad  Hoc"  sites.  Without  a  computer-based  controller  of  some  sort, 
the  data  is  recorded  without  corresponding  date/time  information.  The  data 
is  obtained  at  some  integer  multiple  of  0.0996  seconds  so  that,  ideally,  one 
could  reconstruct  the  times.  If,  however,  power  :s  temporarily  interrupted, 
or  if  equipment  malfunctions  occur,  we  cannot  reconstruct  the  date  and  time 
of  the  measurements  unless  an  operator  re-initializes  the  system  and  is 
meticulous  about  manually  recording  the  date  and  time  of  re-initialization. 
Unfortunately,  we  do  not  have  technically  skilled  operators  at  the  Ad  Hoc 
equipment  sites  and  through  this  "improper  initialization"  mechanism,  have 
lost  a  considerable  amount  of  data.  In  many  cases,  the  Loran-C  data  is 
available  -  "time  of  measurement  unknown"  -  and  this  could  be  used  to  make 
some  statements  about  "raw  Loran-C"  performance.  Given  the  availability  of 
data  from  Massena,  however,  we  will  be  able  to  say  all  we  need  to  about  "raw 
Loran-C"  so  this  is  redundant  information.  The  true  value  of  data  from 
sites  other  than  Massena  is  that  it  could  allow  measurements  of  differential 
Loran-C  performance.  Accurate  knowledge  of  the  date  and  time  of  the 
observation,  however,  is  critical  for  this  analysis  and,  in  many  cases,  this 
vital  information  has  been  lost. 

Table  4-2  provides  a  breakdown  of  the  "unusable  data"  causes  for  each 
site.  From  this  we  can  see  there  are  so  many  problems  with  the  Ad  Hoc  data 
sites  that  we  must  seriously  consider  upgrading  the  sites  with  more  reliable 
equipment.  The  major  problem  with  the  buoy  auditing  equipment  sites  are  all 
related  to  untimely  collection  and  review  of  the  data.  Serious 
cons ider at  uxi  should  be  given  to  upgrading  the  equipment  at  these  sites  to 
allow  remote  (via  phone)  access. 


sit# 


Cowrit# 


Prescott /Ft  Wellington  a.  Buoy  Auditing  Bgulpasnt  Sst  installed  at  Port  Wellington  2/10/82. 

(2/10/82  -  8/8/82)  b.  2/10/82  -  3/ll/82i  Data  not  valid  dua  to  goundlng  problems. 

o.  1/11/82  -  l/29/82i  Data  la  valid  but  faaturaa  savaral  1-1/2 
to  2  day  gaps  whenever  tapaa  ran  out. 

d.  1/10/82  -  4/20/82t  No  valid  data  dua  to  equipment  and  power 
failures. 

a.  4/20/82i  aquipsMtnt  Moved  team  Port  Nall lngton  to  CCG  Baaa  at 
Prescott.  TO  offsets  recorded. 

f.  4/20/82  -  4/24/82:  valid  Data. 

g.  4/25/82  -  6/8/82)  data  not  valid  dua  to  receiver  probleaa. 

h.  6/8/82<  receiver  replaced,  only  9960-W  and  9960-X  data  being 
recorded. 

&uanary>  Insufficient  data  for  a  Meaningful  low  density  analyaia. 

Of  tiia  data  collected  during  tna  period  2/10/82  -  6/8/82, 

14%  is  usable. 


a.  Installed  1/17/82  aa  an  Ad  Hoc  type  site. 

b.  3/17/82  -  4/20/82)  good  data  with  only  one  7-hour  period  of  loet 
data  -  due  to  a  'cycle  slip.* 

c.  4/21/82)  Buoy  Auditing  Sat,  with  a  signal  splitter,  Installed. 

d.  4/21/82  -  S/1/82)  Good  data  with  the  eignal  splitter. 

a.  S/1/82  -  S/14/82 f  Data  lost  due  to  lack  of  re-lnitlal istlon 
after  a  power  loss. 

f.  About  6/1S/82  >  Splitter  r  a  Moved  upon  USCC  recoaseendation. 

Offsets  recorded. 

SUMwaryi  Iroquois  Lock  now  features  a  single  receiver  using  one 
antenna.  9940-W  and  9960-X  data  being  recorded. 

Of  the  data  collected  during  the  period  3/17/82  -  4/14/82, 
78%  la  usable. 


a.  Installed  12/21/81  aa  an  Ad  Boo  type  eite. 

b.  12/22/81  -  4/2/82)  Data  not  valid  due  to  poor  grounding. 
9960-X  TD'a  and  SNA  readings  unstable. 

C.  4/2/82)  laproved  grounding. 

d.  4/2/82  -  5/5/82 >  About  75%  Of  the  data  la  Invalid  dua  to  a 
'cycle  slip"  on  the  9960-X  baseline. 

Susswry)  Insufficient  data  for  a  Meaningful  low  density  analyaia. 
Of  the  data  collected  during  the  period  12/21/82  - 
S/S/82,  6%  la  usable. 


a.  Installed  12/22/81  aa  an  Ad  Hoc  type  site. 

b.  12/22/81  -  4/1/82)  Date  Is  not  valid  due  to  poor  grounding. 
TO'e  and  SNA's  unstable. 

c.  4/1/82)  Grounding  ieproved.  Subsequently  discovered  that 
the  receiver  "went  bed*  ecasetia*  between  4  and  11  Pwb.  1982. 

Suawaryi  None  of  the  data  le  valid. 


a.  S/4/82)  Installed  as  an  Ad  Hoc  type  alts. 

b.  S/S/82)  No  valid  dace  due  to  grounding  probleaa.  9960-X  TO'e 
end  SNA's  unstable. 

SUNMry)  No  further  date  available  by  'press  ties.* 


Table  4-2  Comments  or  Site  Performance 


Specific  examples  of  the  data  obtained  duriny  the  problem  periods 
>ot lined  above  are  provided  in  Appendix  !•'. 


at  jsaard 

4/82  -  5/5/82) 


■  *  ate.  Catherine 
..  22/81  -  5/5/82) 


Beauhemola  book 

<12/21/81  -  5/5/82) 


Iroquois  book 
<3/17/82  -  5/14/82) 


4-3 


4.2  Evaluation  of  Massena  Data. 


We  begin  the  evaluation  by  considering  the  data  collected  at  the 
Massena,  N.Y.  site.  Plots  of  the  time  difference  readings  are  provided  In 
figures  4-1 .a.  and  b.  The  data  Is  plotted  from  Julian  Day  1  to  Julian  Day 
365  to  be  consistent  with  the  plots  Introduced  In  Section  3  and  the 
Appendices.  Thus,  the  first  (read  from  left  to  right)  155  days  of  data 
shown  below  comprise  1982  data  whereas  the  data  plotted  to  the  right  is  from 
1981. 


i  ie  6B  ae  ite  is*  its  21a  t*s  t?»  sea  im  ]«e 

Julian  Day  # 

n««n  -  14721.363  Sign*-  .133 


Julian  Day  * 

Ham  -  27*31.320  Stgiia  -  240 


Figure  4-1 


9960-W/X  Loran-C  TD  Plots  -  Massena 


There  are  Many  things  we  can  say  about  the  data  records  but  we  should  start 
by  comparing  the  statistical  results  to  the  year-round  predictions  for  "raw 
Loran-C"  of  Section  3.  For  the  9960-W  baseline,  we  predicted  a  standard 
deviation  of  0.144  usee.  We  see  the  results  -  for  almost  8  months  of  data 
show  a  "sigma"  of  0.133  usee.  For  9960-X,  we  predicted  a  standard  deviation 
of  0.255  usee  and  observe  a  "sigma"  of  0.248  usee.  We  will  have  to  wait 
several  months  for  the  final  judgment  on  the  predictions  but  note  we  expect 
both  data  records  to  remain  fairly  stable  throughout  the  Summer  and  early 
Fall.  Thus  we  can  "dub  in"  a  sequence  of  data  points  which  follow  the  basic 
trend  established  by  mid-June  and  compute  the  statistics  of  the  results. 

They  turn  out  to  yield  "sigmas"  of  0.137  usee  and  0.256  usee  for  9960-W  and 
9960-X,  respectively. 

From  this  we  can  conclude  the  time  differences  at  Massena  agree  with 
the  pattern  suggested  by  other  N.E.U.S.  data  collection  sites  within  5%.  In 
short,  there  are  no  surprises  at  Massena  regarding  the  year-round  standard 
deviations.  The  data  points  of  figures  4-1. a.  and  b.  are  converted  to 
position  errors  to  produce  the  "scatter  plot*  shown  in  figure  4-2. 


■■  -iia  ■  in  miifi 
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Figure  4-2  9960-W/X  Loran-C  Error  Fix  Scatter  Plot  -  Massena 


In  the  process  of  generating  the  ellipse  that  is  included  in  figure 
i  2,  nave  computed  the  correlation  coefficient  of  the  two  time  difference 
■cords.  At  0.942,  the  observed  value,  as  expected,  is  higher  than  the  0.90 
"  ie  we  used  in  Section  3  for  the  "raw  Loran-C"  predictions. 


Prat  the  results  presented  thus  far,  we  can  draw  several  important 
conclusions : 

a.  The  correlation  coefficient  prediction  was  low.  As  noted  in 
Section  3,  this  was  the  result  of  a  conscious  attempt  to  generate 
conservative  predictions.  The  standard  deviations  we  have  observed  are  not 
different  frcn  the  predicted  values  in  any  statistically  significant  sense. 
This  is  a  critical  result.  It  suggests  that  with  the  simple  "land/seawater 
modification"  to  the  uniform  propagation  assumption  described  in  Appendix  B, 
the  double-range  difference  model  applies  very  well  to  observations  in  the 
vicinity  of  Massena,  N.Y.  Thus,  we  have  no  evidence  of  propagation 
anamolies  thus  far. 

b.  The  error  pattern  of  the  actual  fixes  does  not  perfectly 
follow  the  "pure  ellipse”  pattern.  This  can  be  considered  primarily  a 
result  of  the  non-gaussian  nature  of  the  processes  involved.  To  a  "very 
useful  degree,"  however,  the  pattern  is  as  predicted  and  the  appropriateness 
of  the  "elliptical  analysis"  is  verified. 

c.  In  figure  4-3. a.  and  b.f  we  compare  observed  results  with 
predicted  results.  Table  4-3  summarizes  the  comparison. 
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Figure  4-3  Comparison  of  Predicted  and  Observed  Error  Ellipses 

Massena 


Max  era 

Max  ATE 

Predicted  95%  Ellipse 

157  m 

136  m 

95%  Ellipse  Determipeg 

by  Obsecved  Statistics 

144  m 

130  m 

Actual  Observations 

135  m 

96  m 

Table  4-3  Comparison  of  Predicted  and  Observed  Errors  -  Massena 


The  observation  here  is  that  the  95%  ellipses,  for  "raw  Loran-C," 
year-round,  provide  predictions  which  exceed  the  largest  errors  actually 
observed.  This  verifies  the  contention  of  Section  3  that  these  predictions 
closely  approximate  100%  error  contours. 

The  Massena  data  can  be  presented  in  several  other  ways.  Indeed,  the 
nunber  of  ways  the  data  can  be  presented  and  analyzed  is  seemingly  limitless. 

This  last  statement  leads  to  an  aside  comment  regarding  the  importance  of 
first  obtaining  the  database.  Stated  differently,  the  authors  can  think  of 
at  least  57  different  ways  to  present  and  analyze  the  observations  made  at 
Massena  in  January  1982  (in  a  similar  study  made  about  a  decade  ago  for  the 
old  9930-2  baseline,  we  seemed  to  "make  a  living"  for  several  years  by 
providing  different  presentations  of  the  same  data.  We  call  this  the  "porno 
movie  syndrome.")  If,  several  years  from  now,  we  think  of  another  analysis 
technique,  we  can  immediately  apply  it.  There  will  never  again,  however,  be 
an  opportunity  to  collect  January  1982  data.  Thus,  the  first  priority  is  to 
obtain  the  data  base  -  however  modest  it  may  at  first  appear. 

Herein  we  will  provide  only  some  of  the  possible  presentations  (leaving 
roam  for  "further  study.")  Figures  4. 4. a.  and  4.4.b.  show  the  cross-track 
and  along-track  errors  as  a  function  of  time.  The  seasonal  nature  of  the 
variations  is  clearly  illustrated. 
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Pigure  4-4  illustrates  another  important  concept:  how  does  one  pick 
the  reference  value  for  the  plots?  ftius  far  we  have  chosen  the  reference 
point  as  the  average  value  of  all  available  data.  If  we  were,  for  example, 
only  interested  in  the  month  of  May,  we  might  examine  figure  4-4  and 
erroneously  conclude  we  are  "stuck"  with  a  seemingly  consistent  cross-track 
error  of  about  70  meters.  The  situation  therefore  is  misleading  for  if  we 
were  truly  interested  only  in  May,  we  would  have  chosen  the  average  position 
observed  (or  expected  from  previous  studies)  in  May  as  the  reference.  A 
note  of  caution  should  be  sounded  that  it  is  not  yet  clear  at  this  point  how 
this  average  value  can  be  obtained  for  all  waypoints.  It  is  claimed, 
however,  that  this  is  a  problem  that  can  be  addressed  in  a  systematic 
fashion  at  a  later  date  in  the  study.  First  we  should  attempt  to  determine 
if  consideration  of  such  details  as  this  is  warranted:  perhaps  there  is 
some  other  fundamental  limitation  which  will  rule  out  the  use  of  Loran-C  in 
the  Seaway.  Only  if  no  such  limitations  are  found  should  we  turn  our 
attention  to  such  details  as  this  which  is  simply  noted  in  passing. 

Having  said  about  all  that  need  be  said  about  "raw  Loran-C”  performance 
at  Massena,  we  should  proceed  with  an  examination  of  hew  the  data  can  be 
used  to  test  other  predictions  of  Section  3.  Before  attempting  to  construct 
differential  Loran-C  observations,  we  can  examine  the  "frequency"  components 
of  the  Massena  data  records.  Figure  4-5  shows  the  total  data  record  for  the 
9960-W  baseline  observations  at  Massena,  the  "low  frequency"  component 
(obtained  by  applying  a  2-week  sliding  average  filter  to  the  raw  data)  and 
the  residual  "shorter-term"  component.  The  plots  are  analogous  to  those  of 
figure  3-4.  TTie  same  components  for  the  9960-X  baseline  are  shown  in  Figure 
4-6. 

If  we  try  to  compare  the  statistics  of  the  residuals  plotted  in  figures 
4-5  and  4-6  with  the  20-nanosecond  "shorter-term"  variation  standard 
deviation  used  for  the  predictions  of  Section  3,  we  see  substantial 
disagreement.  The  residuals  for  the  9960-W  baseline  have  a  standard 
deviation  of  47  nanoseconds  and  the  residuals  for  the  9960-X  baseline  have  a 
standard  deviation  of  82  nanoseconds. 

Actually,  the  problem  is  only  that  we  are  mixing  models.  This  is  an 
understandable  thing  to  do  given  the  small  data  base,  and  we  will  be  tempted 
to  do  it  throughout  this  section.  We  are  now  seeing  the  effects  indicated 
in  the  note  of  section  3.5  regarding  the  classification  of  various 
components  of  the  Loran-C  time  difference  variations.  The  need  for  an 
in-depth  discusion  is  now  indicated. 
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The  20  nanosecond  figure  used  in  Section  3.5  (and  Appendix  C)  came  from 
a  differential  Lor an -C  model  and  represented  that  component  of  the 
variations  assumed  independent  from  site  to  site.  Thus  we  cannot  obtain  a 
direct  measure  of  it  by  examining  the  Massena  data  alone.  Ideally ,  we  would 
obtain  an  estimate  of  it  by  subtracting  the  Massena  data  record  from  a  data 
record  obtained  at  a  nearby  site.  Thus,  we  would  simulate  differential 
Loran-C  performance  at  the  other  site  and,  by  breaking  the  resulting  data 
record  down  into  a  "seasonal"  and  "shorter-term"  (more  properly  called 
"site-independent")  component,  could  test  the  bases  of  the  prediction 
model.  We  would  hope  to  find  the  "seasonal”  component  properly  represented 
by  the  double  range  difference  and  the  "residual"  component  with  a  standard 
deviation  of  about  20  V?  nanoseconds.  This  then  is  the  approach  we  would 
take  ^f  we  had  a  reasonable  database  from  another  nearby  site. 

Since  we  do  not  have  this  other  database,  we  will  spend  more  time 
scrutinizing  the  Massena  data.  Although  a  subopt iraal  occupation,  we  can 
obtain  valuable  insight  and  the  exercise  is  deemed  worthwhile.  Care  must 
simply  be  taken,  as  we  have  already  begun  to  emphasize,  that  we  are  not  yet 
doing  a  direct  differential  Loran-C  analysis. 

In  the  preliminary  discussion  of  the  Loran-C  time  difference  variation 
model  in  Section  3,  we  referred  to  Appendix  B  where,  almost  as  an  aside,  we 
noted  that  the  temporal  and  spatial  nature  of  the  TO  variations  were 
interrelated.  We  must  now  expand  upon  this  concept  and  begin  to  introduce 
more  detailed  consideration  of  the  variations.  Use  of  the  reasonably 
complete  Massena  data  base  will  prove  invaluable  in  illustrating  the 
concepts.  We  begin  by  recalling  Section  3.4  claimed  there  were  4  major 
frequency  components  ye  should  consider.  Having  described  the  four 
components,  we  proceeded  in  Section  3.5  with  predictions  based  on  two 
components:  one  that  is  highly  correlated  from  site  to  site  and  one  that  is 
independent  from  site  to  site.  These  two  classifications  of  components  (one 
having  4  constituents,  the  other  having  2)  are  not  directly  related.  Thus, 
we  should  call  the  four  component  classification  the  "temporal"  one,  and  the 
two  component  classification  the  "spatial"  one  to  emphasize  the  distinction. 

At  the  low  frequency  end  of  the  temporal  classification,  we  find  the 
so-called  "seasonal"  component  is  highly  correlated  from  site  to  site  and 
from  baseline  to  baseline.  Thus  it  fits  in  the  highly  correlated  portion  of 
the  spatial  classification.  We  see  significant  evidence  of  this  correlation 
by  comparing  the  plot  of  figure  4-5. b.  to  that  of  figure  4-6. b.  and  even  to 
3-4. b.  At  the  other  end  of  the  temporal  variation  classification,  we  find 
the  so-called  "near  instantaneous"  frequency  component.  It  is  claimed  that 
this  component  is  also,  to  a  large  extent,  highly  correlated  from  site  to 
site.  Small  data  record  discontinuities  resulting  from  local  timing 
corrections  (corrections  entered  by  the  SAM  that  produce  identical  effects 
throughout  the  service  area)  comprise  a  large  portion  of  this  temporal 
component.  Thus,  when  we  chose  the  terminology  "seasonal"  and 
"shorter-term"  components  in  Section  3.5,  we  were  being  somewhat 
misleading.  Indeed,  of  all  four  temporal  classification  components, 
we  expect  only  a  small  portion  of  the  so-called  "short-term"  effects 
and  the  portion  of  the  "near-instantaneous"  effects  produced  by  atmospheric 
noise  to  show  substantial  de -correlation  from  site  to  site.  The  prime 
contributors  are  the  effects  of  weather  fronts  which  pass  between  the 


sites  being  compared  along  the  signal  path.  These  effects  happen  too  slowly 
to  be  called  "near-instantaneous"  but  too  rapidly  to  be  called  "medium-term" 
or  "seasonal"  (as  defined  on  pages  3-7  and  3-8).  Thus,  for  the  remainder  of 
this  report,  we  will  be  more  detailed  in  our  modeling  and  clear  up  the 
terminology  by  using  the  terms  "spatially  correlated"  and  "spatially 
independent”  as  the  components  of  the  spatial  classification. 

The  application  of  the  above  discussion  is  that  whereas  we  expect 
differential  Lor an -C  corrections  to  remove  much  (depends  on  the  A  drd)  of 
the  "low  frequency”  error  shown  in  figures  4-5. b.  and  4-6. b.,  we  also  expect 
a  considerable  portion  of  the  "residual"  errors  shown  in  figures  4-5. c.  and 
4-6.C.  to  be  removed.  By  using  the  2-week  filter  in  figures  4-5  and  4-6,  we 
have  merely  separated  the  component  bearing  the  temporal  classification 
"seaonal"  from  the  other  components.  He  expect  true  differential  Loran-C  to 
do  better  so  that  figure  4-5. c.  and  4-6. c.  do  not  provide  a  good 
differential  Loran-C  "simulation." 

Although  we  have  just  belabored  the  point  that  manipulation  of  "Hassena 
only"  data  is  a  suboptimal  approach  to  testing  for  expected  differential 
Loran-C  performance,  there  are  some  further  steps  we  can  take  in  trying  to 
get  a  good  "simulation."  One  reason  for  doing  this  is  that  Massena  is  our 
largest,  highest  quality  data  base.  Another  is  that  we  can  use  it  to  lead 
into  the  next  section  where  we  can  illustrate  more  detailed  characteristics 
of  the  variation. 

He  proceed  by  noting  our  previous  claim  that  we  expect  the  "medium 
term"  frequency  component,  as  well  as  the  "seasonal"  frequency  component  to 
be  highly  correlated  from  site  to  site.  Thus,  we  expect  differential 
corrections  bo  remove  a  certain  portion  (again,  proportional  to  a  DRD)  of 
this  component  also.  Importantly,  just  as  with  the  seasonal  component,  we 
can  remove  this  "spatially  correlated”  component  through  time  domain 
filtering.  Specifically,  we  can  apply  a  sliding  3-day  filter  to  the  Massena 
data  record  and,  by  subtraction,  remove  both  the  seasonal  and  the  medium 
term  components  (i.e.,  the  components  claimed  to  be  almost  entirely 
spatially  correlated)  from  the  other  components.  The  process  is  illustrated 
for  the  9960-H  baseline  in  figure  4-7  and  for  the  9960-X  baseline  in  figure 
4-8. 


We  see  that  the  standard  deviations  of  the  "residuals"  have  now  been 
reduced  to  26  and  43  nanoseconds  for  the  9960-W  and  9960-X  baselines 
respectively,  this  represents  a  substantial  improvement  over  the  4?  and  82 
nanosecond  residual  standard  deviations  after  the  application  of  the  2  week 
filter  and  is  approaching  the  figure  we  expect.  The  remaining  disagreement 
is  not  an  immediate  cause  of  concern  because  we  know  that  part  of  the 
remaining  short  term  and  near-instantaneous  components  are  spatially 
correlated.  Unfortunately*  for  these  components,  the  spatial  and  temporal 
nature  of  the  variations  are  significantly  related  so  we  cannot  remove  any 
more  spatially  correlated  components  by  time  domain  filtering  alone.  Thus 
we  have  gone  about  as  far  as  we  can  go  with  "Massena  only"  data 

We  should  not  conclude  this  section  on  such  a  negative  note:  the 
Massena  data  has  proved  extremely  useful  in  verifying  much  of  what  was 
presented  in  Section  3.  For  "raw"  Loran-C,  the  predicted  statistics  agree 
extremely  well  with  those  observed  at  Massena.  Thus,  we  can  verify  the 
conclusion  of  Section  3  that  raw  Loran-C,  year-round,  is  inadequate  for  the 
St  Lawrence  Seaway  requirements.  Although  we  could  not  obtain  proper 
verification  of  the  differential  Loran-C  statistics,  this  is  a  limitation  to 
be  encountered  anytime  data  is  available  only  from  one  site.  From  the 
Massena  data,  we  saw  no  indications  that  our  predictions  are  substantially 
out  of  line. 


4.3  Evaluation  of  Data  From  Tibbetts  Point 


Our  prism  concern  right  now  is  to  find  some  other  site  with  reasonable 
data  to  further  pursue  the  verification  of  our  differential  Loran-C 
performance  predictions.  Since  we  have  a  fairly  complete  data  record  from 
Tibbetts  Point  after  early  February  1982,  we  will  start  with  that  site.  The 
extreme  range  from  Massena  makes  this  a  fortunate  site  to  have  winter/spring 
data  from. 

The  time  difference  records  from  Tibbetts  Point  are  presented  in 
figures  4-9  and  4-10  which  include  the  "raw"  data,  the  filtered  (2-week) 
data  and  the  "residuals." 
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The  first  comment  we  can  make  is  that  the  peak-to-peak  values  of  the 
seasonal  component  is  about  900  nanoseconds  for  both  baselines  -  9960-W 
showing  a  slightly  larger  swing  than  9960-X.  Using  the  sinusoidal 
approximation,  we  would  expect  "sigmas”  of  about  320  nanoseconds.  Hie 
computed  "sigmas"  turn  out  to  be  about  10%  higher  -  a  situation  we  might  not 
encounter  if  we  had  data  in  January.  For  the  most  part,  given  we  have  only 
about  4-1/2  months  of  data,  we  conclude  this  is  good  agreement  with  the 
method  used  to  predict  "seasonal  sigmas."  (We  are  fortunate  to  have  this 
particular  4-1/2  month  period.  Although  not  optimal,  it  allows  a  better 
estimate  of  the  peak-to-peak  variation  than  we  would  obtain  from  almost  any 
other  4-1/2  month  part  of  the  year.) 

We  see  a  problem  in  that  we  predicted  "sigmas"  of  242  and  325 
nanoseconds  for  9960-W  and  9960-X,  respectively,  at  Tibbetts  Point.  Whereas 
the  prediction  error  for  9960-X  is  less  than  10%,  it  is  over  50%  for  the 
9960-W  baseline.  This  observation,  in  conjunction  with  the  observation  that 
the  Massena  data  agreed  so  well  with  the  predictions  based  on  other  N.E.U.S. 
site  data,  suggests  a  significant  violation  of  the  uniform  propagation 
(modified  for  "land/seawater")  model.  This  is  a  significant  finding  (a 
"camel")  and  illustrates  the  utility  of  the  low  density  analysis.  Before 
considering  implications  of  the  problem  for  the  differential  Loran-C 
analysis,  we  should  finish  the  discussion  of  "raw"  Loran-C  performance  at 
Tibbetts  Point. 

Figure  4-11. a.  is  a  scatter  plot  of  the  fixes  obtained  by  using  the 
Tibbetts  Point  data  along  with  a  superimposed  error  ellipse  generated  by 
using  the  observed  statistics.  The  predicted  error  ellipse  is  provided  as 
figure  4-11. b.  to  facilitate  side-by-side  comparison. 
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Saw  comments  are  in  order.  First,  we  should  note  that  the  basic 
elliptical  pattern  is  generally  a  correct  representation  of  the  actual 
distribution  of  errors.  Next  we  note  that  the  actual  ellipse  is  rotated 
about  30°  relative  to  the  predicted  ellipse.  This  rotation  is  a 
consequence  of  the  incorrect  estimate  of  the  ratio  of  the  variations  for  th< 
two  tine  differences  (we  predicted  a  "sigma -W/sigma-X  ratio"  of  about  0.7S 
and  observed  one  of  about  1.07).  Interestingly,  the  rotation  (along  with 
the  effects  of  a  much  higher  correlation  coefficient  than  our  conservativ 
estimate)  served  the  purpose  of  reducing  the  maximum  cross-track  error. 

We  should  now  move  to  a  consideration  of  differential  Loran-C 
performance  at  Tibbetts  Point  with  corrections  being  applied  from  a  monitor 
at  Massena.  We  begin  by  subtracting  the  Massena  data  records  from  the 
corresponding  records  at  Tibbetts  Point.  The  results  are  shown  in  figures 
4-12  and  4-13.  Again,  we  provide  the  "raw"  data,  a  filtered  (2-week) 
version  of  the  data,  and  the  residuals  of  the  filtering  process. 

We  see  that  the  peak -to -peak  value  of  the  seasonal  component  is  about 
600  nanoseconds  for  9960-W  and  about  270  nanoseconds  for  9960-X.  By  our 
sinusoidal  approximation,  we  would  thus  estimate  seasonal  standard 
deviations  of  212  and  95  nanoseconds,  respectively.  Itiese  are  "close”  to 
the  computed  values  of  238  and  112.  As  with  the  comments  on  raw  Loran-C  at 
Tibbetts  Point,  we  would  expect  better  agreement  if  we  had  data  from  January. 

The  problem  already  noted  in  the  discussion  of  the  raw  data  at  Tibbetts 
Point  becomes  even  more  pronounced  when  we  compare  these  observations  with 
our  predictions  of  the  differential  variations.  We  predicted  a  *sigma-W"  of 
101  nanoseconds  and  a  "sigma-X"  of  76  nanoseconds.  Thus  we  have  a  136% 
error  for  9960-W  and  a  47%  error  for  9960-X. 
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Breakdown  of  9960-X  Differential  Loran-C  Frequency 
Components:  Tibbetts  Point  Minus  Massena  Using  15-Day 
Smoothing 


Figure  4-14  Comparison  of  Predicted  and  Observed  Differential  Ix>ran-C 

Error  Ellipses:  Tibbetts  Point  Minus  Massena 

We  must  recall,  however,  that  the  prime  use  we  need  to  make  oi  the 
Tibbetts  Point  data  is  to  determine  how  well  we  can  predict  what  happens 
between  Tibbetts  Point  and  Massena.  In  this  regard,  we  recall  the  Section  \ 
conclusion  that  it  appears  we  need  at  least  one  more  differential  correction 
station  -  most  likely  one  in  the  vicinity  of  Wellesley  Island.  We  expect 
the  establishment  of  such  a  site  (or,  perhaps,  additional  sites)  will  r<*duo. 
a  considerable  portion  of  the  seasonal  variation  shown  in  fiqures  4-12. b. 
and  4-13. b. 


If  we  imagine,  for  the  moment,  that  we  can  remove  most  of  the  seasonal 
component  from  the  Tibbetts  Point  data  by  the  application  of  differential 
corrections  from  some  "yet-to-be-determined"  site,  we  can  focus  our 
attention  on  the  "residuals"  shown  in  figures  4-12. c.  and  4-13. c.  We  se<* 
standard  deviations  of  69  and  33  nanoseconds  for  996C-W  and  9960-X, 
respectively.  Although  roles  have  been  reversed  as  t which  TD  varies  mor°, 
we  note  we  are  left  with  standard  deviations  smaller  than  those  for  the 
"analogous"  "Massena  only"  situation  shown  in  figures  4-5. c.  and  4-6. r. 

This  is  in  spite  of  the  fact  that  Tibbetts  Point  is  as  far  southwest  of 
Massena  as  one  can  go  and  still  be  in  the  Seaway.  This  illustrates  the 
claims  of  the  previous  section  that  much  of  the  variation  remaining  a* t<~- 


the  seasonal  frequency  component  has  been  removed  is  still  spatially 
corrc'idied .  Unlike  with  the  "Massena  only"  situation,  we  now  have  a  means 
to  "get  at"  this  other  dimension  of  the  variations.  A  very  important 
observation  is  that  the  33  nanosecond  standard  deviation  for  9960-X 
approaches  the  20  (i.e.,  28)  nanosecond  standard  deviation  we  used  for 

the  "spatially  independent"  component  in  our  predictions. 

At  this  point  we  should  focus  our  attention  on  this  33  nanosecond 
statistic  for  the  9960-X  baseline  and  review  how  it  was  obtained.  First  we 
subtracted  the  Massena  data  record  from  the  Tibbetts  Point  data  record  to 
simulate  differential  Loran-C.  In  doing  so,  we  fully  expected  the  result  * 
have  all  four  frequency  compixients  remaining.  The  magnitude  of  the 
components  was  expected  to  be  in  accordance  with  the  following  breakdown: 

a.  Seasonal  frequency  component.  We  expected  this  component  of  the 
Tibbetts  Point  and  Massena  data  records  to  have  almost  perfect  spatial 
correlation.  Note  correlation  does  not  imply  a  one-to-one  relationship 
(i.e.,  when  the  Massena  data  slowly  "rises,"  we  expect  to  find  the  Tibbettc 
Point  data  slowly  rising  -  but  not  necessarily  by  the  same  amount).  mhus , 
we  expect  the  differential  corrections  will  remove  only  a  certain  percentao 
of  this  component  of  the  variation  -  that  dictated  by  Adrd  considerations. 
Thus  we  expect  bo  find  a  seasonal  frequency  component  remnant  in  the 
"Tibbetts  minus  Massena"  data. 

b.  Medium-Term  frequency  component.  As  with  the  seasonal  frequency 
component,  we  expect  almost  perfect  spatial  correlation  in  this  component  o 
the  Tibbetts  Point  and  Massena  variations.  Also,  we  expect  to  be  able  to 
remove  only  the  percentage  of  this  component  dictated  by  the  Adrd.  Thus 
this  component  is  expected  to  be  found  in  the  "Tibbetts  minus  Masr.-  '-i’  data 

c.  Short-term  frequency  component.  We  expect  this  component  •>  h a v < 
two  significant  sub-components: 

(1)  A  spatially  correlated  part  which  occurs  for  the  same  reason 
the  seasonal  and  medium  term  components  occurred.  Again,  the  percentage 
removed  will  be  dictated  by  ADRD  so  the  "Tiboe  ts  minus  Massena”  data  wi'l 
contain  some  of  this  component. 

(2)  A  spatially  independent  part  due  to  the  passage  of  fronts, 
etc.  We  do  not  expect  to  be  able  to  remove  any  of  this  so  it  will  be  found 
in  toto ,  in  the  "Tibbetts  minus  Massena"  data. 

d.  Near-instantaneous  frequency  component.  We  expect  this  component 
to  have  thr eg  significant  sub -components : 

(1)  A  spatially  correlated  part  due  to  chain  timing  correct  ions , 
most  of  the  effects  of  chain  equipment  changes,  equipment  problems  at  the 
SAM,  and  rapid  weather  effects  that  are  localized  to  the  vicinity  of  SAM. 
etc.  We  expect  this  sub-componen t  to  be  common  to  all  sites  and  thus  no*-  h 
present  in  the  "Tibbetts  minus  Massena"  data. 

(2)  A  spatially  correlated  part  that,  like  the  seasonal  and 
medium  term  components,  is  related  to  double  range  difference.  This  will  1 
present  in  the  "Tibbetts  minus  Massena”  data  because  only  the  percentage 


dictated  by  ADRD  considerations  is  removed  through  differential  corrections. 

(3)  A  spatially  independent  sub-component  due  to  local  weather 
effects,  local  equipment  reactions  to  chain  equipment  changes,  and  local 
noise,  etc.  This  cannot  be  removed  or  reduced  through  differential 
corrections  and  will  thus  be  present,  in  toto,  in  the  "Tibbetts  minus 
Massena"  data. 

By  the  time  we  had  looked  at  only  the  "raw"  Loran-C  data  at  Tibbetts 
Point  we  knew  we  had  encountered  an  anomaly  in  the  double  range  difference 
which  would  impact  on  the  worthiness  of  our  predictions.  The  anomaly  became 
even  more  noticeable  when  we  actually  accomplished  the  data  record 
subtraction  to  simulate  the  effects  of  differential  corrections  from  Massena 
being  applied  to  a  "site  of  interest"  at  Tibbetts  Point.  Since  we  had 
already  concluded  in  Section  3  that  we  probably  could  not  "make  it"  with 
only  the  monitor  at  Massena,  this  was  not  an  immediate  cause  for  alarm. 

Since  the  seasonal  frequency  component  is  hypothesized  (based  on 
considerable  corroborating  observations  at  other  sites)  to  be  almost 
perfectly  spatially  correlated,  a  comparison  of  this  component  in  the 
"Tibbetts  minus  Massena"  data  to  the  predictions  gave  a  measure  of  the  error 
in  the  assumed  ADRD.  We  assume  we  can  remove  this,  and  any  other  spatially 
correlated  component  or  sub-component  to  any  extent  we  find  necessary 
(actually,  at  Tibbetts  Point,  because  of  the  wide  channel,  no  improvement  is 
necessary)  by  appropriate  selection  of  a  closer  differential  correction 
station . 

Putting  this  "assumed  solvable"  issue  aside  for  the  moment,  we  would 
like  to  attempt  to  "get  at"  an  observed  measure  of  the  spatially  independent 
component.  The  criticality  of  this  quest  is  seen  if  we  recall  the  final 
discussions  in  Section  3.5:  because  of  geometry,  it  is  the  spatially 
independent  component  -  the  one  that  we  cannot  remove  no  matter  where  we 
locate  the  differential  station  -  that  appears  to  be  the  limiting  factor  in 
"tight  reaches"  such  as  #61  and  #62.  With  our  assumed  statistics  for  the 
spatially  independent  component,  we  were  already  "too  tight  for  comfort"  in 
these  reaches  so  that  a  critical  goal  of  the  data  analysis  should  be  to  seek 
ways  to  confirm  or  refute  these  assumptions.  If  they  are  correct,  or 
optimistic,  we  have  severe  problems.  If  they  are  pessimistic,  we  may  be  in 
good  shape. 

Thus,  we  removed  the  seasonal  component  in  the  "Tibbetts  minus  Massena" 
data  in  our  first  attempt  to  get  at  the  "spatially  independent  only"  portion 
of  the  variation.  This,  then,  is  how  we  arrived  at  the  33  nanosecond  fiqure 
for  9960-X  (and  the  69  nanosecond  figure  for  9960-W) . 

This  review  complete,  we  can  appreciate  the  question:  are  the  33/69 
nanosecond  figures  representative  of  the  spatially  independent  components  - 
i.e.,  the  best  we  can  expect?  From  the  discussion  we  can  see  we  expect  to 
Find  (if  we  somehow  could)  that  the  "real  bottom  line"  is  possibly  better 
than  the  33/69  nanosecond  figures  would  indicate  -  but  we  cannot  yet  say  how 
much  better.  We  say  this  because  sub-components  decribed  above  to  be  "of 
*•  -  pe"  b.,  c.(l),  and  d.(2)  will  be  reduced  as  we  move  the  differential 
control  station  closer  to  Tibbetts  Point.  By  use  of  a  time  domain  filter  we 
cannot  "get  at"  sub-component  c.(l)  without  removing  sub-component  c.(2) 
'erroneously),  nor  can  we  qet  at  sub-component  d.(2).  We  can,  however, 
'emove  sub-component  b.  by  use  of  a  3-day  filter.  We  do  this  in  figures 
4-15  and  4-16. 
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Figure  4-16 


Breakdown  of  9960-X  Differential  Loran-C  Frequency 
Components:  Tibbetts  Point  Minus  Massena  Usinq  3-Day 
Smoothing 


We  see  the  significant  effect  this  operation  has  had:  the  standard 
deviation  for  9960-W  has  been  reduced  to  33  nanoseconds  and  the  standard 
deviation  for  9960-X  has  been  reduced  to  20  nanoseconds.  Thus,  we  have 
reason  for  hope.  Before  considering  the  implications  of  these  results  we 
should  pause  to  compare  these  results  with  the  analogous  situation  for  the 
"Massena  only"  analysis  reflected  in  figures  4-7  and  4-8.  Note  that  by 
subtracting  the  "3-day-filtered"  component  in  both  cases  we  removed  the 
seasonal  and  medium  term  components.  The  difference  in  the  plots  (which 
resulted  in  a  standard  deviation  reduction  from  43  to  33  nanoseconds  for 
9960-W  and  from  26  to  20  nanoseconds  for  9960-X)  was  obtained  through 
"differential  action."  The  difference  represents  how  much  of  the  "residual" 
variation  was  common  to  both  sites  (and  thus  removed  in  toto)  and  how  much 
was  related  to  ADRD  (and  thus  removed  in  the  percentage  indicated  by  the 
ADRD)  . 


(Note:  A  direct  comparison  of  the  43  nanosecond  figure  to  the  33  nanosecond 
figure  does  indeed  confirm  improvement.  Likewise  a  comparison  of  the  26  and 
20  nanosecond  figures.  There  is  more  going  on  here,  however,  than  the 
direct  comparison  indicates.  Recall  the  formulation  of  Section  3  (and 
Appendix  C)  shows  that  when  we  implement  differential  Loran-C  when  there  is 
no  spatially  correlated  component,  we  actually  expect  the  standard  deviation 
to  increase  -  by  a  factor  of  JT.  Thus,  our  predictions  had  accounted  for  an 
increase  in  the  standard  deviation  if  there  were  no  spatially  correlated 
components.  In  view  of  this  consideration,  the  reductions  are  even  more 
significant  in  their  implications.) 


"Tantalizingly, "  we  cannot  distinguish  between  the  two  sources  of 
improvement  at  present.  At  worst,  all  the  improvement  was  due  to  the 
removal  of  common  error  terms.  More  likely,  some  of  the  improvement  was  du*5 
to  the  removal  of  the  common  error  terms  and  some  was  removed  because 
Massena  is  closer  to  Tibbetts  Point  than  Sandy  Hook,  N.J.  or  Cape  Elizabeth, 
Me,  are.  In  this  latter,  more  likely  case,  we  can  expect  even  further 
reduction  in  the  residual  variations  if  we  locate  a  differential  correction 
station  closer  to  Tibbetts  Point  (which,  it  should  be  clear  by  now,  we  must 
do)  . 

Bearing  this  "hoped  for"  further  improvement  in  mind  (besides  recalling 
we  are  temporarily  ignoring  the  question  of  how  to  get  rid  of  the  majority 
of  the  large  spatially  correlated  components) ,  we  should  turn  to  the 
all-important  consideration  of  how  the  TD  variations  in  figures  4-15. c.  and 
4-16. c.  translate  to  positional  errors.  Those  TD  records  are  converted  to 
position  variations  in  the  scatter  plot  of  figure  4-17. 
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Figure  4-17  Error  Ellipse  for  Differential  Loran-C  High  Frequent 
Components:  Tibbetts  Point  Minus  Massena 

From  the  plot  we  see  further  cause  for  optimism  -  a  cause  that  bas  brer 
evident  in  previous  scatter  plots.  The  increased  "sigma-W/sigma-X  •'at’e 
has  again  rotated  the  error  ellipse.  This  is  of  absolutely  no  signif  arv 
in  the  "wide  open"  channel  near  Tibbetts  Point  but  is  of  immense  concern 
the  tighter  reaches  such  as  #61  and  #62  where,  as  indicated  ’n  Sectior  we 
were  stymied  by  an  almost  "worst  case"  orientation  of  the  major  (larger) 
axis  of  the  ellipse.  Although  the  data  was  not  obtained  in  the  vicinity  of 
these  reaches,  we  can  simulate  what  we  now  have  reason  to  expect  will  occur 
in  these  reaches  -  with  a  nearby  differential  control  station  -  by  applying 
the  geometry  at  these  reaches  to  the  Tibbetts  Point  data  of  figures  4-1 5.<~, 
and  4-16. c.  This  is  not  a  "reckless"  thing  to  do  since,  because  these 
reaches  are  closer  to  Massena  than  Tibbetts  Point  is,  we  expect  even  sma’W 
variations  would  be  observed  at  these  sites.  This  simulation  is 
accomplished  in  figure  4-18. 
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Figure  4-18  Application  of  the  Loran-C  Data  of  Figure  4-17  to  the 
Geometry  at  Reach  162 
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Again,  we  see  the  "benign"  change  in  the  orientation  of  the  ellipse. 

The  significant  reduction  in  the  variation  of  the  9960-X  TD  is  responsible 
for  this.  The  phenomenon  is  illustrated,  somewhat,  in  figure  4-19  where  we 
show  the  actual  LOP's  at  Reach  #61.  We  see  immediately,  of  course,  that  the 
crossing  angles  of  the  LOP's  are  far  from  optimal.  More  importantly,  we  see 
that  the  W  LOP  is  nearly  perpendicular  to  the  course.  Thus,  variations  in 
the  9960-W  TD's  produce  substantial  along-track  errors  but  not  much 
cross-track  error.  The  9960-X  LOP  is  by  no  means  parallel  to  the  course  but 
it  does  have  a  larger  "cross-track-error-axis-projection"  than  the  9960-W 
LOP.  Thus,  a  significant  reduction  in  the  9960-X  variations  from  what  we 
predicted  (i.e.,  a  reduction  from  28  nanoseconds  to,  at  most,  20 
nanoseconds)  has  a  significant  effect  on  the  cross-track  error  near  reaches 
#61  and  62. 


It  is  Interesting  at  this  point  to  note  how  the  95%  percent  contour  is 
beginning  to  truly  approach  that  gaussian  characteristic  -  a  small  but 
noticeable  percentage  of  the  fixes  of  figure  4-18  are  outside  the  contour 
(the  data  record  shows  15  out  of  just  over  200  points  are  outside  the 
contour).  In  figure  4-20,  we  show  the  same  scatter  plot  as  in  figure  4-18 
but  now  superimpose  the  99.9%  probability  contour. 
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Figure  4-20  Same  Data  Points  as  in  Figure  4-18  With  99.9%  Probability 
Error  Ellipse 


Here  we  see  direct  evidence  of  applicability  of  the  gaussian 
characteristics.  We  note  that  the  one  point  which  is  outside  the  ellipse  is 
"way  out."  This  fix  is  a  result  of  the  data  which  shows  up  as  a  "spike"  in 
figures  4-15. c.  and  4-16. c.  just  before  1  April  1982  (Julian  Day  91).  The 
first  thing  we  should  say  is  that  the  large  TD  variations  did  not  result  in 
an  abnormally  large  cross-track  error  so  we  are  interested  in  the  "anomaly" 
only  in  a  statistical  sense  (i.e.,  "this  time"  the  direction  of  the  large 
positional  error  was  benign  but  what  about  next  year?).  This  anomalous 
ooint  will  be  discussed  in  greater  detail  in  the  next  section.  For  now,  we 
wi':l  simply  point  out  that  the  "wild  behavior"  at  about  1  April  1982  is 
evident  in  all  data  records  presented  thus  far  -  i.e.,  in  both  the  9960-W 
and  9960-X  data  records  for  both  Massena  and  Tibbetts  Point  (check  figures 
4-5. c.,  4-6. c.,  4-9. c.,  and  4-10. c.).  This  indicates  a  portion  of  the 
anomaly  is  spatially  related  (specifically,  of  type  c . ( 1 ) )  and  would  have 
been  further  reduced  by  a  closer  differential  correction  station. 

Before  finishing  up  the  discussion  of  the  "Tibbetts  minus  Massena"  data 
wo  should  return  to  the  question  of  the  large  seasonal  and  medium  term 
f-'-'quency  component  that  we  temporarily  put  aside.  Even  though  we  do  not 
vo*  know  the  "bottom  line"  on  the  spatially  independent  variations  (though 
wo  >"iow  we  must  be  approaching  it)  we  see  from  figure  4-20  that  we  have  some 
*>rcMthing  room"  even  after  we  allow  for  17  meter  vessel  half-widths  and  as 
murh  as  10  meters  for  guidance  error.  Thus  we  can  allow  for  "mild" 
spatially  correlated  components  to  be  added  to  the  spatially  Independent 
i  omponents  -  particularly  for  the  9960-W  Daseline  In  the  vicinity  of  Reaches 
*6’  and  62  -  and  still  be  safe.  Thus,  we  could  re-do  the  differential 
Loran-C  predictions  of  Section  3  based  on  the  new  evidence  about  the  size  of 
the  spatially  independent  components.  We  claim,  however,  that  such  an 
update  is  premature  at  present. 
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If  we  were  to  attempt  an  update,  we  trould  have  to  revise  our  estimation 
of  the  spatially  correlated  components  as  well  as  the  spatially  independent 
ones  since  available  data  also  indicates  the  need.  Given  only  two  data 
sites  however,  our  only  recourse  would  be  the  application  of  another  linear 
model  -  i.e.,  create  straight  lines  of  different  slopes  than  those  shown  in 
figure  3-5:  the  basis  for  the  predictions  of  Section  3.  Hie  updated  slopes 
would  have  to  be  determined  by  forcing  straight  lines  to  pass  through  the 
Tibbetts  Point  and  Massena  data  points. 

Although  at  first  glance  this  might  seem  like  a  very  "artificial"  thing 
to  do  it  would  be  consistent  with  the  moet  popular,  classical  approach  to  a 
non-linear  problem.  In  the  approach  it  is  argued  that,  although  a  generally 
nan-linear  effect  is  conceded,  a  linear  approximation  to  the  variation  is 
satisfactory  for  a  small  area  about  the  operating  region.  The  arguments  are 
based  on  Taylor  series  considerations.  If  we  were  to  use  them,  we  would,  in 
essence,  be  claiming  that  terms  of  higher  order  than  the  linear  ones  (there 
may  be  several  independent  variables)  are  negligibly  small.  For  the  case  in 
point,  we  would  be  conceding  that  a  simple  linear  relation  does  not  hold  all 
the  way  from  the  "origins"  (Sandy  Hook/Cape  Elizabeth)  to  Tibbetts  Point, 
but  is  very  nearly  what  happens  over  the  small  area  of  the  Seaway.  Thus,  we 
would  be  arguing  the  "bulk"  of  the  non-linear  effects  have  occurred  before 
we  get  to  the  Seaway.  Herein  lies  the  reason  we  cannot  use  this  approach: 
we  already  know,  from  section  4.2,  that  the  Massena  data  agrees  very  well 
with  what  a  simple  linear  extension  of  N.E.U.S.  coast  data  suggests. 

Thus,  we  know  that  if  Tibbetts  Point  is  "way  off  the  curve,"  the  bulk 
of  the  non-linear  effects  are  occurring  in  the  Seaway.  Functionally,  this 
means  we  have  evidence  that  our  relationships  have  significant  higher  order 
spatial  derivatives  -  which  translate  to  significant  higher  order  Taylor 
series  coefficients  -  in  the  Seaway.  Hie  practical  implication  is  in 
regards  to  number  of  monitor  stations  we  need  and  our  estimates  as  to  the 
safety  margin  the  resulting  Bystem  provides.  To  Bee  this,  suppose  we  use  a 
linear  approximation  between  Massena  and  Tibbetts  Point.  The  assumed 
propagation  velocity  changes  (the  slopes  determine  by  the  regression 
analyses)  would  be  greater  than  those  used  in  Section  3.  This  would  make 
our  predictions  about  how  far  we  can  move  away  from  the  differential 
correction  station  without  encountering  unacceptably  large  errors  less 
optimistic.  This  is  a  step  in  the  right  direction. 

However,  we  would  still  be  exhibiting  too  much  optimism.  Suppose  that 
the  propagation  conditions  that  indicate  the  linear  regreesion  analysis 
slope  that  apparently  holds  true  from  the  origins  to  Massena  also  hold  true 
considerably  further  southwest  along  the  Seaway.  Worst  case,  it  would  hold 
true  all  the  way  to  Reach  #61.  In  this  case  we  would  need  an  even  steeper 
slope  to  account  for  the  variations  from  Reach  #61  to  Tibbetts  Point.  Thus 
our  results,  were  we  to  revise  the  predictions  of  Section  3  based  only  on 
the  data  available  to  us  now,  would  be  much  too  optimistic.  In  summary, 
therefore,  we  desperately  need  information  from  other  sites  -  ideally  in  the 
region  between  Wellesley  Island  and  Reach  #61. 
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In  the  following  section,  we  will  attempt  to  examine  data  from  other 
sites.  Before  doing  so,  we  should  review  the  conclusions  of  this  section. 
With  data  from  Tibbetts  Point  as  well  as  Massena,  along  with  a  substantial 
amount  of  data  processing,  we  have  reason  to  challenge  the  underlying 
assumptions  of  the  predictions  of  Section  3  -  for  the  southwest  portion  of 
the  Seaway.  We  have  concluded  we  were  too  optimistic  about  the  spatially 
correlated  component  of  the  variations  but  too  pessimistic  about  the 
spatially  independent  component.  Given  the  discussion  at  the  end  of  Section 
3.5  regarding  the  implications  of  the  assumed  magnitude  of  the  spatially 
independent  component  at  Reaches  #61  and  62,  this  type  of  finding  is  the 
"next  best  thing"  to  finding  we  were  pessimistic  in  all  regards.  There  is 
no  (Loran-C)  solution  if  the  spatially  independent  component  is  too  large. 
The  spatially  correlated  component  can  be  reduced  by  proper  selection  of  the 
differential  control  station (s) . 

The  concern  now  is  whether  or  not  we  have  adequate  data  from  other 
sites  to  proceed  with  the  selection  process. 


4.4  Evaluation  of  Data  from  Other  Sites. 


Examination  of  Table  4-1  does  not  give  us  much  reason  for  hope  in  our 
further  investigations.  As  it  turns  out,  however,  we  are  partially  in 
Luck.  We  have  a  very  nearly  complete  data  record  from  the  Ad  Hoc  site  at 
Iroquois  Lock  for  the  period  from  Julian  Day  80  to  Julian  Day  108  (21  March 
to  18  April).  Coincidentially,  this  is  one  of  the  two  periods  we  are  most 
concerned  about  -  spanning  a  period  which  includes  the  opening  of  the 
shipping  season.  This,  however,  is  not  the  real  reason  for  the  luck.  The 
real  reason  can  be  seen  by  examining  plots  of  the  data  records  for  this 
period  from  Massena  and  Tibbetts  Point.  Hie  records  are  provided  in  figures 
4-21  and  4-22. 
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Figure  4-21 


Massena  9960- W/X  Data  for  Jul ‘ 


Days  80  Through  108 
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Figure  4-22  Tibbetts  Point  9960-W/X  Data  for  Julian  Days  80 
Through  108 


We  see  variations  that  are  almost  as  large  as  the  total  variation  seen 
throughout  the  year.  This  is  not  too  surprising  since  we  are  concentrating 
on  the  "h’-  h  slope"  portion  of  the  seasonal  component  waveform,  i.e.,  we  are 
viewing  the  dramatic  change  from  winter  to  summer.  The  comparison  is 
illustrated  in  table  4-4. 


Site /Baseline 


"Seasonal  Pk-Pk" 


"30-Day  Pk-Pk" 


Massena  -  w  350  nsec 

Massena  -  X  750 

Tibbetts  Point  -  W  900 

Tibbetts  Point  -  X  850 


350  nsec 
590 
850 
700 


Table  4-4  Tabulation  of  30-Day  and  Seasonal  Peak-to-Peak  Variations 

Before  proceeding,  a  few  comments  are  in  order.  Ttie  "seasonal  pk-pk" 
readings  were  obtained  by  examining  the  "filtered  out"  seasonal  component. 
Thus,  they  were  somewhat  smaller  than  the  total  variation  in  the  data 
record.  The  "30-day  pk-pk"  reading  represents  true  "max  minus  min" 
readings.  With  such  a  short  period,  we  do  not  have  the  luxury  of  examining 
all  frequency  components  (e.g.,  a  2-week  average  is  meaningless  for  the 
first  or  last  week  of  the  period.  This  makes  a  large  impact  on  a  30-day 
data  record).  TTiis  illustrates  the  thought  that  we  would  prefer  to  have 
more  data  from  other  sites.  Lacking  this,  we  have  the  "next  best  thing." 
Given  a  choioe  of  any  30-day  period,  March/April  is  preferable  to  other  less 
eventful  periods.  An  auxiliary  benefit  is  that  the  shipping  season  begins 
on  about  1  April  so  that  our  data  record  straddles  the  period  of  most 
critical  concern. 
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Figure  4-23  Iroquois  Lock  9960-W/X  Data  for  Julian  Days  80 

Through  108 

Figure  4-23  shows  the  data  of  the  period  at  Iroquois  Lock . 

We  see  strong  correlation  with  corresponding  plots  of  figures  4-?1 
-.-22.  Considering  the  "size"  of  the  variations,  we  see  much  stronger 
agreement  with  the  Massena  data  record  than  with  the  Tibbetts  Point  r* 
suggesting  Iroquois  Lock  is  closer  to  Massena  in  a  "double-range-di  f fe*- 
sense  "  Comparing  the  "sigmas,"  it  appears  the  "closeness"  of  Troquo'  ■ 
to  Massena  is  more  pronounced  for  the  99h0-x  baseline. 


From  Table  C-2  we  see  Iroquois  Lock  (Reach  #48),  relative  to  Massen 
has  a  A  DRD  of  about  75  km  for  9960-W  and  about  50  km  for  9960-X.  Tib*', 
^>oint,  relative  to  Massena,  has  a  ,A  DRD  of  about  320  km  for  9960-W  and  i 
170  km  for  9960-X.  Thus,  we  see  why  Iroquois  Lock  "looks  more  like"  Ma^ 
than  Tibbetts  Point,  and  why  Iroquois  appears  closer  to  Massena  fo*-  *-he 
9960-X  baseline  than  for  the  9960-W  baseline.  If  we  note  the  peak-to-pe 
variations  at  Iroquois  Lock  we  obtain  the  values  of  Table  4-^ 
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Recall  we  finished  the  last  section  with  two  unanswered  questions:  "how 
small  is  the  spatially  independent  component  of  the  variation?"  and  "where 
does  the  non-linearity  in  the  spatially  correlated  component  occur?"  We  are 
not  yet  able  to  address  the  first  question  but  we  can  examine  the  second. 

The  question  affected  the  selection  of  tiro  choices.  The  first  choice  was 
that  a  "local  linear"  approximation  describes  what  happens  between  Massena 
and  Tibbetts  Point.  The  second  is  that  there  is  some  other  site,  closer  to 
Tibbetts  Point,  which  we  should  "pair  with"  Tibbetts  Point  in  defining  the 
slope  of  the  straight  line  assumed  to  describe  conditions  in  the  "tiant" 
reaches  from  #61  to  #74.  Adding  these  new  observations,  we  can  address  tv 
"choice"  question  of  Section  4.3  and  boldly  state:  "it  depends." 

Actually,  the  plots  of  figure  4-24  perfectly  illustrate  the  concepts 
involved  in  the  choices.  In  figure  4-24. a.,  we  see  a  straight  line 
approximation  between  Massena  and  Tibbetts  Point  looks  valid:  for  the 
9960-W  baseline,  the  Iroquois  Lock  data  plots  almost  right  on  the  line. 
Conversely,  the  Iroquois  data  for  9960-X  suggests,  given  we  are  concerned 
most  about  what  happens  between  Iroquois  and  Tibbetts  Point,  we  should  draw 
a  line  of  "steeper  slope"  through  the  Iroquois  and  TibbettB  Point  data. 

But,  suppose  data  from  a  site  closer  to  Tibbetts  Point,  say  Prescott,  showed 
essentially  the  same  variation  as  at  Iroquois?  In  such  a  case,  the  line 
through  the  Tibbetts  Point  and  Prescott  data  points  would  have  a  still 
steeper  slope.  Unfortunately,  we  find  we  have  no  other  substantial  data 
base  to  draw  upon.  Thus,  we  must  make  our  final  statements  about  the 
spatially  correlated  portion  of  the  data  at  this  point. 

TTie  question  of  the  source  of  the  non-linearity  persists.  It  is 
interesting  to  note  that  in  Section  4.3,  we  were  worried  about  the  9960-W 
baseline  in  this  regard.  Given  observations  at  only  two  sites  in  the  area 
of  interest,  it  was  impossible  to  test  for  linear  behavior.  The  only  other 
source  of  information  came  from  observations  outside  the  area  of  interest. 
Since  those  observations  tended  to  agree  with  observations  at  Massena,  we 
feared  a  non-linear  effect  starting  at  some  unknown  location  between  Massena 
and  Tibbetts  Point.  The  criticality  of  some  of  the  reaches  between  these 
sites  was  the  reason  for  the  fear. 

With  some  data  from  Iroquois  Lock,  we  at  last  can  concentrate  on  the 
"operating  region."  It  appears  we  can  feel  somewhat  confident  about  a 
linear  description  ol  the  variations  between  Massena  and  Tibbetts  Point  for 
9960-W.  We  are  tentative  about  these  statements  first  because  of  the 
'imited  amount  of  data  from  Iroquois  Lock  (although  the  fact  that  it 
encompasses  the  meet  questionable  part  of  the  shipping  season  is  a 
oenefit).  Another  reason  is  that  Iroquois  Lock  is  a  bit  too  far  away  from 
the  area  of  critical  concern  for  our  purposes. 

Regarding  the  9960-X  baseline,  the  Iroquois  data  seems  (tentativeness 
works  both  ways)  to  confirm  our  fears.  In  view  of  the  discussions  at  the 
end  of  Section  4.3,  tiiis  is  particularly  disturbing  -  the  9960-X  baseline  i<-- 
the  larger  contributor  to  cross-track  error.  Thus,  we  have  still  further 
reason  bo  obtain  measurements  near  Wellesley  Island. 

Before  turning  to  an  examination  of  what  the  Iroquois  Lock  data  tells 
us  about  the  spatially  independent  component,  a  final  word  is  in  order. 
Recall  in  Section  4.3  we  concluded  we  were  finding  the  "next  best  thing." 

Our  predictions  regarding  the  spatially  independent  component  were  proving 
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pessimistic.  This  is  a  joyful  finding:  there  is  nothing  we  can  do  to 
improve  upon  this  component.  The  spatially  correlated  component  was  larger 
than  we  felt  would  be  the  case.  This  situation,  fortunately,  is  "fixable." 
Thus,  we  know  a  solution  exists,  we  are  simply  seeking  a  path  to  the  optimal 
one.  We  have  not  yet  found  sufficient  evidence  to  state  what  the  answer  is 
but  are  getting  closer.  Note  we  should  avoid  the  temptation  to  simply  "plan 
on  two  or  three  differential  control  stations  between  reaches  61  and  74." 

If  one  site  (in  addition  to  Massena)  can  be  strategically  located  to  solve 
the  problem,  this  is  the  solution  we  should  seek.  The  simplest  solution 
typically  proves  to  be  the  most  reliable  and  the  most  successful. 

We  begin  consideration  of  what  the  Iroquois  data  can  tell  us  about  the 
spatially  independent  component  by  generating  the  "Iroquois  minus  Massena" 
data  plots  shown  in  figures  4-25  and  4-26. 
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Figure  4-25 


Breakdown  of  9960-W  Differential  TD  Frequency  Components: 
Iroquois  Lock  Minus  Massena  Using  3-Day  Smoothing 
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Figure  4-26  Breakdown  of  9960-X  Differential  TD  Frequency  Components: 

Iroquois  Lock  Minus  Massena  Using  3-Day  Smoothing 


Several  comments  immediately  come  to  mind  regarding  the  "raw  data" 
plots  in  figures  4-25. a.  and  4-26. a.  First,  we  note  that  the  9960-X  "sigma" 
has  been  reduced  to  20  nanoseconds  -  even  without  having  to  remove  the 
"3-day  filtered*  component.  Thus,  we  have  confirmation  of  our  observation 
in  Section  4.3  that,  for  the  9960-X  baseline,  the  spatially  independent 
component  is  smaller  than  we  predicted  in  Section  3.  In  examining  figure 
4-26, a.  it  does  not  appear  there  is  much  room  for  improvement  except  for  cn<r 
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period  of  time  surrounding  Julian  Day  87.  The  "sigma"  of  the  data  in  figure 
4-26. c.  can  be  reasonably  concluded  to  be  an  excellent  stab  at  the  "bottom 
line  estimate"  for  the  standard  deviation  of  the  spatially  independent 
differential  Loran-C  "remnant."  Assigning  a  value  of  11  nanoseconds  to  this 
component  "before  differential  corrections"  gives  us  the  desired  result 

11  /T  *  16. 


X 

v. 


In  considering  the  9960-W  data,  we  note  that  the  plot  of  figure 
4-25. a.,  the  differential  "residual"  bears  a  striking  resemblance  to  the 
original  data  records  for  9960-W  at  Iroquois  Lock  and  Massena.  This  both 
illustrates  and  confirms  the  contention  of  Section  4.3  regarding  the 
"remnant”  spatially  correlated  components  after  the  differential 
corrections.  These  components  are  only  removed  to  a  certain  percentage 
(depends  on  ADRD  -  although  we  are  not  yet  sure  how  to  compute  A  DRD  for 
both  baselines).  Noting  the  "sigma"  in  figure  4-25. c.,  we  see  the  standard 
deviation  for  the  spatially  independent  component  of  9960-W  baseline  finally 
agrees  with  that  of  the  9960-X  baseline. 


It  is  important  to  note  that  the  a1®0  from  Reach  #61  to  Reach  #74  is 
81  km  for  9960-W  and  40  km  for  9960-X  (see  Table  C-3).  These  reaches  span 
the  "critical  area  of  concern"  identified  previously.  We  note  that  the 
ADRD  between  Massena  and  Iroquois  Lock  is  75  km  for  9960-W  and  50  km  for 
9960-X.  Thus,  if  "things  were  linear"  between  Massena  and  Tibbetts  Point 
(or  at  least  as  far  south  as  Reach  #74) ,  we  would  expect  TD  records  at  Reach 
#61,  with  corrections  from  Wellesley  Island  being  applied,  to  be  essentially 
the  same  as  shown  in  figures  4-25. a.  and  4-26. a.  Resulting  cross-track 
error  plots  are  shown  for  both  sites  (using  Iroquois  Lock  data  in  both  cases 
-  the  geometry  and  courses  are  different)  in  figure  4-27. 
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Figure  4-27 


CTE  Plots  Generated  by  Applying  the  Data  of  Figures  4-25  and 
4-26  to  the  Geometry  at  Iroquois  Lock  and  Reach  *61 


4-40 


Some  disclaimers  are  in  order.  The  first  is  a  follow-up  on  the  comment 
in  Section  4.2s  "how  do  we  pick  the  reference  point?"  For  the  plots  of 
figure  4-25  and  4-26  and,  hence,  figure  4-27,  the  average  value  over  the 
30-day  period  being  considered  defined  the  reference  point.  Since  we  expect 
a  slight  variation  in  the  spatially  correlated  component  from  March/April  to 
mid-summer,  these  figures  are  optimistic:  in  practice  we  only  have  one 
"shot"  at  defining  the  waypoints.  Presumably,  since  the  shipping  season 
starts  on  1  April,  we  would  pick  reference  points  halfway  between  the 
average  values  for  March/April  and  those  for  July/August.  Since  we  do  not 
yet  have  that  data,  we  must  simply  note  the  results  indicated  by  figure  4-27 
are  optimistic.  Here,  near  the  end  of  our  comments  on  the  data,  we  see  an 
important  illustration  of  the  utility  of  having  a  model.  We  do  not  yet  have 
enough  data  to  draw  the  final  conclusion  about  performance  anywhere. 
Nevertheless,  we  have  a  framework  in  which  to  view  what  we  do  have. 
Specifically,  the  spatially  independent  component  is  cheerfully  small  and 
the  spatially  correlated  component  is  still  something  to  be  nailed  down.  We 
know  the  performance  implications  of  both  of  these  statements. 

The  second  comment  is  that  it  is  reasonable,  based  on  what  we  have  seen 
thus  far,  to  assume  what  we  see  at  Iroquois  Lock,  referenced  to  Massena,  is 
representative  of  what  we  would  see  at  Reach  #61,  referenced  to  Wellesley 
Island,  for  the  9960-W  baseline.  For  the  9960-X  baseline,  however,  what  we 
see  at  Iroquois  Lock  is  assumed  to  be  optimistic.  Although,  as  in  Section 
4.3,  it  is  still  premature  to  try  to  update  the  predictions,  we  have  gone 
about  as  far  as  we  can  go  in  this  report.  Thus,  an  update  in  one  specific 
case  is  in  order. 

In  figure  4-28,  we  provide  a  plot  of  the  9960-X  baseline  TD  at  Tibbetts 
Point,  referenced  to  Iroquois  Lock.  We  see  a  standard  deviation  of  83 
nanoseconds  -  accrued  over  a  A  DRD  of  124  km.  If  we  "remove"  a  16  nsec 
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Figure  4-28  9960-X  Differential  TD  Record:  Tibbetts  Point  Minus 

Iroquois  Lock 
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spatially  independent  component  from  this  "total  sigma"  and  assume  a 
perfectly  linear  variation  between  Tibbetts  Point  and  Iroquois  Lock  for  the 
spatially  correlated  component,  we  would  estimate  a  26  nsec  component  at 
Reach  #61,  reference  to  Wellesley  Island  (a  A DRD  of  40  km).  Adding  the  16 
nsec  component  back  in,  we  compute  a  29  nsec  total  sigma  at  Wellesley  Island 
for  9960-X.  Combining  this  with  the  expected  37  nsec  total  sigma  for 
9960-W,  and  using  the  math  of  Appendix  C,  we  compute  a  "rho"  of  0.89.  Thes>-- 
parameters  can  be  used  to  generate  the  99.9%  error  ellipse  plots  of  figure 
4-29:  our  current  best  prediction  of  the  performance  at  Reaches  #61  and  62 
with  differential  correct  ions  from  Wellesley  Island. 
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Figure  4-29  Predicted  99.9%  Differential  Error  Ellipses  at  Reaches  »61 
and  62:  Corrections  From  Wellesley  Island 


Allowing  17  meters  for  the  vessel  "half-width,"  we  are  left  with  21 
(Reach  #61)  and  19  (Reach  #62)  meters  for  guidance  error.  Although  these 
results  look  encouraging,  we  must  recall  the  "reference  point"  problem.  We 
have  just  provided  predictions  which  used  the  mean  of  the  March/April  data 
as  the  reference  point.  If  we  concede  our  99.9%  error  ellipses  are 
representative  of  "as  bad  as  it  gets,"  we  can  obtain  a  better  estimate  of 
the  year-round  situation  at  these  reaches  (we  cannot  do  a  direct  computation 
because  we  have  no  means  of  computing  "March  to  July"  statistics) . 


First  we  examine  figures  4-12  and  4-13.  We  see  the  9960-W  baseline  td 
looks  like  it  will  "steady  out"  at  about  15691.860  usee  at  Tibbetts  Point, 
with  Massena  corrections  applied.  This  is  0.258  usee  higher  than  the 
reference  value  of  figure  4-30. a.  -  Tibbetts  Point  minus  Massena  for  the 
30-day  period.  A  similar  comparison  for  the  9960-X  baseline  shows  the 
reference  in  the  summer  should  be  about  0.74  usee  above  that  of  figure 
4-30. b.  Since  the  shipping  season  starts  in  April,  the  reference  values  in 
figure  4-30  should  be  "as  lew  as  we  expect"  (actually,  somewhat  lower). 

Ideal  selection  of  the  reference  point  (for  now)  can  be  argued  to  be  midwav 
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Figure  4-30  9960-W/X  Differential  TD  Record:  Tibbetts  Point  Minus 

Massena 


between  these  extremes:  leaving  us  with  a  "reference  point  drift*  of  0.129 
usee  for  9960-W  and  0.037  usee  for  9960-X. 

To  be  conservative,  however,  let  us  stay  with  the  values  of  0.258  and 
0.074  usee.  This  is  particularly  applicable  since  we  do  not  yet  have  a  way 
to  pick  the  ideal  reference  point.  At  worst,  we  would  pick  the  largest 
value  observed  in  the  summer. 

Of  course,  we  do  not  expect  to  use  differential  corrections  from 
Massena  all  the  way  to  Tibbetts  Point.  We  are  considering  data  from  these 
sites  because  these  are  the  only  ones  with  data  records  which  span  the 
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timeframe  required  for  the  desired  computation.  Data  available  thus  far 
supports  the  appropriateness  of  a  linear  description  of  the  spatially 
correlated  component  of  the  variation  between  Massena  and  Tibbetts  Point  for 
the  9960-W  baseline.  Since  the  reference  point  drift  is,  in  full,  a 
spatially  correlated  phenomenon,  we  can  use  this  linear  model  to  estimate 
the  drift  at  Reach  #61,  with  corrections  from  Wellesley  Island,  for  9960-W. 

In  a  DRD*  sense,  Reach  #61  is  81  km  from  Wellesley  Island  whereas 
Tibbetts  point  is  321  km  from  Massena.  Thus,  we  expect  the  reference  point 
will  drift  81(258  nsec)/321  =  65  nsec  from  late  March  to  mid-sumner  at  Reach 
#61,  with  corrections  from  Wellesley  Island. 

For  9960-X,  we  have  concluded  we  do  not  have  linear  behavior  from 
Massena  to  Tibbetts  Point.  Let  us  again  assume  that  Iroquois  Lock  "stays 
with"  Massena  throughout  the  year.  If  so,  we  would  argue  the  74  nsec  9960-X 
drift  should  be  apportioned,  in  a  1  i n  r,  DR0X  sense,  over  the  distance 
from  Iroquois  Lock  to  Tibbetts  Point  (124  km  -  as  before).  Reach  #61  is  40 
km  from  Wellesley  Island  in  the  .\DRDX  sense.  Thus,  we  expect  the  9960-X 
reference  point  will  move  40(74  nsec)/124  =  24  nsec  from  late  March  to 
mid-sunnier  at  Reach  #61,  with  corrections  from  Wellesley  Island. 

Using  the  math  of  Appendix  A  and  appropriate  courses,  the  65  nsec 
aTD*  and  the  24  nsec  aTDx  values  are  transformed  to  a  ACTE  of  8  meters 
at  Reach  #  61  and  of  6  meters  at  Reach  #  62.  These  figures  "cut  into"  the 
21  and  19  meters  we  have  left  for  guidance  error  as  shown  in  figure  4-29. 

We  conclude,  however,  we  are  still  have  room  for  the  somewhat  arbitrary 
10-meters  limit  we  selected  as  part  of  the  total  error  budget. 

These  results,  of  course,  are  still  projections  from  a  suboptimal  data 
base:  we  do  not  have  have  data  from  the  desired  number  of  sites  and  we  do 
not  have  data  from  a  long  enough  time  period.  They  may  prove  to  be 
optimistic  but  they  stand  as  our  best  guess  at  present.  We  should  note 
there  are  two  ways  to  improve  upon  the  situation  should  the  need  be 
indicated.  First,  we  can  locate  a  monitor  halfway  between  Reaches  #61  and 
#74.  Next,  we  can  be  a  bit  more  careful  in  the  selection  of  our  reference 
point  -  i.e.,  try  for  an  approximation  to  the  "ideal"  suggested  above  (i.e., 
halfway  between  extremes).  Either  of  these  steps  will  cut  the  reference 
point  drift  in  half  (together  the  drift  is  reduced  by  a  factor  of  4).  As 
the  figures  of  the  previous  paragraph  show,  however,  the  total  reduction 
through  this  mechanism  may  only  be  on  the  order  of  6  meters.  A  larger 
contribution  is  made  by  locating  the  monitor  halfway  between  the  two 
extremes:  it  would  also  reduce  the  size  of  the  ellipses  in  figure  4-29. 

For  now,  we  will  simply  state  that  the  results  look  encouraging  and  note  we 
have  "untapped"  methods  of  improvement,  should  future  data  prove  we  have 
been  too  optimistic. 

The  final  statement  we  should  make  is  in  regard  to  the  "anomalous"  data 
point  indicated  in  figure  4-20.  Recall  the  data  came  from  the  period  just 
before  1  April  1982.  It  was  claimed  this  was  primarily  a  spatially 
correlated  effect  which  would  have  been  removed  by  a  closer  monitor 
station.  Examination  of  the  "Iroquois  minus  Massena"  data  for  the  same 
period  fully  supports  this  contention. 
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4.5  Summary  of  Data. 


Observations  made  at  Massena  confirm  the  year-round  predictions  made  by 
using  data  obtained  at  N.E.U.S.  coast  sites.  Importantly,  the  elliptical 
pattern  of  the  fixes  -  oriented  with  the  larger  axis  perpendicular  to  the 
course  -  is  confirmed.  By  time-domain  filtering,  we  are  able  to  remove  part 
of  the  spatially  correlated  component  to  examine  the  remnant.  The  remnant 
begins  to  approach  the  assumed  values  used  in  Section  3  for  the  spatially 
independent  component  but  too  much  of  the  remaining  "higher  frequency" 
components  appear  to  be  spatially  correlated  for  us  to  confirm  our 
assumptions.  For  this,  we  must  have  data  from  another  "nearby"  site. 

In  comparing  the  Tibbetts  Point  data  with  the  Massena  data  we  can 
remove  more  of  the  spatially  correlated  components.  More  importantly,  we 
can  remove  those  components  which  are  identical  from  one  Internav  404 
receiver  site  in  the  St.  Lawrence  Seaway  to  another.  Time-domain  filtering 
is  again  applied  to  the  "Tibbetts  Point  minus  Massena"  data  to  remove  the 
percentage  of  the  spatially  correlated  components  that  simply  could  not  be 
reioved  because  of  the  large  distance  from  Massena  to  Tibbetts  Point.  The 
result  is  promising:  the  remnant,  known  to  still  have  a  spatially 
correlated  portion,  is  smaller  than  that  assumed  in  the  predictions  of 
Section  3  for  the  9960-X  baseline  and  almost  as  small  for  the  9960-W 
baseline.  Ttie  tie-in  between  the  time-domain  filtering  and  the  spatial 
filtering  afforded  by  differential  action  is  not  yet  fully  justified  but  we 
begin  to  see  its  applicability. 

The  criticality  of  the  promising  observations  regarding  the  spatially 
independent  component  is  seen  by  recalling  Section  3  showed  the  assumed 
values,  given  the  adverse  geometry,  yielded  performance  "too  close  for 
comfort"  in  some  reaches.  One  abnormal  data  point  is  noted  but  it  is  argued 
this  is  a  spatially  correlated  effect. 

The  Tibbetts  Point  data,  particularly  for  the  9960-W  baseline,  also 
shows  variations  in  the  spatially  correlated  components  which  are  much 
larger  than  assumed  in  Section  3.  This  is  not  a  critical  problem  since  it 
can  be  overcome  by  appropriate  selection  of  the  differential  control  station 
site.  The  contrast  in  the  disagreement  with  the  DRD  model  at  Tibbetts  Point 
with  the  agreement  at  Massena,  however,  suggests  a  "local  non-linearity." 
Thus,  although  a  solution  to  the  spatially  correlated  component  variation 
problem  exists,  it  cannot  be  identified  from  consideration  of  just  the  data 
from  Tibbetts  Point  and  Massena. 

Additional  data  is  available  for  a  critical  time  period  at  Iroquois 
Lock.  By  analyzing  the  data  as  before,  we  finally  get  agreement  between 
baselines  as  to  the  proper  estimate  for  the  spatially  independent  component 
-  about  an  11  nsec  sigma.  An  important  observation  is  that  the  previously 
noted  "anomalous"  data  point  is  confirmed  to  be  spatially  correlated  (thus, 
"reducible").  We  see  evidence  that  the  bulk  of  the  non-linearity  in  the 
spatially  correlated  component,  for  the  9960-W  baseline,  has  occurred 
"before"  Massena.  Additional  data  from  the  vicinity  of  Wellesley  Island 
would  be  useful  in  confirming  this.  For  the  9960-X  baseline,  however,  we 
see  evidence  of  local  non-linearity  in  the  spatially  correlated  component. 
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This  is  a  critical  finding  and  makes  data  from  the  vicinity  of  Wellesley 
Island  crucial. 

Finally,  we  update  our  predictions  at  the  most  critical  reaches  usinq 
all  available  data.  The  results  look  promising  -  given  monitors  at 
Wellesley  Island  and  Massena.  We  note  however,  the  final  conclusions  were 
based,  in  part,  on  only  30-days  of  data  and  from  not  enough  sites  as  would 
be  realistically  required  to  draw  final  conclusions. 


s 


Considerations  on  Implementation  and  Alternatives 


^.1.  General  Statement. 


Throughout  the  report,  it  has  been  tacitly  assumed  that  full 
Differential  Loran-C  is  a  system  that  could  possibly  satisfy  precision 
navigation  requirements  in  the  St.  Lawrence  Seaway.  The  prime  emphasis  of 
the  Steering  Committee  to  date,  as  well  as  this  report  to  this  point,  has 
been  to  examine  Loran-C  system  capabilities  and  Seaway  requirements.  These 
areas  should  continue  to  receive  emphasis  but  it  is  time  to  introduce  some 
feasibility  considerations. 

Differential  Loran-C  is  a  concept.  Beyond  that  we  can  say  it  is  a 
concept  that  has  been  discussed  many  times  over  the  years  (dinners  during 
symposia,  for  example,  are  ideally  suited  for  the  discussions).  The  fact 
remains  that  Differential  Loran-C  has  never  been  implemented  as  a  general 
navigation  system  anywhere.  Prospects  for  implementation  by  the  U.S.  Coast 
Guard  are,  presently,  non-existent.  One  reason  for  this  is  the  burden  of 
simply  establishing  the  communications  network  to  support  the  concept.  It 
is  unlikely  that  implementation  in  one  small  area  (as  opposed  to  all 
continental  U.S.  coast  waters)  would  receive  acceptance.  Thus,  for  the  U.S. 
Coast  Guard,  the  first  step  towards  government  sponsored  Differential 
Loran-C  implementation  would  involve  implementation  "wherever  necessary." 
Funding  to  support  this  substantial  increase  in  service  is  unlikely. 

This  does  not,  however,  rule  out  implementation  by  other  agencies  with 
less  direct  ties  to  the  basic  system.  Given  the  inherent  "localized  area  of 
interest"  of  the  St.  Lawrence  Seaway  Development  Corporation  and  the  St. 
Lawrence  Seaway  Authority,  the  "proliferation"  problem  is  minimized. 
Moreover,  the  existing  communications  network  and  protocol  makes  things 
easier.  Feasibility  tests  can  at  least  be  conducted  without  major 
installation  costs.  Having  thus  alluded  to  feasibility,  two  major  areas 
requiring  further  consideration  should  be  mentioned. 

To  make  the  first  point,  we  should  point  out  that  Loran-C  is  a  "new" 
system.  We  should  explain  the  use  of  the  word  "new"  when  applied  to 
Loran-C.  There  is  a  dichotomy  involved  that  is  of  concern  to  us  as  well  as, 
for  example,  those  who  would  advocate  the  termination  of  Loran-C  service  as 
soon  after  the  initial  availability  of  GPS  as  possible.  Loran-C  is  "old"  in 
that  signals  have  been  on-air  for  over  20  years.  A  benefit  of  this  long 
"evaluation"  period  is  that  many  engineering/scientific  analyses  could  be 
conducted  over  the  years.  The  "refined"  system  that  replaced  Loran-A  thus 
benefitted  from  these  years  of  study  (3adly,  in  some  instances,  not 
enough).  For  the  general  user  community,  however,  Loran-C  is  new.  It  is 
just  now  5  years  that  Loran-C  has  existed  at  all  on  the  West  Coast  and 
Loran-C  without  "coverage  holes"  or  "charting  anomalies"  is  an  even  more 
recent  phenomenon  on  the  West  Coast.  Loran-C  for  the  Great  Lakes,  Gulf 
Coast  and  Canadian  East  Coast  is  much  less  than  5-years  old.  Although  the 
old  "East  Coast  Loran-C"  chain  (9930)  was  operated  for  about  20  years,  its 
widespread  utility  for  general  navigation  over  even  a  moderate  coverage  area 
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was  arguable.  Thus,  the  opportunity  for  a  "good  look"  at  the  system  by  a 
large  group  of  practical  users  is  a  recent  occurence. 

Several  factors  are  becoming  Increasingly  evident  about  this  "new" 
system.  The  major  one  for  our  purposes  is  that  a  prime  advantage  of 
Loran-C ,  when  and/or  where  Loran-C  performance  is  good,  is  that  it  is  a 
"stand-alone"  systesi.  This  is  true  about  the  most  modest  of  user  equipment 
but  even  more  so  with  those  featuring  "waypoint  navigator"  capabilities. 
These  observations  cone  frcn  general  users  and  can  be  expected  to  be  of  even 
more  concern  to  pilots.  Differential  Loran-C  is  a  direct  contradiction  of 
the  "stand-alone"  concept . 

The  second  comment  involves  the  degree  of  difficulty  involved  in 
operating  and  maintaining  high  quality  monitor  stations.  Me  can  draw  upon 
experience  gained  in  the  stability  study  of  this  report  to  underscore  the 
concern.  One  can  envision  a  "first  cut"  at  planning  this  study  proceding 
along  the  following  lines.  Note  the  digital  storage  capacity  required  for 
each  sample  to  be  obtained.  Assume  a  sampling  rate  of  once  per  second. 

Thus  aultiply  by  the  number  of  seconds  in  a  year.  Multiply  by  the  number  of 
data  collection  sites.  Observe  that  the  resulting  storage  capacity  is  well 
within  the  capabilites  of  small  systems  these  days.  Observe  the  fact  that 
Loran-C  receivers,  of  reasonable  capabilty,  are  available  on  the  market. 

Note  that  data  collection  sites  are  available.  Thus,  conclude  it  is  a 
straightforward  matter  to  conduct  a  Loran-C  stability  study. 

If  this  report  shows  nothing  else,  it  should  show  that  the  above 
"plans"  are  significantly  removed  from  what  is  truly  required  to  execute  a 
successful  study.  The  harbor  monitors  have  provided  reliable  data  - 
although  not  enough  of  it  according  to  original  plans.  A  reason  for  the 
success  is  that  they  were  developed,  over  a  period  of  several  years  and 
during  the  course  of  a  study  which  is  similar  to  this,  for  the  exact  purpose 
of  characterizing  Loran-C  stability.  Particularly  in  the  case  of  the  Ad  Hoc 
equipment  sets,  we  see  proof  that  the  marketplace  is  just  not  ready  to 
directly  serve  our  purposes.  As  good  as  the  harbor  monitor  sets  are,  they 
are  much  less  than  what  would  be  required  of  an  operational  system.  If  a 
harbor  monitor  set  fails,  there  is  no  backup.  The  failure  may  go  undetected 
for  several  days  (e.g.,  if  a  weekend  is  involved).  Even  after  discovery  of 
the  failure,  repair  is  not  completed  for  several  days.  This  type  of 
performance  represents  a  proper  commitment  of  resources  for  the  purposes  of 
the  stability  study  but  is  not  acceptable  under  operational  conditions. 

The  above  comments  can  be  viewed  as  implying  a  negativism  regarding 
Differential  Loran-C.  The  authors,  however,  intended  a  simple  note  of 
caution  which  is  deemed  appropriate  to  put  the  preceding  sections  of  the 
report  in  perspective.  The  results  of  the  stability  study  -  those  presented 
herein  and  those  to  follow,  simply  address  the  capabilty  of  Differential 
Loran-C  to  support  St.  Lawrence  Seaway  requirements.  A  test  implementation 
will  provide  insights  into  "the  rest  of  the  story."  With  the  above  notes  of 
caution  we  encourage  the  tests  to  begin. 


5.2  User  Generated  Corrections. 


Guidance  equipment  developed  under  USCG  R&D  studies,  as  well  as  many 
commercially  available  sets,  provide  a  mechanism  for  a  user  to  enter 
corrections  for  TD  biases.  The  corrections  can  be  entered  manually  or 
arrangements  can  be  made,  by  adding  appropriate  equipment,  to  have 
corrections  being  broadcast  from  some  shore  monitor  automatically  entered. 
Guidance  equipment  with  such  features  are  capable  of  operating  in  the 
Differential  Loran-C  mode.  Ttiere  are,  however,  alternate  means  of  utilizin'! 
this  correction  capability. 

A  vessel  at  the  Eisenhower  bock  at  Massena,  for  example,  has  an 
opportunity  to  observe  the  TD’s  over  a  reasonably  long  period  of  time  at  a 
known  (tied-into  the  previously  conducted  waypoint  survey)  location.  By 
comparing  the  observed  value  to  the  expected  value,  the  operator  obtains  a 
differential  correction  every  bit  as  good  as  that  to  be  obtained  from  some 
shore-based  monitor.  Indeed,  in  applying  corrections  based  on  observations 
made  by  his  own  receiver,  the  operator  also  removes  any  TD  biases  that  may 
occur  from  receiver  to  receiver.  Thus,  these  user  generated  corrections, 
for  a  time,  are  better  than  those  obtained  from  a  shore-based  station. 

In  many  areas,  this  method  of  correcting  for  systematic  Loran-C 
variations  is  entirely  adequate.  What  is  typically  missing,  however,  is  the 
opportunity  (or  willingness)  to  "sit"  in  one  known  location  for  a  roasonahh' 
period  of  time.  With  an  appropriate  investment  of  equipment  and  allowed 
system  complexity,  however,  alternate  "self-calibration"  means  can  be 
provided.  There  are,  for  example,  many  small  area  electronic  positioning 
systems  which  can  provide  high  accuracy.  Like  Differential  Loran-C,  none  or 
these  have  ever  been  demonstrated  as  a  practical  navigation  system.  For 
some,  the  problem  is  simpl-  *  at  they  cannot  cover  a  large  area  in  a 
practical  way.  They  may,  nowever,  be  useful  for  allowing  a  periodic 
"calibration  range."  Avoiding  a  deep  discussion  of  practical  implementation 
considerations,  let  us  for  the  moment  simply  assume  periodic  calibrations 
can  be  obtained,  at  specific  locations,  thus  making  shore-based  differential 
monitor  stations  unnecessary. 

Under  this  assumption,  it  becomes  important  to  ask  "for  how  long  and 
how  far"  are  these  "user  generated"  corrections  valid.  With  the  presently 
available  data  base,  we  have  a  means  to  address  these  questions.  We  take 
the  Massena  data  records  for  9960-W  and  9960-X  and  subtract  corresponding 
records  "advanced"  1-day  (2  sample  periods).  This  gives  us  a  measure  o'  -nn. 
well  a  user  who  stays  in  the  vicinity  of  Massena  will  do  with  "24-hour  old" 
corrections  obtained  at  Eisenhower  Lock.  Resulting  time  difference  records 
are  provided  in  figure  5-1.  A  scatter  plot  of  the  resulting  f.xes  is  shown 
in  figure  5-2. 
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Figure  5-1  9960-W/X  TD  Records  at  Massena  With  24-Hour  Old  Corrections 


Figure  5-2  Fix  Scatter  Plot  for  Data  of  Figure  5-1 


We  see  that  most  of  the  "seasonal  structure"  of  the  variations  has  been 
removed.  Canparing  the  statistics  of  the  "residuals"  with  those  of  the 
original  data  record  (from  figure  4  —  1 ,  "sigma"  is  0.133  usee  for  9960-W  and 
0.248  usee  for  9960-X) ,  we  see  an  improvement  of  about  a  factor  of  3  for 
both  baselines.  Although  the  scatter  plot  indicates  we  need  more 
improvement,  we  note  that  the  TD  variations  become  small  shortly  a '  *  <■ :  »►., 
shipping  season  starts. 

We  can  generate  one  level  of  improvement  over  t*e  ah  we  t.- 
comparing  the  Massena  data  with  'Massena  minus  12  *  .,!«■  to  a  ■* 
are  provided  in  figures  5-3  and  5-4. 
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Figure  5-3  9960-W/X  TD  Records  at  Massena  With  12-Hour  Old  Corrections 
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Figure  5-4  Fix  Scatter  Plot  for  Data  of  Figure  5-3 


Here  we  see  substantial  improvement  over  the  "day  old  corrections" 
situation.  Particularly  with  regard  to  the  9960-W  baseline,  the  variations 
are  approaching  the  levels  which  we  might  consider  usable.  With 
"low-density"  data,  we  can  go  no  further  in  this  sequence.  At  present  this 
is  not  a  significant  drawback  since  we  really  need  to  consider  what  happens 
someplace  away  from  Massena  with  the  "stale"  corrections  obtained  at 
Eisenhower  Lock.  For  the  30-day  period  for  which  substantial  data  is 
available  from  Iroauois  Lock,  we  can  answer  the  question.  From  the  Iroquois 
data,  we  subtract  "Massena  minus  12-hours"  data  to  obtain  the  plots  of 
figures  5-5  and  5-6. 
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Figure  5-5  9960-W/X  TD  Records  at  Iroquois  Lock  With  12-Hour  Old 

Corrections  From  Massena 
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Figure  5-6  Fix  Scatter  Plot  for  Data  of  Figure  5-5 


Here  we  see  that  there  is  a  small  amount  of  "spatially  correlated 
effects"  that  has  occurred  over  the  period  of  12  hours.  If  we  believe  what 
was  presented  in  Section  4,  we  can  say  this  is  what  we  expected.  Notice 
that  by  comparing  the  statistics  of  figure  5-3  with  those  of  figure  5-5,  we 
see  that  this  has  resulted  In  a  substantial  Increase  in  the  variations  for 
9960-W  but  only  a  mild  increase  for  9960-X.  Again  from  Section  4,  this  is 
the  situation  we  expected  for  Massena  and  Iroquois. 
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Without  comparing  Massena  or  Iroquois  data  to  Tibbetts  Point  data,  this 
is  about  as  far  as  we  can  go  with  this  type  of  analysis.  Since  Tibbetts 
Point  is  so  far  away  and,  since  the  variations  seem  substantial  even  at 
Iroquois,  reference  to  "stale"  Massena  data,  there  is  no  reason  for  further 
analysis  at  this  point.  We  should  note,  hcwever,  that  this  is  a  useful 
correction  technique  -  it  might  even  be  all  that  would  be  required  if  the 
geometry  were  not  so  poor.  We  should  continue  to  examine  this  form  of 
correction  in  future  studies. 


5.3  Changes  to  Existing  Chain  Geometry. 


Installation  and  operating  costs  for  Loran-C  transmitting  stations  are 
substantial.  Nevertheless  we  should  examine  the  implications  of  our  best 
estimate  of  what  would  happen  if  the  geometry  could  be  improved  by  the 
addition  of  another  secondary  transmitting  station.  In  considering  the 
orientation  of  the  9960  chain  LOP's  in  the  St.  Lawrence  Seaway,  as 
illustrated  in  figure  4-19,  we  see  what  we  really  need  is  a  stable  LOP 
oriented  at  about  a  course  of  225°T.  If  "a  benefactor"  were  to  install 
and  operate  a  station  in  the  vicinity  of  North  Bay,  Ontario  (near  Lake 
Nipissing)  we  would  have  an  LOP  with  about  the  desired  orientation.  If  we 
were  able  to  control  the  baseline  by  use  of  a  receiver  located  at  Massena, 
we  would  have  the  stability  we  seek. 

Appendix  G  contains  error  ellipses  which  were  obtained  by  use  of  the 
prediction  techniques  of  Section  3  -  applied  to  the  9960-W  and  9960-N  (North 
Bay)  baselines.  TO  be  conservative,  we  stayed  with  a  spatially  independent 
component  "sigma*  estimate  of  20  nsec.  Additionally,  we  used  the  9960-W 
"pk-pk  vs  DRD"  regression  line  slope  as  the  basis  of  9960-N  spatially 
correlated  component  predictions.  In  the  Appendix,  we  begin  by  assuming  no 
differential  corrections  are  applied.  In  figure  5-7,  we  show  the  resulting 
error  ellipses  for  several  reaches  of  interest. 

As  in  Section  3,  we  scale  the  95%  statistics  up  to  obtain  the  99.9% 
predictions  which  we  will,  again,  call  "conservative."  The  resulting 
cross-track  error  predictions  for  each  reach,  along  with  the  error  margin 
and  asterisks  are  shown  in  Table  5-1. 
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Me  aee  asterisks  clustered  after  Reach  #61.  Partially,  these  are  due 
to  the  suboptiaal  courses.  Primarily,  however,  they  are  due  to  the  large 
9960-W  variations  which  combine  with  even  the  slightest  9960-N  variations  to 
produce  large  ellipses  in  the  southern  reaches.  The  problem  areas  before 
Reach  #61  are  siaply  the  result  of  suboptiaal  courses.  It  is  laportant  to 
emphasize  that  these  results,  though  they  aay  not  be  fully  satisfactory,  are 
for  "raw"  Loran-C  -  with  the  new  baseline.  The  geometry  has  improved  so 
much  that  there  aay  be  satisfactory  alternatives  to  full  Differential 
Loran-C.  Before  exploring  the  alternatives,  we  should  examine  the  results 
we  could  expect  if  we  were  to  apply  full  Differential  Loran-C  corrections 
from  Massena  -  for  the  9960-W  baseline  only.  In  all  cases,  the  results  will 
be  better  than  indicated  in  Table  5-1.  Thus,  we  need  only  consider  the 
reaches  which  are  aarked  with  asterisks. 

Resulting  error  ellipses  are  also  provided  in  Appendix  G.  Conservative 
(99.9%)  cross-track  error  predictions  and  error  margins  for  the  problem 
reaches  are  listed  in  Table  5-2. 
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99.9%  Probability  Cross-Track  Error  Predictions:  9960-W/n 
With  Differential  Corrections  from  Massena  Applied  to  9960-W 


Here  we  see  we  have  "solved  all  problems"  (if  we  believe  in 
predictions):  the  largest  cross-track  error  is  27  meters.  Examining  the 
error  margins,  which  achieve  a  minimum  value  of  40  meters  at  Reach  #62, 
shows  that  for  17  meter  "half-width”  vessels,  23  meters  is  allowed  for 
guidance  error.  This  significantly  exceeds  our  assumed  requirements.  He 
have  achieved  this  performance,  however,  at  the  expense  of  requiring  a 
differential  monitor  station.  This  is  a  disadvantage.  An  advantage  over 
the  situation  provided  by  the  existing  chain  is  that  we  only  need  one 
monitor  station  and  we  have  considerable  room  for  error.  Careful 
consideration  of  the  error  margin  suggests  that  with  the  new  station,  there 
may  be  some  intermediate  step  between  the  "raw  9960-W"  and  "full 
differential  for  9960-W"  situations  which  would  meet  our  needs. 

Specifically,  we  now  recall  the  "factor  of  3"  improvement  over  the  "raw” 
data  we  were  able  to  obtain  in  the  vicinity  of  Massena  with  "24-hour  old” 
corrections  from  Massena  and  the  even  greater  improvement  if  the  corrections 
were  only  12-hours  old. 

In  Section  5.2,  we  found  that  12-hour  old  corrections  from  Massena 
reduced  the  standard  deviation  of  9960-W  to  66  nsec  -  over  a  "worst  case"  30 
day  period  at  Iroquois  Lock.  Such  a  small  variation,  combined  with  the 
expected  9960-N  variation  would  provide  acceptable  performance  throughout 
the  Seaway.  At  this  point  in  the  experiment,  unfortunately,  we  do  not  have 
a  database  to  support  further  pursuit  of  this  concept.  It  is  introduced, 
however,  as  an  important  thought  for  future  examination. 

One  final  comment  should  be  made  to  encourage  this  line  of  investi¬ 
gation  into  the  implications  of  the  improved  performance  afforded  by  the  new 
station.  One  advantage  has  already  been  suggested:  we  would  not  be 
critically  dependent  upon  real-time  corrections.  Perhaps  a  more  important 
advantage  is  that  we  would  not  be  critically  dependent  upon  data  from  any 
one  shore-based  system.  In  this  regard,  recall  we  are  leaning  towards  a 
conclusion,  for  the  existing  chain,  that  a  monitor  at  Wellesley  Island  may 
not  be  good  enough  to  allow  safe  passage  of  Reach  #61.  Presumably,  a 
monitor  right  at  Reach  #61  would  be  satisfactory.  If,  however,  there  is  an 
equipment  failure  at  that  site,  the  whole  system  fails. 

If  we  had  more  room  for  error,  as  with  the  new  baseline,  we  could 
locate  several  monitors  along  the  Seaway  and  deduce  an  "interpolated"  (in, 
of  course,  a  DRD  sense)  correction  to  be  used  in  any  reach  of  interest. 

(For  future  reference,  such  corrections  provide  what  is  refered  to  as 
"proportional  control.")  The  result  would  not  be  as  good  as  if  we  had  a 
work ing  monitor  right  at  the  reach  of  interest.  With  the  geometry  provided 
by  the  new  station,  however,  it  does  not  have  to  be.  With  the  new  baseline 
and  this  technique,  it  is  hypothesized  that  we  could  "lose"  all  but  two, 
perhaps  even  one,  of  the  monitor  sites  and  still  provide  acceptable 
corrections  at  any  reach  along  the  Seaway.  This  system  "reliability" 
feature  is  the  significant  advantage  brought  about  by  the  new  station.  With 
existing  chain  geometry,  we  are  simply  too  sensitive  to  any  increase  in 
errors  and  thus  must  demand  higher  (probably  unachievable)  reliability. 

Future  efforts  should  strive  to  obtain  adequate  data  bases  at  Tibbetts 
Point,  Wellesley  Island,  Iroquois  Lock  and  Massena  to  further  examine  the 
perforrance  of  "proportional  control  corrections"  for  9960-W. 


5-11 


6 


CONCLUSIONS  AND  RECOMMENDATIONS 


We  have  commented  extensively  on  data  set  reliability,  the  adequacy  of 
the  sampling  strategy,  the  nature  of  Loran-C  variations,  expected  system 
performance,  and  actual  observations.  In  summary,  the  following  conclusions 
can  be  drawn. 

— —  The  reliability  of  the  harbor  monitor  sets  and  the  quality  of  the 
data  base  exceeds  that  obtained  in  earlier  studies  in  the  St.  Marys  River. 

-  Except  for  a  valuable  30-day  period  at  Iroquois  Lock,  no  useful 

data  was  obtained  from  the  high  density  data  collection  sites.  Part  of  the 
problem  could  be  remedied  if  the  data  could  be  reviewed  in  a  more  timely 
fashion.  Ttiis  is  not  practical  without  some  sort  of  phone  line  access. 
Another  source  of  the  problem  is  the  enormous  size  of  the  data  base  -  errors 
can  be  generated  almost  as  fast  as  they  can  be  detected  by  the  experimenters. 

-  Without  a  reasonable  amount  of  high  density  data,  we  cannot  make 

an  in-depth  comparison  of  the  results  of  the  high  density  and  low  density 
data  collection  approaches.  Careful  review  of  the  model  presented  herein, 
however,  suggests  the  approaches  should  yield  essentially  identical 
results.  The  strong  agreement  between  the  observations  and  the  predictions 
which  were  made  using  the  model  stands  as  the  only  evidence  to  support  this 
suggestion . 

-  The  latest  version  of  the  harbor  monitor  sets  (type  D)  could  be 

used  to  collect  "higher  density"  data  than  is  used  herein. 

-  Predictions  based  on  the  model  described  herein,  along  with  data 

from  Northeast  U.S.  coast  sites,  indicate  "raw"  Loran-C  from  the  existing 
9960  chain  will  not  yield  satisfactory  year-round  performance  in  the 
Seaway.  Data  from  Massena  and  Tibbetts  Point  confirms  this. 

-  Predictions  further  show  that  Differential  Loran-C,  even  with 

several  shore  sites,  will  be  hard  pressed  to  meet  performance  requirements. 
Near  "worst  case"  geometry  is  the  reason. 

-  Data  to  directly  confirm  the  Differential  Loran-C  predictions  is 

not  available.  Using  an  elaborate  series  of  analysis  techniques,  however, 
we  are  able  to  conclude  the  spatially  independent  component  of  the 
variations  is  smaller  than  predicted.  The  spatially  correlated  component  is 
slightly  larger  than  predicted. 

-  Available  data  indicates  the  spatially  correlated  component  of  the 

9960-W  baseline  variations  increases  as  an  approximately  linear  function  of 
double  range  difference  between  Massena  and  Tibbetts  Point. 

-  Data  indicates  the  spatially  correlated  component  of  the  9960-X 

baseline  variations  increases  in  a  non-linear  fashion  between  Massena  and 
Tibbetts  Point.  With  the  available  data,  we  cannot  adequately  character iz» 
this  non-linearity.  This  baseline  is  the  prime  contributor  to  cross-track 
error  -  the  critical  performance  parameter. 
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—  Because  of  the  low  spatially  Independent  component  variations,  the 
data  Indicates  a  flawlessly  operated  and  surveyed  Differential  Loran-C 
system  could  provide  satisfactory  performance.  Given  the  uncertainty 
regarding  the  9960-X  spatially  correlated  component  and  the  scant  room  for 
error,  we  cannot  yet  say  where  the  shore  stations  should  be  located. 

Another  data  site  in  the  vicinity  of  Wellesley  Island  Is  needed. 

—  Full  Differential  Loran-C  has  never  been  implemented.  Since  the 
geometry  leaves  scant  room  for  error,  this  would  be  a  rugged  area  for  a 
first  operational  deployment.  A  determination  of  practical  limitations 
Incurred  In  actual  Implementation  is  needed  to  state  the  "bottom  Line." 

—  If  a  station  In  the  general  vicinity  of  North  Bay,  Ontario  were 
added  to  the  9960  chain,  the  geometry  would  be  improved  considerably.  “Raw" 
Loran-C  still  would  not  yield  satisfactory  performance  in  some  reaches. 

—  The  new  station,  however,  would  leave  enough  room  for  error  so 
that  alternatives  to  full  Differential  Loran-C  look  promising.  The  periodic 
application  of  proportional  control  corrections  is  an  example  of  such 
alternatives. 


Based  on  these  conclusions,  we  offer  the  following  recommendations. 


—  Type  D  harbor  monitor  sets  should  be  deployed  at  Brossard, 
Beauharnois,  Massena,  Iroquois  Lock,  Wellesley  Island,  and  Tibbetts  Point. 

—  m  increases  the  density  of  the  data  collection  routine 
should  be  tried  so  that  continuing  fears  regarding  the  inadequacy  of  the  low 
density  data  base  can  be  allayed.  If  the  percentage  of  data  available  for 
low  density  analysis  decreases  as  a  result  of  this  attempt  to  "strain  for 
gnats,"  however,  the  density  should  be  decreased. 

—  Future  data  analysis  should  concentrate  on  refining  the  model 
presented  herein  and  "hardening"  the  performance  estimates. 

—  A  test  implementation  of  Differential  Loran-C  in  the  Seaway  should 
be  undertaken. 
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Purpose. 

There  is  a  large  body  of  literature  on  geometry  considerations 
pertaining  to  lor an.  Much  of  it  is  not  applicable  to  our  purposes  for  a 
variety  of  reasons .  These  include: 

a.  Developed  for  "standard  lor an"  (Lor an -A)  not  Loran-C 

b.  Developed  for  "synchronized"  chain  operations 

c.  Developed  for  chain  planning  purposes  (addresses  "global"  vie 
specific  questions,  uses  no  "a  posteriori"  information,  etc) 

d.  Developed  to  compress  two-dimensional  information  (e.g.,  along- 
track/cr oss -track  error)  into  one  number  (e.g.,  drms,  CEP,  etc) 

etc. 

The  purpose  of  this  appendix  is  to  present  the  mathematics  used  for  the 
evaluation  of  critical  lor  an  errors  in  situations  such  as  the  St  Lawrence 
Seaway  study.  The  concepts  are  not  new  -  simply  arranged  in  the  order 
suited  to  this  report. 


We  begin  by  describing  the  "small  area"  transformation"  for  Loran-C 
tune  difference  errors  to  position  errors.  From  this  we  can  show  how  to 
present  the  positional  variations  implied  by  time  difference  fluctuations. 
We  will  consider  how  to  present  the  results  in  terms  of  "north/south"  (x,y) 
or  along- track/cross -track  error  (OTE/CTE)  coordinates. 

Next,  we  hypothesize  a  statistical  description  of  the  time  difference 
variations  and  derive  a  resultant  statistical  description  of  the  positional 
variations.  Again,  the  results  will  be  in  terms  of  (x,y)  or  (ATC/CTE) 
coordinates.  What  results  is  sensitive  only  to  the  assumptions  about  the 
time  difference  statistics  as  discussed  in  the  main  fcxxty  of  the  report. 


Snail  Area  Transformation 

We  recognize  that  the  world  is  "round  and  lumpy"  and  that  loran  LOP’s 
are  "warped  hyperbolas."  O/er  a  small  area,  however,  it  is  appropriate  to 
consider  a  so-called  "flat  earth,  linear  grid"  model  (a  related  discussion 
may  be  fotnd  in  reference  1  ).  In  such  a  case  we  can  represent  positional 
variations  in  Cartesian  coordinate  terms  -  ’  x  and  Ay.  Because  of  the 
linear  grid  assunption,  we  need  only  consider  the  slopes  of  the  loran  lines 
at  a  point.  Small  loran  time  difference  changes  are  thus  described  by  the 
following  equations: 
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In  matrix  form,  this  is 
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where  r^  «  range  to  secondary  station  1  transmitting  antenna 

rn  «  range  to  master  station  transmitting  antenna 
v  *  speed  of  signal 
Y-^  *  baseline  emission  delay. 


Both  the  emission  delay  and  v  are  time  varying  -  indeed,  such  variations 
cause  the  need  for  stability  studies  and  form  the  basis  of  the  double  range 
difference  model  as  discussed  in  the  report.  Emission  delay,  however,  is 
not  related  to  position  of  the  observer.  Hence 
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The  speed  of  the  signal  1,8  a  function  of  position  (indeed,  this  is  what 
makes  a  survey  necessary).  Over  the  small  areas  considered  here,  however, 
(e.g.,  about  50  meters)  we  can  say 
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To  consider  the  relationship  of  an  observer  (x,y)  to  a  station  of  a  Lor  an -C 
baseline  we  draw  the  following  diagram.  He  show  the  observer  at  point  (x,y) 
and  the  transmitting  antenna  at  location  (x^y^). 
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(A- 3)  gives  us  the  ability  to  express  TD  variations 
coordinates  (  Ae  *  Ax,  AN  ■  Ay)  -  assuming  we  know 
off  angles  from  the  observer  to  the  stations  of  the 
algorithm  to  compute  these  angles  will  be  presented 
with  that  description,  however,  we  should  introduce 
feature.  Besides  considering  the  mapping  from  TD's 
consider  the  related  coordinate  system: 


in  terms  of  cartesian 
the  great  circle  take- 
Loran-C  chain.  A  sinple 
later .  Before  proceding 
another  general 
to  x-y,  we  will  want  to 


ATE 

1 

*  w  * 

CTO 

The  relationship  of  concern  is  w  *  L  v  vhere 


The  -  matrix  Is  specified  by  the  vessel  course  and,  as  shown  In  the 

diagram.  Is  a  standard  "axis-rotation"  transformation. 

/  % 


w2  (CTE) 


v2  (y) 


w1  (ate) 


\  Yw 


vx  (x) 


s1nYt 


-slny. 


cos  y 


cosyw  -slnY 
s1nYw  COSY 


j  -i  nr 

J  k 


Course  Is  traditionally  specified  In  terms  of  9W  *  90  -  Yw  so  we  have 


Sln9w  cos9w 
-cos  9W  sin  9w 


(A-4) 


Thus,  given  knowledge  of  the  transmitting  station  locations  and  the  location 
of  the  observer,  we  simply  need  to  compute  great  circle  bearings,  as 
described  In  the  next  section,  to  transform  TD  variations  to  position 
variations.  With  B  described  by  (A-3)  and  L  described  by  (A-4),  we  have: 


WTiTM 

Li 


Great  Circle  Bearing  Calculations. 


The  method  for  computing  great  circle  bearings  is  fairly  common  and  a 
derivation  need  not  be  presented  herein  (see,  for  example,  reference  (,  )  . 
Since  it  is  a  critical  ingredient  for  computing  the  Loran-C  geometry  matrix, 
however,  a  simple  algorithm  for  making  the  calculation  will  be  presented. 

The  derivation  assumes  the  earth  is  a  perfect  sphere  -  an  entirely  adequate 
assumption  for  our  purposes. 


Due  to  the  form  of  the  spherical  triangles  considered  in  the 
derivations,  along  with  conventions  for  specifying  angles,  there  are  four 
cases  to  be  considered: 


Case  I 
Case  II 
Case  III 
Case  IV 


observer  is  south  and  west  of  the  transmitter 
observer  is  north  and  west  of  the  transmitter 
observer  is  south  and  east  of  the  transmitter 
observer  is  north  and  east  of  the  transmitter 


Define  the  following  quantities: 


Lj  =  latitude  of  the  observer  (Lj  is  positive  for  north 
latitudes) 

Xl  *  longitude  of  the  observer  (Xi  is  positive  for  east 
longitudes) 

1*2  x  latitude  of  the  transmitter  antenna  (same  sign  convention 
as  for  Lj) 

X  2  =  longitude  of  the  transmitter  antenna  (same  sign  convention 

as  for  Xi) 

3  =  great  circle  bearing  from  the  observer  to  the  transmitter 

antenna  (measured  clockwise  from  north) 

Compute: 

E  *  cosl>2  sin(X2  -  X i ) 

F  *  sinL2  cos  Iq  -  sinL]  cosL2  cos(\2  -  X 1 ) 

H  *  tan-1  (E/F) 


In  case  I,  3  *  H; 
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In  case  II,  3 


H  +  180°; 


In  case  III,  P  -  H  ♦  360°; 

In  case  IV,  P  -  H  ♦  180° 

A  popular  algorithm  recognizes  that  F  Is  positive  In  case  I  and  III  and  that 
H  Is  positive  In  cases  I  and  IV. 

Thus,  the  algorithm  features  the  following  steps: 


1.  Compute  E  and  F 

2.  Compute  H  using  a  single-argument  tan-1  function  (l.e., 

-900  <  h  <  900) 

3.  Let  P  -  H 

4.  If  F  <  0  (l.e..  If  case  II  or  IV  ), 

replace  p  with  p  ♦  1800 

5.  If  P  <  0  (originally  possible  only  In  case  II  or  III.  After 

step  4,  possible  only  In  case  III) 
replace  P  with  £  +  360° 

6.  Done 


Statistics. 

With  the  mathematics  presented  above,  we  have  the  ability  to  transform 
observations  of  two  time  differences  Into  useful  2-dlmenslonal  position 
references .  These  transformations  see  heavy  use  In  the  body  of  the  report. 
We  will  also  want  to  arrive  at  a  statistical  description  of  the  positional 
variations  as  a  function  of  time  difference  statistics.  If  we  can  then 
determine  a  reasonable  model  for  the  time  difference  statistical  parameters, 
and  how  th$y  vary  from  site  to  site,  we  can  extend  *he  results  to  locations 
between  monitor  sites.  The  first  step  Is  to  develop  a  relationship  between 
time  difference  statistics  and  position  statistics.  We  will  use  the 
notation: 


•  1  we 

"uf 

'ATDf 

U 

«2_ 

atd2 

As  before. 


— 

VI 

f  A  W~\ 

AX 

V  - 

■ 

v2 

W.  — 

-Ay. 

A- 7 


The  assumption  we  will  make  about  the  statistics  of  the  TD  variations  is 
that  they  ares 

jointly  distributed,  zero  mean,  correlated  (in  general)  gaussian 
randan  variables  with  a  probability  density  function: 


r 


Puru2<ul'U2> 


- o~r 

2TTaia2^1_P^ 


exp 


2(1-, .2) 


(  hl  .  ,  a 

'  O  rt  rX 


°la2 


or,  in  vector  form, 


pu(U)  *  -Jl— r  exp  (-i  ut  k"1  u) 

2tt  |K| 1 


where  K  *  the  covariance  matrix  = 


P  0^2 


p  oxo2 


2  2 

E  (U^  1  *  On 


Elu^uj]  *  P  0^2 


,  2 
E[u2  ] 


=  On 


We  have  the  probability  density  function  for  u  and  want  to  find  the 
probability  density  function  for:  v  =  g ( ij)  .  The  procedure  is  a  standard 
one  found  in  many  textbooks  (e.g.,  reference  7  )•  Applying  the  procedure 
with 

v  =  Bu  yields 


so  that 

- i - r  exp  i  -A  (B-1  V)T  (K-1)  (B-1  V)  ' 

2n  |B|  iK |  ■’  ' 
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Pv  (V) 


1 


2rr|Bj|K| 


1  exp 


-7  vVV  k'b'v 


or  pv(V) 


2tt|H|  ■ 


-T-  exP  (-2  XT  H'1  V) 


where  H'1  -  (B-l)TlH  B"1 


Interpretation 


h 

s 


t 


[ 


Having  determined  the  probability  density  function  for  the  random  vector 
(subject  to  the  validity  of  the  assumed  statistical  properties  of  u)  we  ne 
to  find  how  best  to  use  this  Information.  Typical  uses  (e.g.,  as  per  the 
main  body  of  this  report)  require  us  to: 

1.  find  arbitrary  equal  probability  contours  py(Y)  ■  k  (so, 

for  example,  we  can  plot  the  contour)  ~ 

2.  determine  the  probability  of  falling  within  that  contour  (so, 
for  example,  wo  know  what  to  call  what  we've  plotted  -e.g.,  50*  ?n>bab111ty 
contour,  95*  probability  contour,  etc) 


First  note  that  for 


we  have 


or 


Pv(V)  *  k  «  — - — r-  exp  (-|  VT  H-1  V) 

-  2tt  |Hj J 


2tt  kj  H| 2  =  exp[-(1/2)  VT  ]Hv) 

yT  h-1  V  -  -21n  (2tt  k|H| ' ) 


m|< 


yT  h~1  V  ■  c2  Is  the  equation  of  an  ellipse.  Thus,  we  know  what  to 
plolf.  Tn  general,  there  will  be  a  Vj-Vg  "cross-product"  term  (due  to 
geometry,  a  TD  correlation,  or  both)  resulting  In  a  rotation  of  the  axes  of 
the  ellipse  about  the  V1-V2  axes. 

Determining  what  contour  we  have  Is  another  textbook  problem.  As  shown  In 
reference  8  for  a  2x2  case,  the  probability  that  an  error  vector  lies 
Inside  an  ellipse  whose  equation  Is 


Aj  D  Ap  *  c2  ("concentration  ellipse") 

P  *  1-  exp[-(l/2)c2]  . for  normally  distributed  random 

variables. 


We  will  be  primarily  Interested  In  P  *  0.95 

0.95  «  1  -  exp[-(l/2)c2] 
or  exp[-(l/2)c2]  •  0.05 
c2  -  -2  1 n (0.05 ) 

Other  concentration  ellipses  are  described  by: 


=  0.50 

c2  = 

-2  In  (0.50) 

-  0.90 

c2  * 

-2  In  (0.10) 

=  0.99 

c2  ■ 

-2  ln(0.01) 

=  0.999 

c2  « 

-2  ln(O.OOl) 

etc . . . 

The  main  body  of  the  report  discusses/illustrates  the  applicability  of  the 
gausslan  assunptlon  and  how  to  obtain  estimates  of  and  •  •  Given  the 


position  of  the  observer  (Lat/Long),  knowledge  of  the  chain  stations,  we 
now  have  everything  we  need  to  compute 

H-l  .  (H)T  K-l  H 

Thus  we  have  all  we  need  to  plot  the  error  ellipses. 

Uses. 

We  can  use  the  math  developed  herein  as  follows: 

Scatter/CTE/ATE/Radlal  Error  Plots 

1.  Store  the  coordinates  of  tne  transmitter  antennas. 

2.  Store  the  coordinates  of  the  monitor  sites  (or  any  other  sites  of 
Interest). 

3.  Identify  the  site  of  Interest  and  compute  the  great  circle  bearings 
to  form  the  B  -  matrix. 

4.  Access  the  time  series  TO  data  for  the  site  of  Interest. 

5.  For  each  element  of  tne  time  series  vector,  compute: 


Plotting  the  resulting  series  (E  *  x,  N*y)  yields  the  scatter  plot 

6.  For  each  element  of  the  series  resulting  In  step  5,  above,  compute 

d  a  [  AE^  ♦  AN^]2  .  Plot  the  result  as  a  function  of  time  to  obtain 
the  radial  error  plot. 

7.  Identify  the  course  of  Interest.  Compute  the  (.-matrix.  For  each 
element  of  the  series  computed  In  step  5,  above,  compute: 

"ate"  'a  e! 

«  L 

CTE  “  (AN 

From  the  resulting  vector,  a  time  series  of  CTE  or  ATE  can  be  plotted. 


a 


Statistics  Plots. 

1.  For  a  site  of  Interest,  compute  the  B-matrlx  as  above. 

2.  Given  estimates  of  01  ,  0  ,  and  »»  ,  form  the  K-matrlx. 

3.  Compute  the  H-l  matrix  where 
H-l  »  (B-l  )T  H  H 

4.  Specify  a  desired  probability  contour,  P,  and  compute: 
c2  .  -2  ln(l-P). 

5.  The  error  ellipse  Is  specified  by  the  equation 
VT  H-1  V  -  c2 

To  put  tills  In  more  familiar  form,  let 


H-l 

and  note  v 
Thus, 

c2 


(fl2  +  f21>  Xy  + 


A-12 


APPENDIX  B 

DOUBLE  RANGE  DIFFERENCE  MODEL 


The  so-called  double  range  difference  model  represents  a  classic 
approach  to  explaining  the  nature  of  Loran-C  time  difference  variations.  It 
was  the  prime  model  used  In  the  St  Marys  River  Loran-C  Mini-Chain  Stability 
Studies  reported  In  reference  u  and  several  earlier  studies.  Rather  than 
presenting  a  full  bibliography  on  the  subject,  we  simply  note  It  has  been 
discussed  In  references  2,9,  and  10.  Recent  extensions  of  model 
applications  to  areas  outside  the  St  Marys  River  have  shown  the  need  to 
modify  the  model.  The  modifications  will  be  discussed  in  this  appendix 
along  with  a  tutorial  on  the  basic  model. 

The  Basic  Mooel 

To  begin  to  examine  Loran-C  time  difference  (TD)  variations  we  should 
examine  both  the  spatial  and  temporal  nature  of  the  processes  Involved.  As 
will  soon  become  evident,  these  tw>  aspects  are  Interrelated.  They  will  be 
separated  In  the  next  few  paragraphs  only  for  ease  of  presentation. 

Consider  the  situation  Illustrated  In  figure  B-l  below. 


p 


SAM 


Figure  B-l 

Define  the  following  quantities: 

Rm-p  =  range  from  the  master  station  to  a  point  of  Interest 

Rm-sam  ~  ran9®  from  the  master  station  to  the  system  area  monitor 

Rs-p  =  ran9«  from  the  secondary  station  to  the  point  of  interest 

Rs-sam  *  range  from  the  secondary  station  to  the  system  area  monitor 


vm-p 

vm-sam 

vs-p 

vs-sam 


ED 


average  signal  speed  over 
average  signal  speed  over 
average  signal  speed  over 
average  signal  speed  over 
baseline  emission  delay 


the  path 
the  path 
the  path 
the  path 


from  m  to  p 
from  m  to  sam 
from  s  to  p 
from  s  to  sam 


From  standard  Loran-C  theory,  we  have 


TDsam  *  ED  +  R$-sam/vs-sam  -  Rm-sam/vm-sam 
TDp  ■  ED  +  Rs-p/vs-p  '  Rm-p/ys-p 


Uniform  Propagation  Assumption 


Let  us  assune  the  speed  of  propagation  of  the 
throughout  the  service  area  of  the  baseline.  This 


vs-sam  *  vm-sam 

=  vs. 

•p  =  vm-p  *  v 

TDSam  * 

ED  + 

(Rs-sam  "  Rm-sam 

TDP  * 

ED  + 

(Rs_p  -  Rm-p^v 

signals  is  uniform 
means: 

that 


Now  let  us  assume  that  the  speed  of  propagation  changes  from  time  to  time 
but  In  a  uniform  manner.  By  this  we  mean  that  If  the  speed  of  propagation 
IF~time  Ti  is  V]  over  some  portion  of  the  service  area,  it  is  vi  over 
all  of  the  service  area  at  Ti.  At  some  other  time,  T?,  we  have  some 
other  speed,  V2,  everywhere  in  the  service  area.  Under  these 
cl rc instances  we  would  have: 


At  T1f 


™l,sam  =  ED  +  (Rs-sam  "  Rm-sam)/vl 
TDi  fp  *  ED  +  (Rs_p  -  Rm.pl/V! 
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at  T2, 

T^2,sam  a  ED  +  (Rs-sam  ~  ^-sam^/v2 
^2,p  =  ED  +  (R$-p  -  Rm-p^/v2 

The  above  equations  may  seem  to  Imply  that  the  emission  delay  Is  held 
constant  and  that  the  TD  at  the  SAM  varies  from  time  to  time.  As  a 
practical  matter,  Loran-C  chains  are  not  operated  this  way.  Control  is 
exerted  to  adjust  the  emission  delay  so  that  the  TD  at  SAM  remains  constant 
(over  an  averaging  period  of  several  hours).  Rather  than  choosing  some 
notation  such  as  ED(t)  to  Imply  emission  delay  Is  a  time-varying  function, 
we  will  let  the  quantity  ED  represent  some  fixed  value  which  is  the  true 
emission  delay  at  some  arbitrary  point  in  time  and  incorporate  the  term 
ED  +  A  ED  in  the  equations  for  any  other  point  in  time.  Thus,  taking  into 
account  the  realities  of  chain  control,  but  still  under  the  assumption  of 
uniform  propagation  variations,  we  have: 


At  Tlt 

TDl  ,sam  “  ED  +  (Rs-sam  "  Rm-sam^vl 

at  Tj, 

TD2,sam  ED  +  aED  +  (RS-sam  '  Rm-sam^v2 

Since  chain  control  action  forces  TD!>sam  =  TD2>sam  -  a  constant, 
we  have 


ED  +  (Rs-sam  ~  Rm-sam^vl  =  ^D  +AED  +  (Rs-sam  "  Rm-sam)/v2 


or, 

aED  =  (Rs-sam  "  Rm-sam^  H/V]  -  l/v2)  (B-l) 


Meanwhile,  the  situation  at  the  point  of  Interest  is: 


TD] fp  =  ED  +  (Rs_p  -  Rm-p)/vl 
TD2,p  =  ED  +  aED  +  (Rs.p  -  Rm-p^/v2 
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The  TD  change  at  the  point  of  Interest  Is  A TDp 


ATOp  «  T02,p-TDi,p 

-  AED  +  (Rs.p  -  Rw-p)  (l/v2  -  1/vi) 

aEO  a  -(Rs-sam  "  Rm-sarO  ^/v2  ~  l/v-|) 


ATDp  *  C(Rs-p  ”  Rm-p)  _  (Rs-sam  "  ^m-sam^  H/V2  -  1/vi) 

The  quantity  In  brackets  Is  a  function  of  two  range  differences  and  Is  the 
so-called  double  range  dl fference  {DRDp).  Thus, 


ATDp  *  DRDp  (l/v2  -  1/v-j ) 


( B-2 ) 


where 


DRDd  =  (Rs-d  "  Rm-p)  '  (Rs-sam  "  Rm-sam^ 


We  should  note  an  approximation  which  simplifies  the  math: 


let  V]  *  v 


v2  3  v  +  Av 


Then, 


l/v2  -  1/V!  =  1 /{ v  +  Av)  -  1/v 


-  Av/(v2  +  V  Av) 


A  large  data  base  shows  that  for  Loran-C,  I ^ vmax I  ^  v  so  that. 


v2  +  v  av  ~  v2 


1/k  (  .  i  I  I  V<  I  i'll'.  I  .III!  ) 


aTDr  ■=  -k  av  DRDp 


(B-3) 


I 


We  should  note  that  if  the  point  of  interest  is  on  the  same  hyperbola  as  the 
SAM, 


Rs-p  -  Rm-p  *  Rs-sari  -  Rm-sam 


so  that 


DRDp  =  0 

and,  with  uniform  changes  in  the  speed  of  signal  propagation, 

ATOp  *  0. 

In  practice,  the  uniform  propagation  assumption  is  not  entirely  valid. 
Variations  on  the  basic  theme,  however,  can  prove  extremely  useful  for 
explaining  a  high  percentage  of  observed  Loran-C  TD  variations.  Before 
exploring  observations  of  the  actual  variations,  however,  we  should  finish 
the  theme  by  presenting  some  basic  rules. 

In  general,  when  it  gets  cold,  the  signals  travel  faster.  Thus  If  we 
establish  the  summer  as  our  reference  time,  Av  is  positive  in  the  winter. 

Under  normal  circumstances,  the  location  of  the  SAM  is  fixed  so  that 

DRDp  =  (Rs-p  -  Rm-p)  ‘  (Rs-sam  *  Rr-sam) 


Rs-p  "  Rm-p  "  M 


A  basic  loran  feature  is  that  for  a  fixed  speed  of  propagation,  the  maximum 
TD  reading  Is  achieved  at  the  master  station  or  on  the  baseline  extension 
"behind"  the  master  station.  The  same  logic  that  leads  to  this  result  can 
be  applied  to  the  above  equation  to  show  that  the  maximum  double  range 
difference  is  achieved  at  or  "behind"  the  master  station  where: 


DRDmax  *  (Rs-m  -  0)  -  k-|  =  B  -  M 

(rfiere  $  *  baseline  length  in  standard  loran  notation. 


Thus,  DRDffl  is  positive  and  represents  the  largest  (most  positive) 
value  achieved  anywhere  in  the  service  area. 

Similarly,  at  the  secondary, 

DRDS  ■=  DRD„,tn  *  -g-ki- 

This  is  the  smallest  (most  negative)  value  achieved. 
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Since: 

aTD  ■  -k  Av  ORD  and 

DRD  ts  positive  at  the  nster,  negative  at  the  secondary,  and 
Av  Is  positive  In  the  winter, 

we  expect  the  TO  at  the  master  station  to  be  smaller  In  the  winter  than  It 
Is  In  the  simmer.  Similarly,  at  the  secondary  station,  we  expect  the  TD  to 
be  larger  In  the  winter  than  It  Is  In  the  summer.  Thus,  we  would  expect 
year  long  plots  of  the  TD's  at  either  end  of  the  baseline  to  follow  the 
pattern  shown  In  figure  B-2. 
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Figure  B-2 


Of  course,  temperature  Is  not  the  only  factor  affecting  signal  speed  and 
even  If  It  were,  we  would  not  expect  to  see  such  "smooth"  variations  as 
shown  above.  Ac tual  TO  observation  at  either  end  of  a  baseline,  showing  a 
more  "ragged"  pattern,  will  be  presented  later.  Before  considering  these, 
however,  we  should  examine  some  other  features  of  the  simplified  situation 
shown  In  figure  B-2. 

The  TD  variations  presented  for  the  two  ends  of  the  baseline  In  figure 
B-2  are  "equal  and  opposite."  This  Implies  that  TDm  +  TDS  ■  a  constant 
value  throughout  the  year.  If  the  uniform  propagation  assumption  were  true. 
It  Is  claimed  this  would  imply  the  SAM  Is  located  on  the  perpendicular 
bisector  of  the  baseline  as  shown  in  figure  B-3.  Let  us  check  this  claim. 


Figure  B-3 

From  (B-l) , 


^ED  *  (RB-sam  '  I  Vvi  ~  i/vj) 

-  0  for  any  velocity  since  RB_BBIB  -  Rt,-,an 

Thus,  under  uniform  propagation  conditions,  locating  SAM  on  the  baseline 
perpendicular  bisector  results  in  a  constant  emission  delay. 


should 

note 

that 

TDb 

-  ED  + 

<Rs-m 

-  Rji-m)/v 

■  ED  +£  ED  ♦ 

* s-m /v 

and 

TDg 

-  ED  ♦ 

(Rs-S 

*  Rm-s>/v 

■  ED  +  AED  - 

R«-s/v 

so  that 

(TDfc  +  TDa)/2  -  (ED  +  A ED  +  RB-n/v  ♦  ED  ♦  AED  -  RB_m/v)/2 

-  ED  +  AED 

Thus  the  average  of  TDa  and  TDB  is  a  direct  measure  of  the  emission 
delay.  Independent  of  v,  at  any  point  in  time.  When  SAM  is  on  the  baseline 
perpendicular  bisector,  A  ED  ■  0,  so  that 

(TDb  +  TDgl/2  ■  a  constant  and  any  TD  variations  at  the 
master  and  secondary  must  be  equal  in  magnitude  and  of  opposite  sign  as 
previously  claimed. 

If  the  SAM  is  not  on  the  perpendicular  bisector  of  the  baseline,  a  situation 
such  as  shown  in  figure  B-4  occurs. 
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Figure  B-4 
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In  generating  this  example,  we  placed  the  SAM  on  a  hyperbola  which  crosses 
the  baseline  at  a  location  331  of  the  way  from  the  master  to  the  secondary 
as  shown  In  figure  B-5. 


In  figure  B-5, 


A^Dm  »-lcAvDRDm  ■  -kA  v  [(Rs«ib  -  -  (Rs-sam  "  Rm-sam^ 

«-kA  v  (Rm.s  ♦  Rffl-s/3) 

-  -(4/3)  k A  V  Rh,_s 

A  TDS  =  "  k  A  v  DRDS  «-kAv  C(RS„S  -  Rn_$)  -  (Rs-sam  "  Rm-sam^ 

«-kAV  ( -Rm- s  +■  Rm-s^^ 

*  (2/3)  kA  V  Rm.s  -  -  ATUm/2 

In  this  case, 

A  ED  *  (ATDin  ♦  ATDs)/2  -  (1/2)  (-4/3  +  2/3)  kAvRm_s  =  (-1/3)  kAv 

At  this  point  we  can  summarize  the  conclusions  drawn  from  consideration  of 
smoothed  uniform  propagation  TD  plots: 

1.  By  averaging  the  TD  observations  at  each  end  of  the  baseline,  we 
obtain  a  measure  of  the  baseline  emission  delay, 

ED+AED  «  (TDi„  +  TDs)/2 

2.  With  the  SAM  located  on  the  perpendicular  bisector  of  the  baseline, 
the  emission  delay  Is  held  constant: 

A  ED  -  0  when  Re_sw,  «  R„,_sam  and  the 
variations  at  either  end  of  the  baseline  are  equal  and  opposite, 

ATtfo  *  *  ATDs 

3.  With  the  SAM  located  off  the  perpendicular  bisector  of  the 
baseline,  the  emission  delay  varies,  and  variations  at  the  end  of  the 
baseline  closer  to  the  SAM  are  proportionately  smaller  than  the  variations 
at  the  other  end  of  the  baseline: 

for  Rs-sam  '  Rm-sam  =  h 

ATI^n  - -(1-h)  k  AV  Rm.j 
and  ATDS  *  {l+k-j)  k  Av  R^j 
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so  that 


ato*  - 

-(1  -  ki)  ATDS/(1  ♦  ki ) 

(B-4) 

and 

AED  - 

(ATD„,  ♦  aTDs)/2 

■  (-1  +  ki  +  1  +  ki)  k  AV  Rm-s/2 

*  ki  k  Rfl.j 

-  -  kT  TDm/d-ki )  -  ki  6TDs/(l+ki) 


Example: 

Suppose  SAM  Is  closer  to  the  secondary  station  and  that  the  difference  In 
distance  from  the  SAM  to  M  and  S  Is  80%  of  the  baseline  length.  Suppose 
over  the  course  of  a  year,  the  "Baseline  Length"  (In  usee)  changes  1.0 
usee.  Then  ki  *  -0.8  and  kAv  Rm_s  *  1.0  usee  so  that 

ATDm  «  -1 .8  usee 

aTDs  *  0.2  usee 

and  a  ED  *  -0.8  usee 


Actual  Observations. 

The  foregoing  discussions  were  aimed  at  presenting  the  basic 
considerations  In  examining  Loran-C  time  difference  variations.  To  be  mc-e 
realistic,  however,  we  must  realize  that  smooth  variations  such  as  were 
shown  In  figures  B-2  and  B-5  will  not  occur  since  weather  variations,  the 
prime  cause  of  TO  fluctuations,  do  not  occur  In  such  a  smooth  fashion  ovc* 
the  course  of  a  year.  Moreover,  If  we  acknowledge  that  the  signal  speed 
changes  from  time  to  time.  It  Is  not  realistic  to  claim  these  changes  occur 
Instantaneously,  In  the  same  amount,  throughout  the  service  area.  This, 
however.  Is  no  reason  to  abandon  the  model  -  we  need  simply  acknowledge  v:s 
shortcomings  and  take  care  when  making  applications.  For  paths  such  as  a*® 
Involved  In  the  St  Lawrence  Seaway,  for  example,  the  model  can  be  very 
useful . 

To  obtain  some  of  the  utility  of  the  model  we  can  simply  consider  what 
happens  during  the  coldest  month  of  the  year  as  compared  with  the  wannest 
month  of  the  year.  We  thus  consider  the  two  extremes  and  do  not  get 
ourselves  hopelessly  confused  trying  to  account,  on  a  day-to-day,  or  eve" 
week-to-week  basis,  for  what  happens  In  between.  Thus,  we  do  not  try  to 
predict  or  model  weather  In  a  deterministic  sense.  We  can  and  should 
examine  the  "In-between"  periods  -  but  only  In  a  statistical  sense. 

We  begin  the  consideration  of  "real  observations"  with  figure  B-6.  T'-' 
data  was  obtained  at  the  master  and  secondary  stations  of  the  9960-X 
baseline  during  calendar  year  1981. 
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Figure  B-6 


The  "raggedness"  of  the  data  records  Is  very  much  In  evidence  -  most 
noticeably  in  the  winter  months.  We  note,  however,  that  the  underlying 
trend  follows  that  as  shown  In  figures  B-2  and  B-5:  the  master  TD’s  are 
larger  In  the  summer  than  In  the  winter  and  the  secondary  TD's  are  smaller 
in  the  summer  than  in  the  winter.  If  we  "mark"  the  January  and  July  average 
TD's,  we  see  that  the  change  from  simmer  to  winter  (summer  minus  winter)  Is 
about  +1100  nsec  at  M  and  about  -250  nsec  at  S.  Thus, 

ATDm/ATDs  =  -110/25  =  -4.40 
From  (B-4)  we  have 

A TDrj/A TDs  =  - ( 1  -k-|  )/(l+k! ) 

We  would  thus  estimate  ki  as  follows, 


1  -  ki  =  4.40  (1  +  k-|) 

1  _  k]  =  4.4  +  4.4  ki 

ki  =  -3. 4/5. 4  =  -0.63 

so  that  we  estimate  Rs-sam  -  Rm-sam  =  -0.63  Rm-s 
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This  suggests  that  SAM  Is  very  close,  hyperbolic ally  to  the  secondary. 
Specifically,  Me  Mould  guess  SAM  Is  on  a  hyperbola  which  passes  through  the 
baseline  18.51  of  the  baseline  length  from  the  secondary  and  81.51  of  the 
baseline  length  frost  the  vaster  as  shown  In  figure  8-7.  (0.185  -  0.815  » 
-0.63). 


Figure  B-7 

As  It  turns  out,  the  actual  situation  for  the  9960-X  baseline,  with 
Master  at  Seneca,  N.Y.,  secondary  at  Nantucket,  Ma,  and  SAM  at  Sandy  Hook, 
N.J.,  Is  as  shown  In  figure  B-8. 


590  Kn 


Seneca 


Nantucket 


Figure  B-8 


In  actuality,  SAM  Is  very  nearly  on  the  perpendicular  bisector  of  the 
baseline: 


Rs-sam  *  ^w-san  ■  0*02 


and  our  model  seems  to  have  failed  us. 


Closer  scrutiny  of  the  situation,  however,  shows  the  path  from  the 
Master  to  the  Secondary  and  the  path  from  the  Master  to  the  SAM  are  almost 
exclusively  over  land.  Moreover,  the  land  paths  Involved  are  Interior  New 
York  State  and  Massachusetts  -  areas  which  normally  feature  sub-freezing 
temperatures  In  mid-winter.  The  path  from  the  secondary  to  the  SAM, 
however.  Is  almost  entirely  comprised  of  seawater  which  has  a  significant 
moderating  effect.  The  particular  effect  Is  that  dry,  sub-freezing  surface 
temperatures  do  not  occur. 

These  observations  suggest  a  modification  to  the  double  range 
difference  model  In  cases  where  paths  Involve  seawater  that  does  not  freeze 
We  have  experimented  with  several  modifications  and,  given  the  presently 
available  data  base,  feel  the  best  modification  to  make  Is  to  not  count  salt 
water  paths  at  all  In  calculating  ranges  (at  least  In  the  temperate  zones). 
If  we  apply  this  modification  to  the  9960-X  baseline,  we  find  the  "modified 
ranges"  of  concern  are: 


Seneca-to-Nan tucket  *  R^j  *  526  km 

Seneca  to  Sandy  Hook  *  Rm_saro  =  332  *cm 

Nantucket  to  Sandy  Hook  *  Rs-sam  = 


Thus  we  have 

(Rs-sam  ~  Rm-sam^m-s  *  -0.60 


which  agrees  very  closely  with  the  estimated  value  of  -0.63. 

Thus,  the  efficacy  of  this  modification  to  the  DRD-model  can  be  seen. 
Further  evidence  can  be  seen  if  we  tabulate  the  observations  at  all 
available  data  collection  sites  for  the  9960-X  baseline.  In  table  B-l  we 
show  the  sunnier- to -win  ter  variations  for  each  site  along  with  both  the 
"simple  DRD"  and  "modified  DRD. " 

The  data  Is  also  plotted  In  figures  B-9  and  B-10.  In  figure  B-9,  the 
"Peak-Peak  seasonal  variation"  Is  plotted  versus  the  simple  double  range 
difference.  In  figure  B-10,  It  is  plotted  versus  the  modified  double  rang® 
difference.  A  linear  regression  analysis  shows  application  of  the  s^mp1® 
double  range  difference  model  reduces  the  data  record  variation  by  about 
factor  of  5.5.  Application  of  the  modified  double  range  difference  mo'V 
yields  an  additional  factor  of  3.4  reduction  so  that  the  standard  dev’a" 
of  the  residuals  is  only  21  nsec. 
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These  results  Indicate  the  modified  double  range  difference  model  Is  a 
very  good  one  for  the  seasonal  (l.e.,  the  largest)  component  of  the  Loran-C 
time  dlt'erence  variations.  This  model  Is  used  In  the  main  body  of  the 
report  to  make  predictions  based  on  both  the  1981  9960-W  and  9960-X  baseline 
observations  made  at  harbor  monitor  sites  and  from  1981/82  observations  made 
at  Massena,  N.Y.  More  detailed  application  of  the  simple  model  was  made  In 
reference  u  .  At  present  we  are  developing  modifications  to  the  model  to 
deal  more  with  the  statistics  of  the  short  term  variations  than  the  actual 
data  variations  as  was  done  In  reference  4  .  The  underlying  concept  of 
this  extension  Is  that  the  "raggedness"  of  the  data  records  Increases  as  a 
function  of  the  double  range  difference  just  as  the  peak -peak  variation  over 
a  year  did. 

At  the  time  of  this  writing,  we  have  obtained  a  full  year  of  Loran-C 
chain  9960  data  In  electronic  form  from  only  three  sites.  Thus  we  feel  It 
would  be  premature  at  this  time  to  present  the  short  term  variation  model. 
Publication  of  this  further  extension  of  the  model  Is  saved  for  a  future 
report  when.  In  particular,  we  hope  to  have  a  substantial  St  Lawrence  Seaway 
region  Loran-C  data  base. 


B-15 


APPENDIX  C 


A  PRIORI  ESTIMATES  OF  ST.  LAWRENCE  SEAWAY  LORAN-C  STATISTICS 


APPENDIX  C 


A  PRIORI  ESTIMATES  OF  ST  LAWRENCE  SEAWAY  LORAN-C  STATISTICS 


Purpose. 

In  Appendix  A  we  developed  how  to  generate  Loran-C  error  ellipses  and 
showed  how  we  must  first  obtain  estimates  of  aw,  ax,  and  p  -  the 
standard  deviations  of  the  two  prime  time  differences  (TD's)  In  the  area  and 
the  correlation  coefficient  between  the  two  TD's.  Ultimately,  we  want  to 
obtain  these  estimates  by  direct  observations  at  data  collection  sites  in 
the  Seaway.  In  section  3  of  this  report,  however,  we  try  to  present  an  a 
priori  assessment  of  the  Loran-C  performance  to  be  expected  In  the  Seaway. 
This  serves  to  focus  our  attention  on  the  vital  statistics  to  obtain  from 
the  stability  study  end  Illustrate  how  these  are  to  be  used.  From  this,  tne 
overall  experiment  strategy  can  be  carefully  examined. 

Since  we  do  not  yet  have  an  adequate  data  base,  the  analysis  of  section 
3  must  draw  statistical  parameter  estimates  from  some  other  source.  The 
double  range  difference  concepts  presented  In  Appendix  B  can  be  applied  to 
N.E.U.S.  Loran-C  station  and  harbor  monitor  sites  to  provide  a  first-cut 
estimate  of  the  standard  deviations.  The  methodology  to  do  this,  along  wit;1 
a  procedure  for  obtaining  first-cut  estimates  of  the  correlation 
coefficients  Is  presented  herein. 

Formulation  -  Non-Differential  Loran-C. 

We  begin  by  noting  that  section  3  explains  that  the  statistical 
parameters  of  Interest  are  time-varying  quantities.  For  our  first  cut  at 
estimating  performance,  however,  we  will  ignore  this  technicality.  We  say 
that  at  any  point  and  at  any  time,  the  time  difference  measurement  is  a 
randan  variable: 


TDX  *  k  +  xi  +  X2 


where  k  =  the  non-varying  component  of  the  time  difference 


x-|  =  the  seasonal  component  of  the  time  difference,  a  zero-mean 

random  variable 


X2  =  the  sum  of  all  shorter  term  components,  a  zero-mean  random 
varlabl  e 


We  will  assum*  independence  between  the  varying  components  so  that  the  tota- 
standard  deviation,  c  x,  Is 
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Data  obtained  at  USCG  harbor  monitor  sites  suggests  that  for  double-range 
differences  such  as  those  encountered  In  the  St  Lawrence  region  for  the 
9960 -W  and  X  baselines,  crxl  is  at  least  3  or  4  times  as  large  as  crx 2- 

If  ax]  =  3  ax2, 

ax  =  axl  [1  +  (l/3)2]l/2  *  1.05  axl  »  <jx1 

For  a  "Raw  Loran-C"  (vice  differential  Loran-C)  analysis,  therefore,  we  can 
Ignore  the  short-term  component  In  generating  a  year-round  performance 
prediction . 


To  obtain  an  estimate  of  ax] ,  we  will  simply  assume  we  have  pure 
sinusoidal  seasonal  components.  Thus, 


<JX]  =  0.707  (peak-peak  var1at1on)/2 
(i.e.,  rms-value  of  a  sinusoid). 


By  using  the  results  of  the  double-range  difference  regression  analysi 
performed  using  N.E.U.S.  chain  and  harbor  monitor  site  data,  we  can  obtain 
estimates  of  crw  and  ax  at  the  center  of  all  the  St  Lawrence  Seaway 
reaches  of  Interest.  Appropriate  double  range  differences  and  resulting 
"seasonal  sigmas"  (a  s)  are  listed  in  Table  C-L.  In  the  computations,  we 
used: 


(pk-pk)w  =  0.993  DRDW  -  45  (In  nanoseconds) 


(pk-pk)x 


1.316  DRDv  +3 


(In  nanoseconds) 


RMCh  # 

000-0 

fr*41ct*4 

SwmmI 

!Slg Lg 

m-i 

Predicted 

SmhmI 

ftMch  « 

00-8 

Predicted 

SmmriI 

*5104-0* 

00-1 

P»9dlctM 

taiuwi’ 

2 

275  It* 

.000  IMC 

401  ka 

.107  iwc 

43 

617 

.10 

10 

.774 

3 

2M 

.004 

301 

.190 

44 

119 

.10 

06 

.27$ 

4 

206 

.000 

414 

.104 

46 

524 

.10 

6B 

.776 

S 

303 

.001 

420 

.196 

46 

B7 

.169 

03 

177 

7 

313 

.004 

423 

.198 

47 

031 

.170 

595 

.278 

0 

321 

.007 

429 

.201 

40 

034 

.17? 

07 

*79 

11 

320 

.099 

4X 

.204 

B 

Ml 

.174 

900 

.<80 

12 

337 

.103 

444 

.toe 

51 

01 

.170 

06 

.283 

14 

352 

.100 

462 

.215 

02 

157 

.10 

610 

.285 

IS 

363 

.100 

463 

.216 

S3 

02 

.101 

612 

.206 

16 

3C2 

.112 

472 

.221 

55 

068 

.103 

616 

.20 

17 

369 

.114 

478 

.224 

K 

073 

.10 

619 

.209 

16 

301 

.110 

486 

.227 

57 

687 

.IB 

628 

.293 

19 

392 

.122 

495 

.232 

50 

B4 

.10 

638 

.218 

20 

396 

.123 

499 

.233 

• 

613 

.10 

•44 

.301 

21 

401 

.125 

B1 

.234 

50 

617 

.201 

644 

.30? 

22 

407 

.127 

504 

.236 

01 

620 

.202 

647 

.30? 

23 

414 

.129 

B* 

.230 

02 

62? 

.20? 

M9 

.304 

24 

419 

.131 

613 

.240 

0 

619 

.206 

65? 

.301 

25 

425 

.133 

519 

.142 

04 

636 

.207 

03 

.305 

06 

431 

.136 

122 

.244 

55 

648 

.212 

69 

.307 

27 

434 

.137 

523 

.244 

67 

05 

.214 

01 

.309 

20 

436 

.117 

527 

.246 

0 

667 

.216 

667 

711 

29 

439 

.138 

529 

.247 

71 

01 

.223 

674 

.3' 5 

31 

443 

.140 

533 

.249 

74 

701 

.230 

683 

.319 

33 

467 

.148 

554 

.259 

77 

713 

.234 

694 

.3?4 

34 

474 

.150 

558 

.280 

70 

710 

.236 

695 

.3?* 

35 

478 

.152 

162 

.263 

79 

734 

.242 

07 

•  J?5 

36 

404 

.154 

566 

.264 

B 

75? 

.248 

706 

.330 

37 

400 

.156 

570 

.266 

01 

7M 

.25? 

71? 

.33? 

30 

492 

.157 

572 

.267 

02 

779 

.20 

718 

.335 

39 

406 

.158 

S74 

.050 

03 

7*6 

.04 

727 

.339 

40 

600 

.in 

176 

.259 

41 

506 

.162 

578 

.270 

42 

$12 

.164 

801 

.27? 

table  C-1 
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All  that  remains  now  Is  for  us  to  obtain  estimates  of  the  correlation 
coefficients.  Note  that  we  have: 


c 


L 


TD  *  k  +  w,  +  w  I'D  =  k  +  x,  +  x„ 

w  w  1  2  x  x  1  2 

E[TD*  TDX]  ■  P Gy  ^x 

Table  C-1  lists  our  estimates  f  or  aw  and  ax  so  that  pis  the  remaining 
unknown.  As  before,  we  assume  the  seasonal  and  short  term  components  are 
Independent  so  that 

E[w]  W2]  *  E[x]  X2l  *  E[w|  X23  *  E[x*|  W2]  =  0 
This  leaves 

p  crx  =  E[TI^  TDX]  “Pi  °wl  °X1  +  P2tJw2°x2  (C-l) 

where 

Pi  awl  axl  *  E£"l  xl3. 
and  p2aw2ax2  *  E[w2  x23 

Since  we  have  assumed  pure  sinusoids  for  wi  and  X2  we  might  as  well 
assume  they  are  perfectly  In  phase  so  that  p 1  ■  1.  (Actually,  harbor 
monitor  data  suggests  this  Is  a  valid  assunptlon).  The  search  now  focuses 
on  estimating  p?.  To  obtain  a  useful  a  priori  estimate  of  p  2  we  should 
note  the  following. 

Let 

signal  time -of -arrival  variation 

M  II  II 

I  II  II 

-  WTOA. 

nT0A)(xT0A  wrOA)] 

■  eCwtqa  xtoa3  -  e[wtoa  ">toa3 

2 

-  ECxjoa  *TQa3  +  ECmjoA  3 

We  assume  that  all  short  term  TOA's  measurements  are  Independednt  so  that. 
ECwfOA  xT0a3  =  E[wtqa  n»T0A3  *  E[xTqa  mT0A3  3  0 


"TOA  * 

short 

term  component  of  master 

*T0A  * 

II 

"  "  "  whisky 

xT0A  * 

II 

"  "  "  xray 

Since  W2 

=  "TOA 

-  mroA  x2  =  xT0A 

E[W2  X23  =  P  2  aw2  ax2  =  E[(wtoa  - 
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Thus, 


E[w2 

X2]  *  P; 

But, 

aw2 

*  ( °wT0A 

and 

CM 

X 

o 

=  (axT0A 

so. 

N 

L 

L 


mTOA 


[(a. 


2  2x  ,  2 

vfTOA  °mTOA  '  ^xTOA 


JmTOA 


:>] ' 


Classical  treatments  assume  that  short  term  TOA  variations  are 
signal -to-nolse  ratio  dominated.  Thus,  If  the  observer  Is  close  to  the 
master  station,  so  the  argument  goes,  %T0A  <<:  °WT0A*  °xT0A  so  That 
we  can  sayP2«  0.  In  the  middle  of  the  chain  coverage  area, 


tfmTOA  ~  tfwTOA  «  <?xT0A 


so  we  have 


a2 

[(2a2)  (2 o2)]  2 


_  x 


B 


These  treatments  thus  Ignore  the  effects  of  imperfections  of  chain 
operations.  Including  timing  correction  granularity,  phase  modulation,  etc. 
In  a  more  recent  treatment  (reference  5  ),  llgon  suggests  that  for  the  short 
baseline  chains  of  Canada  and  CONUS,  unless  one  Is  truly  In  a  "fringe  area," 
(we  are  not).  It  Is  best  to  simply  assume  all  short  term  noise  Is 
"transmitter  noise"  and  thus  Independent  of  distance  to  the  stations.  Thus, 
we  say  amT0A  «  ^wTOA  *  OxTOA  «  a  and  we  get  p 2  =1/2.  This 
assumption  is  suggested  as  the  most  appropriate  a  priori  estimate  and  the 
one  we  will  use. 

If  we  let  aw]  =  kw  aW2 
and  o v-w  =  ky  o x2 , 

we  can  re-write  equation  C-l : 


(>  1  °wl  °xl  *  tJ  2  °w2  ax2 


a  a 

W  X 


C-5 


0w  =  <°wl  *  0w2  >  ’  °x  -  <°xl  +  ax2  >  J  * 


0w2  ^  k  ^  °wl  * 

w 


0x2  ^  k  ^  °xl  , 

x  ’ 


we  have 


p  _  °»1  °xl  CP1  *  P2^kwkx^] _ 

°wl  °xl  !>  +  ^  <l  +  !Ax>] 


Pi  k«  kx  ^2 

[(k»+1>  <kx"1)] 


For(J  i  «*  1  andp  2  *  0.5, 

k  k  +  } 

We  previously  argued  It  was  reasonable  to  assume  kw  -  kx 
least).  For  these  values  we  have 

9  +  1/2 

►J  *  ■  0.95 

9  +  1 

A1  ternatlvely,  for  kw  >  kx  *  2,  we  have 


*  3  (at 


4  +  1/2 
4  +  1 


In  the  St  Lawrence,  the  r  «  0.95  assumption  yields  smaller  errors  r j 
the  p*  0.90  assunptlon.  Thus,  to  be  conservative,  we  will  use  =  0.?r:  * 
all  reaches.  With  this  final  assumption,  we  have  the  ability  to  general 
predicted  error  ellipses  for  each  reach  of  the  seaway  and  estimate  maximur 
cross-track  errors  -  at  some  specified  probability  contour.  For  a  f<-sr 
cut,  we  use  a  probability  of  0.95  -  further  discussion  on  this  subject 
provided  In  section  3.  Plots  of  the  ellipses  are  provided  in  Append- .»  ' 
the  section  called  "Raw  (Non-Differential)  Loran-C  Error  Ellipses  ” 


Formulation  -  Differential  Lor an -C. 


If  we  place  a  differential  monitor  site  In  the  St  Lawrence  Seaway  -  s a 
at  the  Elsenhower  Locks  at  Massena,  N.Y.  -  we  can  considerably  reduce  the 
seasonal  component  of  the  TD  variations  observed  along  the  river.  The 
manner  In  which  the  variations  are  reduced  -  assuming  the  uniform 
propagatl on/double-range  difference  formulation  holds  true  -  Is  very 
straightforward.  From  equation  B-3  In  Appendix  B,  we  have  the  change  -n 
TD's  at  any  site  being  represented  by: 


TD2,p  -  TO]  *  ATDp  *  k  Av  DRDp 


If  we  call  the  differential  monitor  station  point  d,  we  have 
TD2,d  -  TOl  ,d  *  k  Av  DRD<i 
At  Tj  (say,  the  time  of  the  survey),  we  have 


TDi >p  *  kp  +  TDi ,d  (C-2) 

where  kp  Is  some  constant  value  determined  by  the  Loran-C  chain  numbers 
and  the  exact  location  of  the  observer  -  this  Is  all  determined  during  the 
survey.  At  T2,  we  have 


But, 

so  we  have 


™2,P 

-  TDi  ,p  +  k  AV 

•  kp  +  TDi  ,d 

TD1fd 

-  TD2,d  -  k  AV 

T°2.P 

“  kp  +  TD2>d  * 

*  kp  +  TD2  d 

DRDp 

♦  k  Av  DRDp 

DRDd 

k  av  DRDd  +  k  A v  DRDp 

♦  k  av  (DRDp  -  DRDd) 


In  utilizing  differential  Loran-C,  an  observer  at  point  p  computes  pos''Vor 
based  on: 

TD2,p  "  differential  correction 
where  differential  correction  *  TD2,d  -  TDl.d 


C-7 


Thus,  the  fix  Is  based  on 
TD2,p  -  TDa.d  +  TD^d 


«  kp  +  TD2,d  +  kAv  (DRDp  -  DRDd)  -  TD2,d  +  TDl,d 
*  Ckp  +  TDifd]  ♦  kA  v  (DRDp  -  DRDd) 


The  term  In  brackets,  of  course,  can  be  seen  from  equation  C-2  above  tc  <= 
simply  TDifP.  Thus,  we  have 

Corrected  TD2,p  ■  TDi  ,p  +  kA  v  (DRDp  -  DRDd) 


and  the  TD's  at  point  p,  with  differential  corrections  applied,  vary  from 
time  T-|  to  T2  In  the  manner  described  below: 


TOp,  differential  “  ^2,p,  differential  “  TDi  tp, differential 


*  kA  v  (DRDp  -  DRDd) 


Thus,  the  variations  are  of  the  same  form  as  described  In  Appendix  B,  we 
simply  have  to  use  what  we  might  call  the  "differential"  double  range 
difference.  We  pick  the  proposed  monitor  site  and  subtract  its  double  r»n g*- 
difference  from  all  other  sites  of  Interest. 


We  expect  that  the  move  to  differential  Loran-C  will  remove  a  large 
portion  of  the  TD  variations  throughout  the  St  Lawrence  Seaway.  Thus  we  ;?r 
no  longer  simply  assume  the  seasonal  components  are  so  much  larger  ~ha- 
shorter  term  components  that  thqy  serve  as  a  good  estimate  of  the  tc  * 
variation  and  we  must  come  up  with  estimates  for  the  short  term  var-.-v.- 
There  are  numerous  classical  approaches  to  this  problem  which  relate  v’e 
short  term  variation  to  distance  from  the  transmitting  stations.  As 
described  earlier  In  the  discussion  on  short  term  correlation  coefficie-  - 
these  treatments  assimie  atmospheric  noise  predominates  "transmitter  re  .  . 

Again,  we  argue  this  Is  not  an  appropriate  claim  for  the  St  Lawrenc®  S®-w 
vtoere  we  assume  we  are  not  "atmosperlc  noise  limited."  Thus,  for  tf*-  ■ 
cut,  we  will  choose  a  single  figure  to  represent  the  short  term  stan*?- 
deviation  -  for  both  TD’s  at  all  reaches.  From  harbor  monitor  expe*-:  "v 
an  estimate  of  20  nanoseconds  seems  most  appropriate. 

This  short  term  standard  deviation  must  be  modified  before  we  car  n. 
direct  application  to  the  differential  Loran-C  situation.  We  must  ~er  'e- 
that  the  differential  monitor  station  provides  Information  to  the  uv 
It  also  provides  noise.  Thus,  we  assume  the  short  term  noise  at  th*  tc 
Is  Independent  of  the  short  term  noise  the  user  Is  encountering  sc  th-k 


after  the  correction  has  been  applied,  the  user  noise  has  been  Increased: 


2  2  / 

0p,  differential  *  (0p  +  °d  )1/2 

We  assume  that  an  •  *  20  nsec  so  that  a p , dl f f  =  20  */r2~»  Thi  s  *s 

the  figure  we  will  use  for  our  predictions. 

We  further  assure  that  the  short  term  and  seasonal  variations  are 

independent  so  that 

0tota1  *  C*s2  ^  (»/*  )231/2  (C-3) 

-  {  CJS2  +  800)1/2 


where  crs  a  the  seasonal  component  standard  deviation  which  Is  estimated 
by  use  of  the  differential  double  range  difference  formulation. 


Additionally,  with  differential  Loran-C,  we  must  revise  our  approach  tc 
estimating  the  correlation  coefficient.  As  before,  we  will  assume  the 
seasonal  components  are  pure  sinusoids.  Depending  on  the  "sign"  of  the 
differential  double  range  differences  these  sinusoids  could  be  In-phase  or 
1800  out  of  phase.  Due  to  the  nature  of  the  geometry  in  the  St  Lawrence, 
the  differential  double  range  differences  are  always  of  the  same  sign. 

Thus,  we  will  again  use  pi  *  1  and  note  that  If  the  differential  double 
range  differences  were  of  the  opposite  sign,  we  would  have  used  pi  =  -1 . 

To  compute  the  total  correlation  coefficient,  we  re-write  equation  C-l : 

P  =  P1  °wl  axl  *  0  2  aw2  °x2 


Using  =1,  i>2  «0.5  as  before,  and  0*2  =  0x2  =  ZO*'?  *  this 
combines  with  equation  C-3  to  become: 


Jwl  qxl 
+  TOO)  (u 


+  400 


(C-4) 


Using  equation  C-4  and  the  re-computed  double  range  differences,  we  can 
construct  Table  C-2  which  lists  the  predicted  statistics  based  on  a 
differential  monitor  site  located  at  Elsenhower  Locks.  Plots  of  the 
resulting  error  ellipses  are  presented  In  Appendix  D. 
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Table  C-2 


Section  3  discusses  the  situation  to  be  attained  If  we  add  a 
differential  monitor  site  at  Wellesley  Island.  Table  C-3  lists  the 
statistics  for  appropriate  reaches  with  this  site  providing  differential 
corrections.  The  data  was  computed  using  the  same  formulation  used  for 
Elsenhower  Locks.  Plots  of  the  resulting  error  ellipses  are  also  provider 
In  Appendix  0. 
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DESCRIPTION  OF  SEAWAY  REACHES 
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APPENDIX  F 

EXAMPLES  OF  DATA  COLLECTION  PROBLEMS 


4 
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Comments  on  Examples 


Example 


Comment 


Figure  F-l 


Figure  F-2 


Figures  F-3 
and  F-4 


The  printout  is  from  the  Cote  Ste.  Catherine  site 
Ad  Hoc  equipment.  The  column  labeled  "GRI"  provides 
a  measure  of  the  receiver  oscillator  stability.  It 
should  be  very  steady  at  some  value  In  the  range 
99600.00  +/-  0.20.  The  significant  variation  In 
this  measure  of  oscillator  performance  is  reflected 
in  the  unstable  time  difference  readings.  This 
malfunction  was  the  major  cause  of  lost  data  at 
Cote  Ste.  Catherine. 

The  printout  is  from  the  Iroquois  Lock  site  Ad  Hoc 
equipment.  (Note  the  stable  oscillator  readings.) 
Receiver  estimated  slgnal-to-nolse  ratios  for  the  M, 
W,  X,  Y,  and  Z  signals  are  listed  In  the  last  five 
columns,  respectively.  Notice  how  the  SNR's  periodi¬ 
cally  drop  substantially  and  then  return  to  normal. 
Related  fluctuations  In  the  time  differences  (columns 
2  through  5)  can  be  seen.  At  present,  the  cause  of 
these  fluctuations  Is  unknown. 

The  printouts  are  from  the  Prescott  Buoy  Auditing 
equipment  sets.  The  data  of  figure  F-3  was  obtained 
on  Julian  Day  116.  The  column  labeled  "DB" 
represents  the  deviation  of  the  9960-X  average  from 
the  nominal  value  of  28130.93  In  tens  of  nanoseconds. 
The  column  labeled  "SDB"  represents  the  standard 
deviation  of  the  time  differences,  also  in  tens  of 
nanoseconds,  ever  the  period  for  which  the  average 
value  In  column  "DB"  was  computed.  The  fluctuations 
in  the  averages  are  extraordinarily  high  as  are  the 
standard  deviations.  At  this  point.  It  could  have 
been  concluded  there  was  a  receiver  problem. 

Figure  F-4  shows  the  data  from  the  same  site  2 
weeks  later.  The  receiver  Is  still  abnormal.  The 
two  figures,  taken  together,  illustrate  the  major 
characteristic  of  the  data  bases  from  the  high 
density  data  collection  sites:  2  weeks  after  the 
problem  was  first  evident,  It  still  has  not  been 
detected  and  corrected.  (The  receiver  was  eventually 
changed  on  Julian  Day  159  -  more  than  6  weeks  after 
the  problems  were  evident.) 
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Comments  on  Examples  (con't) 


Example  Comment 


Figure  F-5  The  printout  is  from  the  Beauharnois  site  Ad  Hoc 

equipment  set.  Two  problems  are  illustrated. 

Notice  the  third  column*  ‘"he  collect  9960-X  time 
difference  reading  is  about  27579.  usee.  The 
receiver  is  tracking  on  the  wrong  cycle  (thus  the 
reading  is  about  10  usee  low).  Later  it  is  tracking 

on  the  correct  cycle  (the  heading  "INTERNAV . 

....300  BAUD"  means  an  operator  intervened).  Still 
later  it  is  again  tracking  10  usee  low. 

The  second  problem  is  "when  did  this  all  occur?" 

In  the  heading,  the  notation  "00/00/00  00:00"  means 

"it  is  now  midnight  on  the  0th  day  of  the  0th  month 
of  the  0th  year."  (Should  a  star  have  been  visible 
in  the  midnight  clear  or  did  that  happen  a  year 
later?)  After  seven  more  readings  are  taken,  and, 
after  data  collection  was  interrupted  for  some 
indeterminate  amount  of  time,  we  are  again  led  to 
believe  it  is  still  midnight  on  the  same  day. 

These  unknown  times  of  observation  make  the  data 
unusable  for  Differential  Loran-C  analysis  purposes. 


Figure  F-5  The  data  here  was  also  obtained  at  Beauharnois.  The 

two  sets  of  data  show  a  substantial  offset.  This  is 
not  a  "seasonal  effect"  nor  is  it  the  result  of  an 
equipment  problem.  Careful  investigation  by  the 
authors  revealed  the  antenna  was  moved  to  a  new 
location  during  this  period  (Loran-C  works!) 
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